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This  final  report  was  prepared  by  the  Center  for  Regional 
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Kirtland  AFB,  New  Mexico. 

This  report  consists  of  three  volumes:  Volume  I - 

Evaluation  and  Recommendations,  Volume  II  - Biological 
Inventory  1974/1975,  Volume  III  - Environmental  Planning 
System.  Volume  III  was  presented  as  a master's  thesis  to 
San  Diego  State  University  in  partial  fulfillment  of  the 
requirements  for  the  degree,  Master  of  Science  in  Biology, 
by  Richard  M.  Reilly. 

This  report  has  been  reviewed  by  the  Information  Office 
(01)  and  is  releasable  to  the  National  Technical  Information 
Service  (NTIS) . At  NTIS  it  will  be  available  to  the  general 
public,  including  foreign  nations. 
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SECTION  I 
INTRODUCTION 

This  report  contains  a description  and  documentation  of  the  computer- 
based  ecological  evaluation  tools  developed  for  Vandenberg  Air  Force  Base 
(VAFB).  It  is  presented  as  a companion  to  Volume  I (Evaluation  and 
Recommendations)  and  Volume  II  (Biological  Inventory),  (References  1,2). 

This  volume  is  primarily  intended  for  use  by  personnel  involved  in  environ- 
mental planning. 

The  ecological  inventory  developed  for  VAFB  provided  the  basic  data, 
in  computer-compatible  form,  for  the  Environmental  Planning  System  (EPS) 
presented  in  this  volume.  The  EPS  is  designed  to  be  a flexible  tool  usable 
by  base  personnel  having  a minimum  of  computer  experience.  The  EPS  will 
calculate  and  display  the  area  of  each  major  vegetation  type,  or  other 
ecologically  sensitive  unit,  that  is  likely  to  be  altered  by  any  proposed 
\ development  configuration.  It  will  assist  in  the  identification  of  optimal 

locations  for  development  or  operational  actions,  either  for  the  Space 
Transportation  System  (STS)  or  other  elements  of  the  VAFB  mission,  on  the 
basis  of  their  relative  ecological  suitability. 

i This  report  is  functionally  organized  in  the  interests  of  the  user. 

Sections  II  and  III  briefly  describe  the  EPS  and  how  to  use  it.  Section 
IV  is  a series  of  actual  case  studies  using  the  EPS.  Section  V is  a 
comparison  of  methods  for  determining  areas  of  vegetation.  The  data 
presented  in  this  section  provide  a comparison  of  the  accuracy  of  the 
manual  method  versus  the  automated  method  for  determining  areas  of  vegetation. 
Section  VI  is  a detailed  Users  Manual  which  describes  how  to  use  the 
computer  programs  in  the  EPS.  Section  VII  is  a description  of  the  computer- 
**  ized  data  base  in  the  EPS  along  with  instructions  for  its  use  and  updating 

procedures.  The  appendix  contains  complete  program  listings  for  the 
computer  programs,  subroutines  used  in  the  case  studies,  and  the  portion 
of  the  data  base  containing  the  irregular  outline  used  by  the  GRID  computer 
• program  to  map  VAFB. 


SECTION  II 


DESCRIPTION  OF  THE  ENVIRONMENTAL  PLANNING  SYSTEM 

The  EPS  used  for  VAFS  was  based  upon  a study  done  for  the  coastal  plain 
of  San  Diego  County  (Reference  3).  Use  of  this  type  EPS  seemed  to  allow  for 
the  concurrent  experimental  and  observational  design  of  the  ecological 
assessment  of  VAFB. 

The  GRID  computer  mapping  program  was  selected  to  display  the  computer- 
ized data  base  developed  from  the  ecological  inventory.  The  GRID  program 
was  developed  at  the  Harvard  Laboratory  for  Computer  Graphics  and  Spatial 
Analysis.  It  provides  a highly  efficient  means  for  graphic  display  of  large 
quantities  of  information  collected  on  the  basis  of  a rectangular  coordinate 
grid.  The  GRID  program  is  written  in  FORTRAN  IV.  It  is  operational  on  the 
IBM  360  computer  at  San  Diego  State  University  and  the  Burroughs  3500 
computer  at  VAFB. 

The  GRID  computer  mapping  program  is  an  automated  alternative  to 
drawing  and  planimetering  overlays  for  mapping  studies.  The  GRID  program 
produces  a composite  map  output  which  selectively  displays  from  two  to  ten 
classes  of  data.  GRID  calculates  the  frequency  of  each  class  of  data  (the 
number  of  grid  cells),  the  percent  total  by  classes,  and  displays  the 
results  with  a histogram.  GRID  has  two  distinct  advantages  over  manual 
methods : 

1.  Combinations  of  overlays  for  displaying  alternatives  can  be 
produced  almost  without  limit. 

2.  The  cost  is  minimal  once  the  data  base  and  the  computer 
graphics  program  are  established. 

To  apply  the  GRID  program  to  VAFB,  the  area  encompassing  the  base  was 
divided  into  rectangular  grids.  Two  grid  scales  of  analysis  were  used.  A 
square  grid  cell  1000  feet  on  a side  was  adopted  for  standard  processing 
by  the  GRID  program.  This  cell  size  was  chosen  because  it  was  the  same  as 
the  1000-foot  grid  of  the  California  Coordinate  System  used  in  the  Base 
Master  Plan  maps  (Reference  4).  This  provided  a raster  of  4646  cells 
(22.96  acres  per  cell)  covering  the  entire  base.  For  data  requiring  more 
detailed  resolution,  such  as  vegetation  cover,  the  basic  grid  cell  was 
subdivided  into  nine  square  subunits  of  2.55  acres  each.  An  interpretive 
subroutine  in  the  GRID  program  then  provides  a decision  for  the  assignment 
of  each  23-acre  cell  to  a category  determined  by  the  values  of  the  nine 
subunits. 

The  EPS  also  contains  a Search/Count  program.  This  program  is  run 
independently  of  GRID  with  the  computerized  data  base.  It  will  locate  and 
print  out  data  base  variables  for  selected  grid  cells.  By  entering  the 
alphanumeric  coordinates  of  selected  grid  cells,  the  program  will  print 
out  the  data  base  variables  for  each  grid  cell  and  count  the  frequency  of 
each  subvariable.  The  Search/Count  program  can  be  used,  for  example,  to 
compute  the  areas  of  various  vegetation  types  that  would  be  affected  by 
alternate  placements  of  a runway-launch  complex. 
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SECTION  III 

HOW  TO  USE  THE  ENVIRONMENTAL  PLANNING  SYSTEM 

To  effectively  operate  the  EPS,  as  a flexible  tool  for  environmental 
analysis  and  planning,  the  user  must  be  familiar  with  the  components  of  the 
system  and  how  they  operate.  Section  II  gave  an  overview  of  the  EPS.  A 
detailed  Users  Manual  for  operating  GRID  and  the  Search/Count  program,  as 
adapted  for  VAFB,  is  presented  in  Section  VI.  Many  examples  for  use  of 
different  program  options  were  incorporated  into  the  Users  Manual  to  aid  in 
understanding.  A complete  description  of  the  computerized  data  base 
developed  for  VAFB,  and  instructions  for  using  it,  are  given  in  Section  VII. 
The  combination  of  the  Users  Manual  and  the  data  base  description  and 
instructions  will  provide  the  user  with  the  tools  to  operate  the  system. 

A series  of  case  studies  using  the  EPS  are  presented  in  Section  IV. 
Several  of  the  case  studies  are  fully  documented  in  the  appendix  with 
program  listings.  In  addition,  the  appendix  contains  a complete  program 
listing  for  GRID  and  the  Search/Count  program  as  adapted  to  the  Burroughs 
3500  computer.  To  gain  initial  experience  operating  the  EPS,  the  user 
should  practice  by  duplicating  the  computer  output  of  one  or  a number 
of  the  case  studies. 

i 

The  user  now  has  the  tools,  the  data  base,  and  the  instructions  for  a 
computerized  environmental  planning  system.  The  challenge  is  now  to  imple- 
ment the  EPS  for  operational  environmental  analysis  and  planning. 
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SECTION  IV 


CASE  STUDIES  - APPLICATION  OF  THE  ENVIRONMENTAL  PLANNING  SYSTEM 
TO  VANDENBERG  AIR  FORCE  BASE 

The  following  are  a series  of  case  studies  to  illustrate  different 
uses  of  the  EPS. 

1.  DISTRIBUTION  OF  ELEVATION  CLASSES  FOR  VAFB 

Elevation  was  coded  in  200-foot  intervals  for  each  grid  cell.  The 
highest  elevation  in  the  grid  cell  determined  the  code  assigned  to  that  cell. 
Using  this  system,  10  different  codes  were  used  to  describe  the  elevation 
from  sea  level  to  the  highest  point  of  VAFB.  All  areas  1800  feet  and  above 
were  combined  into  one  elevation  class.  The  highest  elevation  on  the  base 
is  2170  feet. 

The  GRID  program  was  used  to  display  the  distribution  of  the  elevation 
classes  on  VAFB.  Table  1 is  a summary  of  the  elevation  classes  based  on 
GRID  display  output.  The  same  information  could  easily  be  obtained  by 
planimetering  the  appropriate  contours  on  a topographic  map.  This  case  is 
too  simple  to  display  the  full  power  of  the  EPS,  but  it  is  desirable  to 
include  for  just  that  reason. 

The  elevation  data,  along  with  other  data  variables,  can  be  used  to 
study  general  hydrologic  features  on  the  base.  A GRID  map  could  be  printed 
to  study  lowland  drainage  areas  by  printing  out  all  the  areas  with  an 
elevation  of  less  than  600  feet.  The  areas  shown  might  also  indicate  areas 
dominated  by  coastal  fog  if  600  feet  represents  the  mean  low  ceiling  (see 
Volume  I,  Figure  3.2.3).  Similarly,  areas  below  1200  feet  may  be  consistently 
below  the  inversion  layer  and  thus  subject  to  potential  effects  of  the 
exhaust  cloud  from  a rocket  launch  (see  Volume  I,  Figure  3.2.4). 

2.  AREAS  SENSITIVE  TO  EROSION  ON  VANDENBERG  AFB 

This  study  uses  the  GRID  program  and  data  base  variables  to  define 
areas  of  high  environmental  sensitivity  to  erosion.  The  soil  and  vegetation 
data  base  variables  were  used  in  the  study.  They  are  described  in  detail  in 
Section  VIII . 

The  erosion  potential  of  the  soil  was  analyzed  using  the  soil  series/ 
phase  or  land  type  with  its  slope  characteristic.  Using  this  procedure, 

62  soil  codes  in  the  data  base  were  defined  to  have  a high  erosion  potential. 
These  included  the  following  groups:  severely  eroded  soils,  soils  with 

steep  slopes  (greater  than  15  percent),  and  unsuitable  land  types  (based  on 
engineering  restraints).  The  term  "land  type"  is  used  to  describe  areas 
where  the  soil  materials  are  too  rocky,  shallow,  severely  eroded,  sandy,  or 
wet  to  be  classified  as  distinct  soil  types.  The  category  "unsuitable  land 
types"  includes  land  types  on  VAFB  judged  to  be  unsuitable  for  construction 
purposes  based  on  engineering  restraints.  Examples  of  unsuitable  land 
types  are:  coastal  beaches,  sedimentary  rock  land,  gullied  land,  etc. 

A GRID  map  (Figure  1)  was  produced  to  display  the  erosion  potential 
of  the  soil  in  5 levels  based  on  the  above  criteria.  They  are: 
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TABLE  1.  SUMMARY  OF  ELEVATION  CLASSES  FOR  VANDENBERG  AIR  FORCE  BASE,  BASED 
ON  GRID  DISPLAY  OUTPUT. 


Elevation  Classes 
(Feet) 

No.  of  Cells 

*Estimated 

Acres 

Percent  of 
Total  Cells 

0-  200 

844 

19378 

18.2 

200-  400 

1339 

30743 

28.8 

400-  600 

963 

22110 

20.7 

600-  800 

528 

12122 

11.4 

800-1000 

406 

9322 

8.7 

1000-1200 

248 

5694 

5.3 

1200-1400 

131 

3008 

2.8 

1400-1600 

88 

2020 

1.9 

1600-1800 

50 

1148 

1.1 

1800  and  over 

49 

1125 

1.1 

TOTALS 

4646 

106670 

100.0 

*Estimated  acres  = no.  of  cells  x 22.96  acres  per  cell. 
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Level  1 = Background  (.)  e.g.,  all  soils  not  classed  in  levels  2 to  5. 
Level  2 = Severely  eroded  soil  (+) 

Level  3 = Soil  with  steep  slopes  (X) 

Level  4 = Severely  eroded  soil  with  steep  slopes  (e) 

Level  5 = Unsuitable  land  types  (■) 

The  Flexin  subroutine  to  produce  the  GRID  map  in  Figure  1 is  listed  in 

Appendix  C.  The  DATA  statements  in  Subroutine  Flexin  list  the  individual 

soil  codes  in  each  of  the  5 levels.  Table  2 is  a summary  of  the  5 levels 
based  on  the  GRID  map  output. 

The  GRID  program  has  the  advantage  of  producing  computer  maps  displaying 
different  combinations  of  variables  almost  without  limit  at  minimal  cost. 

This  is  demonstrated  by  another  GRID  map  (Figure  2).  This  GRID  output  used 
the  soil  erosion  potential  data  (described  above)  combined  with  unstable 
vegetation.  Unstable  vegetation  was  defined  as  annual  grassland  in  this 
study. 

A GRID  map  was  produced  to  display  areas  of  potentially  high  soil 
erosion  and  unstable  vegetation  in  10  levels.  They  are: 

Level  1 = Background  (.)  e.g.,  all  soils  and  vegetation  types  not 
included  in  levels  2 to  10. 

Level  2 = Severely  eroded  soil  (,) 

Level  3 = Soil  with  steep  slopes  (■) 

Level  4 = Severely  eroded  soil  with  steep  slopes  (+) 

Level  5 = Unsuitable  land  types  (X) 

Level  6 = Annual  grassland  (0) 

Level  7 = Annual  grassland  and  severely  eroded  soil  (e) 

Level  8 = Annual  grassland  and  soil  with  steep  slopes  (e) 

Level  9 = Annual  grassland  and  severely  eroded  soil  with  steep  slopes 

(8) 

Level  10  = Annual  grassland  and  unsuitable  land  types  (•) 

The  Flexin  subroutine  to  produce  the  GRID  map  in  Figure  2 is  listed  in 
Appendix  D.  Table  3 is  a summary  of  the  10  levels  based  on  GRID  output. 

Other  combinations  of  eroded  soil,  steep  slopes,  and  unstable  vegetation 
can  be  mapped  using  GRID  and  the  data  base  variables.  For  example,  the 
exposure  variable  contains  subvariable  codes  for  multiple  exposures  in  a 
cell  due  to  drainage  areas  or  ridge  lines.  These  subvariable  codes  would 
indicate  sloping  land  with  potentially  eroded  soil.  A GRID  map  could  be 
produced  to  display  the  variables  for  soil,  vegetation,  and  exposure  in  up 
to  10  different  levels  depending  upon  the  design  of  Subroutine  Flexin. 

3.  USE  OF  GRID  TO  DISPLAY  AREAS  OF  BIOTIC  SENSITIVITY 

The  GRID  program  and  data  base  developed  for  VAFB  can  be  used  to  locate 
and  print  out  areas  likely  to  contain  piant  or  animal  species  o*  substantial 
importance  that  must  be  considered  in  an  Environmental  Impact  Statement. 

As  a result  of  the  ecological  assessment  of  VAFB  three  vegetation  areas  on 
the  base  were  identified  to  be  of  prime  ecological  significance  and  were 
cited  as  a potential  environmental  problem  deservinq  special  consideration 
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TABLE  2.  SUMMARY  OF  GRID  LEVELS  DISPLAYING  AREAS  WITH  SOILS 
OF  HIGH  EROSION  POTENTIAL. 


Level 

No.  of  Cells 

Estimated 

Acres 

Percent  of 
Total  Cells 

1 

2648 

60798 

57.0 

2 

25 

574 

0.5 

3 

1152 

24.8 

4 

89 

2043 

1.9 

5 

732 

16807 

15.8 

TOTALS 

4646 

106672 

100.0 

TABLE  3. 
EROSION 

, SUMMARY  OF  GRID  LEVELS 
POTENTIAL  BASED  ON  SOILS 

DISPLAYING  AREAS 
AND  VEGETATION. 

OF  HIGH 

Level 

No.  of  Cells 

Estimated 

Acres 

Percent  of 
Total  Cells 

1 

1674 

38435 

36.0 

2 

8 

184 

0.2 

3 

687 

15774 

14.8 

4 

16 

367 

0.3 

5 

589 

13523 

12.7 

6 

974 

22363 

21.0 

7 

17 

390 

0.4 

8 

465 

10676 

10.0 

9 

73 

1676 

1.6 

10 

143 

3283 

3.1 

TOTALS 

4646 

106671 

100.1 
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for  the  STS  program  (Reference  1).  The  areas  are:  Tanbark  Oak  forest. 

Bishop  Pine  forest,  and  stabilized  sand  dunes. 

Each  of  the  areas  above  is  a vegetation  type  and  is  coded  as  a 
subvariable  in  the  computerized  data  base.  A GRID  map,  Figure  3,  was 
produced  to  display  areas  of  prime  ecological  significance  in  7 levels. 

They  are: 

Level  1 = Bishop  Pine  forest  (X) 

Level  2 = Tanbark  Oak  forest  (0) 

Level  3 = Stabilized  sand  dunes  (•) 

Level  4 = Background  (.)  e.g.,  all  vegetation  types  not  included  in 
levels  1 to  3 or  levels  5 to  7 
Level  5 = Bishop  Pine  forest-partial  occurrence  (+) 

Level  6 = Tanbark  Oak  forest-partial  occurrence  (§) 

Level  7 = Stabilized  sand  dunes-partial  occurrence  (0) 

The  Flexin  subroutine  to  produce  the  GRID  map  in  Figure  3 is  listed  in 
Appendix  E.  The  vegetation  subvariables  were  coded  at  the  2.55-acre  subcell 
level.  Subroutine  Flexin  was  written  to  direct  the  computer  to  read  each 
of  the  9 subcells  for  a grid  cell.  A routine  then  calculated  the  vegetation 
type  of  the  grid  cell  using  the  vegetation  type  occurring  in  a majority  of 
the  subcells.  If  the  vegetation  type  was  one  of  the  key  types,  it  was 
mapped  as  level  1,  2,  or  3 for  the  grid  cell.  If  a key  vegetation  type  was 
present,  but  did  not  occur  in  the  majority  of  the  subcells,  it  was  mapped 
as  level  5,  6,  or  7.  All  other  vegetation  types  were  mapped  as  background. 
Table  4 is  a summary  of  the  7 levels  based  on  the  GRID  output. 


The  GRID  map  displaying  areas  of  prime  ecological  significance  can 
1 also  be  interpreted  to  display  the  distribution  of  plant  species  which  are 

characteristic  of  the  vegetation  type.  For  example,  6 rare  plant 
> species  are  restricted  to  the  stabilized  sand  dunes  (Reference  1). 

Therefore,  the  6 rare  plant  species  have  a high  probability  of  occurring 
within  grid  cells  mapped  as  stabilized  sand  dunes. 


i 


Similarly,  the  probable  distribution  of  some  animals  can  be  mapped 
based  on  their  habitat  preference.  The  California  Legless  Lizard  (Anniel la 
pulchra)  is  a species  regulated  by  the  California  Fish  and  Game  Commission. 
It  occurs  on  VAFB  and  is  considered  relatively  sensitive  due  to  its  habitat 
specificity  (Reference  1). 

A.  pulchra  is  a burrowing  reptile  and  is  relatively  restricted  to  soils 
it  can  penetrate.  Based  on  the  field  studies  done  at  VAFB  during  the 
ecological  assessment  and  Miller's  research  (Reference  5),  suitable  habitat 
for  A.  pulchra  was  described  using  subvariables  in  the  computerized  data 
base.  Suitable  habitat  for  A.  pulchra  included:  a suitable  soil  type 

(sandy  or  loamy  soils),  a preferred  vegetation  type,  or  a combination  of 
both. 


A GRID  map  (Figure  4)  was  produced  to  display  areas  of  suitable  habitat 
in  10  levels.  They  are: 
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LEGEND 

BISHOP  PINE  FOREST  X 

TANBARK  OAK  FOREST  0 

STABILIZED  SAND  DUNES  I 

BACKGROUND 

BISHOP  PINE  FOREST- 
PARTIAL  OCCURRENCE  + 

TANBARK  OAK  FOREST- 
PARTIAL  OCCURRENCE  © 

STABILIZED  SAND  DUNES- 
PARTIAL  OCCURRENCE  0 


FIGURE  3.  AREAS  OF  PRIME 
ECOLOGICAL  SIGN  I FCANCE  ON  VAFB. 
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TABLE  4.  SUMMARY  OF  GRID  LEVELS  DISPLAYING  AREAS  OF  PRIME 
ECOLOGICAL  SIGNIFICANCE  OF  VAFB. 


Level 

No.  of  Cells 

Estimated 

Acres 

Percent  of 
Total  Cells 

1 

20 

459 

0.4 

2 

2 

46 

0.0 

3 

411 

9437 

8.8 

4 

4032 

92575 

86.8 

5 

45 

1033 

1.0 

6 

9 

2G7 

0.2 

7 

127 

2916 

2.7 

TOTALS 

4646 

106673 

99.9 

TABLE  5 
HABITAT 
ON  VAFB 

SUMMARY  OF  GRID 
FOR  THE  CALIFORNIA 

LEVELS  DISPLAYING  AREAS  OF  SUITABLE 
LEGLESS  LIZARD  (ANNIELLA  PULCHRA) 

Level 

No.  of  Cells 

Estimated 

Acres 

Percent  of 
Total  Cells 

1 

14 

321 

0.3 

2 

15 

344 

0.3 

3 

5 

115 

0.1 

4 

39 

895 

0.8 

5 

1360 

31226 

29.3 

6 

68 

1561 

1.5 

7 

374 

8587 

8.0 

8 

4 

92 

0.1 

9 

146 

3352 

3.1 

10 

2621 

60178 

56.4 

ITOTALS 


4646 


106671 


99.9 
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FIGURE  4.  AREAS  OF  SUITABLE 
HABITAT  FOR  THE  CALIFORNIA 
LEGLESS  LIZARD  ON  VAFB. 


I J** 
1 3- 
133 

; it 
131 


080 C 

.00 , 

,600b. . .too o 

€66080 0 

.ee«0««M 

. 66H008UCLUG 0 UOCO 

tibtshbBab. . iac  .uuc uo.c 

883«8!Si- “&£?: : : “ ii8o.a..o<*x1®oSiooo 

: 8SI!!!!!MSSi?n8SS&: : : ?V8S8§o 

b8LeHaeL0««Mt«oouocnouaa.  .aoaaaa oooo 

660000800000000. gum  . . uuuaaauBanaa .on. 

. uc'uoooooocGaacu . ctaoo  .o  .ooo 

.0«feeO«:10086fc6*tBB0..CuC.uUCCCBCaB0UCO010OUOaO 
fc«8088. BBBBPfl. 80,890000.0000.00000* . OBBi. . UOOOBOOOOU. 
brtBB8BinSBC:ietB.«0«BBmOoOU.B.OOOC.O.O.  . .0.000000 CS0UO 
m"bbS6«bbmcbbmnnh^70oi:o  . • . . o. . ooooooooooo«JOCOBoon 

,680C808B8800008000000OCO 000000  . -00000.  omOiJOBBOO 

.b0BBe08B«88C 80. 000800000.  . O.XOO0.O.  . . UOOBOOUIXJO0C0 000 
.ee800B088eUBPf08800U008OliU)U.  *00000000.  q000q(M»)O00000J 
.6b  . BB 6808 8888 8o88Bi;B(jB0(JOBO BO  - . . OOBBOO.  000000000000080 
.e€BtS8B0«ljO(;88OBOOOOf'OBBO*UOO<WOOOOBOOOOBBtOBBaaBBBO8U. 

. 68888  . . BBBBOOcOOUBBOOOOOC  .1)000000000000.  BBOUXBB.  . BBB  xOBD 
. 6BBB6 . . 88.  oBocBoocounccxjcMjLioooouuoooBBB. ■aooaaaBoat.  .oao 
8000000 . 00 . . . . CCOOOOOOOOUO. 000.0000000000000 .0000 . OOQUOBB 

(360086 COO.  00000 00  . . .0000  . . . 0000000)000000 

66MHB0  \ IJO  . 00  ..0.  .0000000.  .BB 

. 6b 08B888 . . uouua .OB.  .a* . .c.c ooooo. . 

6680868 . .0.0. ...CO.  .0.0.  ■ iUCOOOOOOOOOO. OC.GOO  .0 0000. . 

0680008 0.  ...  0.  . b.OU.OU.OOGO.  . . .000.  .00  .QUO ' 

00 *».  .0000 OCU.  .0000  . .0 

008080881' UQOC.  . . 00 . . 000 . . UUUGQG  .000000 . 0 . .... 

HB0BM0O0:  : . 1 . .o!  . .000 OOO.OBCOOOO0OOBOOOBOOOO 

6000000 uO.  00000.  . .000.000000000000*000.0..  . . 0 

0000001  . ........  0. OCCOOOO. . .000 . ..X  ... .0O.OO.X. .0. . .OB 

. uc.oooo.  .ouooc. *. . -u... ...  .oo 

00080 o.ooo. . ooooooc *...  .00X..XJ..0 

8U880 u. aocm.u.uouoouoo oouo.oxx.o. . 

C0H8(  01  o couoooo 0000.  .000.0  . 

060000  ....0 O..UO..U.  QOO  .0.  . ^ .00000 .0 


■ -MfeHS-.M,:  v!\ : :S8ae*6":  fVSi«8i«?88?,i 

♦ BCCOO.O...  *”  r'"  n 

. .occ. 


. . .OOOOOOC. XU. 0000000.0 
. .OOC.B.G..C.. . OB. 

CC  .OOO.  . . .00..  . 

BOOlio OCO.OO.O.  XOOO.B. 

c CLLi o.  .000*  .000.  .OU 

cocoa oooo. . .oooo.#o 

OCCOOOO OOOO.  .BCUO  .00 

....o.oocc 0 cgcc.ocou.oo 

. . 0000000 0 O00O  .0 . 170 

. uooooo o. . .oxxo.oqo 

,U.U.O c CO. C.  . .00000 

8.  .0.0 0 00.0.  . .000 

68..0.00 ...  .0  . .00  ...  .000 

. 6.  . . . JCC.OOO 0.  . . C.OOCO..OC 

.6 OOCUOOCiU.U.  . '.  .COOUB.O 

OOOOM.UCUUU.  t.CGOOOO 

a.  j ooconroooo  .on 

. 0 .uO ocuuoc.oo 

,8.0000 CO 0 

*0.0000 

. ecccaucG 

,80.00. .000 

bo.  .no. . .ocoo 

.BCCLOOCO.  OCCUC.  . . . 

•occcoooca. .on 

6CCCC0GU.. .LOU.U 

. BOCCCuOCC . . CC.CUOOO.  . . 


I?  1 

uo 


10? 
I Oo 
UN 

l j<t 

u. 
102 
101 
I Oo 
09V 
i:98 

jw 

oo< 

ovb 

oSis 


.tcooccnoaoocjcco.  c oooooo. 

B0OOCC.OOOOOCLICCCOO.O. . 

. .BCCC ccucccoccccc cococo 
.socaccc.  norocc  (.).  oooo  o 

. 0OCOOCCOOOCOOLCII.  COOOOOU 
. •BCU.CCCC.OCGOC lu. CCOOOCCL . 

,tooa  pooacoca.  c co. cooouonro.o 

. 80 . . . JUOGUO. OOOCtO . . OUUCUULU00 

08 C.  CCCUO.  .OOO.  . .0 

• CCL 0 


.0. 


.'ololo. 

. .0000. 


A A A \ A 6 A ft  A t.AAAM. 


14 


«WAr.  • • uu 

lot-*-  Cl  Cl  C iufu  Ci.ct 


I 


I 


Level 

Level 

Level 

Level 

Level 

Level 

Level 

Level 

Level 

Level 


1 = Coastal  strand  ( , ) 

2 = Stabilized  sand  dunes  (■) 

3 = Coastal  bluff  (+) 

4 = Oak  woodland  (X) 

5 = Suitable  soil  associations  (0) 

6 = Coastal  strand/suitable  soil  overlap  (9) 

7 = Stabilized  sand  dunes/suitable  soil  overlap  (9) 

8 = Coastal  bluff/suitable  soil  overlap  (8) 

9 = Oak  woodland/suitable  soil  overlap  (•) 

10  = Background  (.)  e.g.,  areas  of  unsuitable  habitat. 


The  Flexin  subroutine  to  produce  the  GRID  map  in  Figure  4 is  listed  in 
Appendix  F.  The  data  statement  in  Subroutine  Flexin  contains  22  soil  codes 
which  were  designated  as  suitable  soil  associations.  Table  5 is  a summary 
of  the  10  levels  used  to  display  suitable  habitat  for  A.  pul chra  based  on 
the  GRID  computer  output.  The  levels  of  suitable  habitat  were  not  ranked 
in  this  study,  but  with  expert  consultation  this  could  be  done. 


4.  EFFECTS  ON  VEGETATION  TYPES  DUE  TO  ALTERNATE  PLACEMENTS  OF  A RUNWAY- 
LAUNCH  COMPLEX  USING  THE  SEARCH/COUNT  PROGRAM 

The  Search/Count  program  developed  for  the  EPS  is  useful  for  making 
detailed  studies  within  the  general  study  area  of  VAFB.  Given  the  alpha- 
numeric coordinates  of  selected  grid  cells,  the  Search/Count  program  will: 
direct  the  computer  to  locate  the  grid  cells  in  the  data  base;  frequency 
count  the  subvariables  of  the  selected  grid  cells  (and  convert  them  to 
acres  if  requested);  print  out  the  data  base  file  for  each  selected  grid 
cell;  and  print  out  the  frequency  counts  of  the  subvariables  for  the  group 
of  selected  grid  cells.  Instructions  for  using  the  Search/Count  program 
are  given  in  Section  VI.  A listing  of  the  Search/Count  program  is  given  in 
Appendix  G. 

In  this  case  study,  the  Search/Count  program  was  used  to  compute 
areas  of  various  vegetation  types  that  would  be  affected  by  alternate  place- 
ments of  a runway-launch  complex  to  be  used  for  the  proposed  Space  Transport 
System  (STS).  The  dimensions  of  the  runway-1 aunch  complex  are  given  in  the 
DOD  STS  Facility  Development  Specification  (Reference  6).  The  grid  cell 
layout  requires  that  analysis  of  anything  other  than  exact  North-South  or 
East-West  layout  include  adjacent  areas. 

Figure  5 is  a diagram  showing  the  relative  positions  of  the  grid  cells 
affected  if  the  runway-launch  complex  was  located  on  the  present  runway 
facility  at  VAFB  with  a new  runway  extension  and  additional  construction  of 
facilities.  The  darkened  grid  cells  represent  the  area  affected  by  the 
runway  extension  only  (34  grid  cells).  The  light  area  in  the  diagram 
represents  the  existing  runway  with  existing  and  proposed  facilities  (43 
grid  cells). 

Table  6 is  a summary  of  the  vegetation  types  affected,  at  the  2.55- 
acre  subcell  resolution,  for  the  proposed  runway  extension  only  (dark  area 
of  Figure  5)  based  on  the  Search/Count  program  output.  Table  7 is  a 
similar  summary  of  the  vegetation  types  affected  in  the  light  area  of 
Figure  5. 
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Shaded  grid  cells  = Area  affected  by  a new  runway 

extension 


Unshaded  grid  cells  = Area  affected  by  existing 

runway  with  facilities  and 
proposed  facilities 


FIGURE  5.  RELATIVE  POSITIONS  OF  AFFECTED  GRID  CELLS  DY  A RUNWAY-LAUNCH 
COMPLEX  LOCATED  ON  THE  PRESENT  RUNWAY  FACILITY  AT  V AND EN BERG  AFB. 


TABLE  6.  SUMMARY  OF  THE  VEGETATION  TYPES  AFFECTED  BY  THE 
PROPOSED  RUNWAY  EXTENSION  (34  GRID  CELLS)  BASED  ON  THE 
SEARCH/COUNT  PROGRAM. 


Vegetation  Type 

No.  of 
Subcells 

Estimated 

Acres 

Grassland-annual 

43 

122 

Chaparral 

130 

332 

Ruderal  vegetation 

72 

184 

Coastal  sage  scrub- 
normal  phase 

8 

20 

Planted  trees 

5 

13 

Chaparral-sparse  phase 

34 

87 

Coastal  sage  scrub- 
stabilized  dune  phase 

1 

3 

Man-made  facilities 

8 

20 

TOTALS 

306 

781 

TABLE  7.  SUMMARY  OF  THE  VEGETATION  TYPES  AFFECTED  BY  THE 
EXISTING  AND  PROPOSED  RUNWAY  FACILITIES  (43  GRID  CELLS)  BASED 
ON  THE  SEARCH/COUNT  PROGRAM. 

Vegetation  Type 

No.  of 
Subcells 

Estimated 

Acres 

Chaparral 

54 

138 

Ruderal  vegetation 

256 

653 

Coastal  sage  scrub- 
normal  phase 

2 

5 

Man-made  facilities 

64 

163 

Grassland-annual 

11 

28 

TOTALS 

387 

987 

Two  alternate  configurations  of  the  runway-launch  complex  are  tested  to 
demonstrate  the  flexibility  of  using  the  Search/Count  program.  Configura- 
tion I shows  the  effect  on  the  various  vegetation  types  if  the  runway-launch 
complex  was  constructed  as  a mirror  image  of  the  previous  cell  configuration 
at  the  same  general  location  of  the  present  runway.  Figure  6 shows  the 
relative  position  of  the  77  grid  cells  for  Conf iguration  I.  Table  8 is  a 
summary  of  the  vegetation  types  affected  by  Configuration  I. 

Configuration  11  demonstrates  the  effect  on  the  various  vegetation 
types  if  the  proposed  runway-launch  complex  was  shifted  approximately  6000 
feet  north  of  the  present  runway  location.  The  configuration  of  the  affected 
grid  cells  is  the  same  as  in  Figure  5.  Table  9 is  a summary  of  the  vegeta- 
tion types  affected  by  Configuration  II. 


In  addition  to  information  on  the  vegetation  types  affected,  as 
described  in  the  above  examples,  the  Search/Count  program  prints  out  sub- 
variables and  frequencies  for  the  soil  types,  exposure  categories,  and 
elevation  classes.  An  infinite  variety  of  configurations  for  the  runway- 
launch  complex  (or  other  developments)  can  be  tested  using  the  Search/Count 
program  depending  upon  the  needs  of  the  user. 


TABLE  8.  SUMMARY  OF  THE  VEGETATION  TYPES  AFFECTED  BY 
CONFIGURATION  I OF  THE  RUNWAY-LAUNCH  COMPLEX  BASED  ON  THE 
SEARCH/COUNT  PROGRAM. 


I 


L 

Vegetation  Type 

No.  of 
Subcel 1 s 

Estimated 

Acres 

[ 

Grassland-annual 

120 

306  1 

i 

Coastal  sage  scrub- 
normal  phase 

108 

l 

275 

K 

Chaparral 

285 

727 

Riparian  woodland- 
sparse  phase 

1 

3 

Riparian  woodland 

1 

3 

\ ;Ruderal  vegetation 

151 

385 

Man-made  facilities 

27 

69 

TOTALS 

693 

1768 

■ * 

1 

TABLE  9.  SUMMARY  OF  THE 
CONFIGURATION  II  OF  THE 
SEARCH/COUNT  PROGRAM. 

VEGETATION  TYPES  AFFECTED  BY 
RUNWAY-LAUNCH  COMPLEX  BASED  ON 

THE 

> 

Veqetation  Type 

No.  of 
Subcel 1 s 

Estimated 

Acres 

Coastal  sage  scrub- 
stabilized  dune  phase 

27 

69 

Grassland-annual 

276 

704 

Coastal  sage  scrub- 
normal  phase 

113 

288 

Ruderal  vegetation 

14 

36 

jChaparral 

250 

637 

Planted  trees 

4 

10 

Riparian  woodland 

2 

5 

i 

Riparian  woodland- 
sparse  phase 

1 

3 

1 

Man-made  facilities 

6 

15 

TOTALS 

693 

1 767  ' 
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SECTION  V 


COMPARISON  OF  MANUAI  AND  AUTOMATED  METHODS  FOR 
DETERMINING  AREAS  OF  VEGETATION 

The  purpose  of  this  study  is  to  test  several  methods  for  manually 
determining  areas  of  vegetation  and  statistically  compare  them  to  automated 
methods  using  the  digitized  data  base.  Three  different  scales  of  analysis 
were  used  to  compare  selected  areas  of  vegetation  --  two  manual  methods  and 
one  automated  method. 

1.  MANUAL  METHODS 

Two  manual  methods  were  used  for  determining  areas  of  vegetation  in  the 
following  study: 

a.  cutting  and  weighing  method 

b.  planimeter  method 

The  overlay  vegetation  maps  made  for  the  C series  of  maps  in  the  Base  Master 
Plan  (Reference  4)  were  used  in  the  study. 

2.  CUTTING  AND  WEIGHING  METHOD 

The  cutting  and  weighing  method  involves  cutting  out  (from  a copy  of 
the  vegetation  map)  each  distinct  area  of  vegetation  and  weighing  the  combined 
pieces  of  each  vegetation  type.  The  weight  of  each  vegetation  type  is  then 
compared  with  the  known  weight  and  area  of  the  vegetation  map  to  calculate 
the  area  of  each  vegetation  type.  A sensitive  balance  (+0.0001  gram)  must 
be  used  to  accurately  weigh  the  pieces  of  the  map. 

The  following  procedure  was  used  with  Vegetation  Map  Sheet  No.  61  using 
the  cutting  and  weighing  method: 

a.  Determine  the  total  area  of  the  map  sheet. 

b.  Weigh  the  whole  map  sheet. 

c.  Cut  out  the  individual  areas  of  vegetation  with  a sharp  razor 
blade  or  knife. 

d.  Weigh  the  combined  pieces  of  each  vegetation  type. 

e.  Calculate  the  areas  of  each  vegetation  type  by  equating  the  area 
and  weight  of  the  whole  map  with  the  individual  weights  of  each 
vegetation  type. 

Table  10  is  a summary  of  the  cutting  and  weighing  method  analysis  for 
Vegetation  Map  Sheet  No.  61.  The  areas  calculated  by  the  cutting  and 
weighing  method  compare  rather  closely  to  those  planimetered  for  the  same 
map  sheet  (see  Table  11). 

The  cutting  and  weighing  method  has  several  disadvantages  when  compared 
to  the  planimeter  method: 

a.  It  takes  longer  (depending  on  the  number  of  areas  to  be  cut  and 
their  irregularity). 

b.  There  are  more  calculations  and  measurements  required  to  estimate 
the  unknown  areas. 
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TABLE  10.  AREAS  OF  VEGETATION  BY  THE  CUTTING  AND  WEIGHING  METHOD 
FOR  VEGETATION  MAP  SHEET  NO.  61. 


Vegetation 

Wei ght 

Estimated  Area 

Estimated  Area 

Code  No.* 

(Grams ) 

(Square  Inches) 

(Acres ) 

25 

0.0052 

0.12 

2 

9 

0.0046 

0.10 

2 

12 

0.0180 

0.41 

7 

118 

0.0338 

0.77 

14 

6 

0.0746 

1.71 

31 

00 

1 . 7802 

40.73 

736 

1 

19 

0.0017 

0.04 

1 

5 

0.0153 

0.35 

6 

7 

1.0607 

24.27 

439 

16 

3.1711 

72.56 

1312 

UND** 

0.0169 

0.39 

7 

TOTALS 

6.1821 

141 .45 

2557 

Weight  of  whole  map  = 6.1825  grams 

i 

i 

Area  of  whole  map 

= 141.47 

square  inches 

1 

Total  time  for  analysis  = 

3.5  hours 

*See  Section  VII 

**Undetermi ned  vegetation 

types 

1 

r 

i * 

i 
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TABLE  11.  PLANIMETER  SUMMARY  FOR  VEGETATION  MAP  SHEET  NO.  61. 


Vegetation 
Code  No.* 

Area 

(Square  Inches) 

Estimated  Area 
(Acres) 

UND** 

0.44 

8 

12 

0.40 

7 

16 

73.05 

1321 

18 

0.65 

12 

5 

0.33 

6 

6 

2.05 

37 

19 

0.02 

0 

9 

0.13 

2 

7 

24.18 

437 

00 

41.21 

745 

25 

0.12 

2 

TOTALS 

142.70 

2577 

Total  time  for 

analysis  = 

3.2  hours 

*See  Section  VII 

**Undetermined 

vegetation 

types 

c.  A copy  of  the  map  must  be  sacrificed  to  the  procedure. 

d.  It  is  extremely  tedious  work. 

Because  of  the  reasons  given  above,  the  cutting  and  weighing  method  is 
unsuitable  for  determining  areas  of  vegetation  in  large  scale  studies. 

3.  PLANIMETER  METHOD 

The  planimeter  method  was  the  second  manual  method  used  for  determining 
areas  of  vegetation  in  this  study.  The  planimeter  used  was  calibrated  to 
give  area  measurements  in  square  inches.  Of  the  two  manual  methods  tested, 
the  planimeter  method  was  easier  and  quicker.  Four  selected  map  sheets  were 
planimetered  for  use  in  later  analyses.  Tables  11,  12,  13,  and  14  summarize 
the  planimeter  data  for  the  four  map  sheets. 

An  advantage  of  using  the  planimeter  method  over  automated  methods  is 
that  the  display  is  in  a form  familiar  to  those  accustomed  to  reading  maps. 
Disadvantages  of  using  the  planimeter  method  over  automated  methods  are: 

a.  Cost  and  time  of  preparing  and  pi animeteri ng  overlay  combinations 
limit  the  number  of  closely  related  alternatives  that  can  be 
evaluated  and  displayed. 

b.  Too  many  overlays  put  on  top  of  one  another  quickly  get  confusing 
unless  irrelevant  elements  are  screened  out,  which  is  difficult 
with  conventional  cartography. 

4.  AUTOMATED  METHOD 

The  GRID  computer  graphics  program  was  the  automated  method  used  for 
determining  and  displaying  areas  of  vegetation  in  this  study.  The  Search/ 
Count  program  and  the  Small  Cell  Count  program  were  used  with  the  data  base 
to  independently  determine  areas  for  computation.  Listings  of  the  Search/ 
Count  program  and  the  Small  Cell  Count  program  are  given  in  Appendixes  G and 
H. 

Advantages  of  the  automated  method  over  manual  methods  are: 

a.  Combinations  of  overlays  for  displaying  alternatives  can  be 
produced  almost  without  limit. 

b.  The  cost  is  minimal  once  the  data  base  and  the  computer  graphics 
program  are  established. 

Disadvantages  of  the  automated  method  include: 

a.  The  display  may  look  crude  and  unattractive  to  the  unaccustomed 
user. 

b.  The  resolution  is  not  as  fine  as  with  the  manual  method. 

5.  INVESTIGATION  OF  THREE  DIFFERENT  SCALES  OF  ANALYSIS 

This  study  compares  areas  of  vegetation  on  four  selected  map  sheets 
using : 

a.  Planimeter 

b.  Small  Cell  Count  program  (2.55-acre  grid  cells) 

c.  GRID  (aggregated  to  23  acres) 
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TABLE  12.  PLANIMETER  SUMMARY  FOR  VEGETATION  MAP  SHEET  NO.  28. 


'Vegetation 
Code  No.  * 

Area 

(Square  Inches) 

Estimated  Area 
(Acres) 

18 

2.96 

54 

UND** 

2.25 

41 

6 

17.79 

322 

9 

0.34 

6 

16 

63.57 

1149 

5 

49.31 

891 

19 

1.10 

20 

4 

3.13 

57 

3 

0.32 

6 

23 

0.83 

15 

22 

0.01 

0 

TOTALS 

141.61 

2561 

Total  time  for  analysis  = 6.0  hours 
*See  Section  VII 
**Undetermined  vegetation  types 
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TABLE  13.  PLANIMETER  SUMMARY  FOR  VEGETATION  MAP  SHEET  NO.  41. 


Vegetation 
Code  No.* 

Area 

(Square  Inches) 

Estimated  Area 
(Acres ) 

12 

0.26 

5 

13 

9.14 

165 

6 

11.94 

216 

16 

57.22 

1034 

14 

0.32 

6 

17 

3.80 

69 

5 

14.91 

270 

25 

0.11 

2 

22 

4.14 

75 

52 

2.06 

37 

4 

15.90 

287 

18 

4.79 

87 

23 

1.93 

35 

9 

0.05 

1 

UND** 

1.35 

24 

20 

13.42 

243 

19 

0.34 

6 

42 

0.41 

7 

TOTALS 

142.09 

2569 

Total  time  for  analysis  = 4.5  hours 
*See  Section  VII 
**Undetermined  vegetation  types 


TABLE  14.  PLANIMETER  SUMMARY  FOR  VEGETATION  MAP  SHEET  NO.  62. 


'Vegetation 
Code  No.* 

(Sq 

Area 

uare  Inches) 

Estimated  Area 
(Acres) 

7 

43.40 

785 

16 

91 .40 

1652 

6 

4.74 

86 

9 

1 .02 

18 

5 

0.78 

14 

52 

0.95 

17 

00 

0.10 

2 

19 

0.05 

1 

UND** 

0.09 

2 

22 

0.13 

2 

TOTALS 

142.66 

2579 

Total  time  for 

analysis  = 

4.0  hours 

*See  Section  VII 

**Undetermined 

vegetation 

types 
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The  map  sheets  used  in  the  study  are  described  in  Section  V #3.  The 
Small  Cell  Count  program  was  simulated  for  the  map  sheets  by  tallying  the 
vegetation  types  in  the  2.55-acre  grid  cells  from  a printout  of  the  data 
base.  GRID  (aggregated  to  23-acre  cells)  was  simulated  using  the  same 
procedure.  The  vegetation  type  of  the  center  subcell  determined  the 
vegetation  type  for  the  23-acre  grid  cell. 

The  grid  cell  counts  were  converted  to  acres  for  comparison  of  areas  of 
vegetation  for  the  three  scales  of  analysis.  Table  15  is  the  pooled  data  of 
the  four  vegetation  map  sheets.  The  vegetation  types  were  ranked  from 
largest  to  smallest  in  numbers  of  acres,  using  the  estimated  acres  planimetered 
for  each  vegetation  type.  The  percent  deviation  (small  cells)  was  computed 
by  the  following  formula: 

percent  deviation  (small  cells)  = 

acres  (small  cells)  - acres  (planimetered) 

acres  (planimetered)  ^ 

Acres  computed  for  the  planimetered  data  were  assumed  to  be  more  accurate  than 
acres  computed  by  the  grid  method. 

The  percent  deviation  (large  cells)  was  computed  by  a similar  formula: 

percent  deviation  (large  cells)  = 

acres  (large  cells)  - acres  (planimetered)  , nn  nprrpn, 
acres  (planimetered)  ^ 

Comparing  the  percent  deviations  for  both  the  small  and  large  cells 
gives  a common  trend  of  less  deviations  between  the  three  scales  of  analysis 
for  estimating  large  areas.  As  the  area  decreases,  the  percent  deviations 
increase  for  both  the  small  cell  and  large  cell  grid  data.  For  this  parti- 
cular study,  there  was  less  than  7 percent  deviation  between  the  three 
different  scales  of  analysis  for  areas  greater  than  1200  acres. 

The  weighted  mean  deviation  was  used  to  compare  the  percent  deviation 
(small  cells)  to  the  percent  deviation  (large  cells)  when  their  values  were 
weighted  by  using  the  planimeter  estimated  acres  for  each  vegetation  type. 

The  weighted  mean  deviation  was  calculated  by  the  following  formula: 

...j.ui.j  _ l(acres  (planimetered)  x ^percent  deviation!) 

weighted  mean  deviation  = — i ^ — — rS — -t — -n 

3 lacres  (planimetered) 

Using  the  data  in  Table  15,  the  following  values  were  calculated: 

weighted  mean  deviation  (small  cells)  = 3.3 

weighted  mean  deviation  (large  cells)  = 7.1 

The  above  values  indicate  that  overall,  there  is  more  variability  when 
estimating  acres  of  vegetation  using  the  aggregated  23-acre  GRID  cell,  as 
compared  to  the  small  cell  count  method. 
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TABLE  15.  POOLED  DATA  OF  FOUR  VEGETATION  MAP  SHEETS  USING  THREE 
DIFFERENT  SCALES  OF  ANALYSIS. 


■Vegetation 

! Type  * 
l(Code  No. ) 

Planimeter 

Estimated 

Acres 

Small  Cells 
Estimated 
Acres 

Percent 
Devi ati on 
Sma  1 1 
Cells 

Large  Cells 
Estimated 
Acres 

Percent 

Deviation 

Large 

Cells 

16 

5157 

5182 

0.5 

5250 

2.4 

7 

1222 

1242 

1 .6 

1240 

1.5 

5 

1181 

1247 

5.6 

1263 

6.9 

00 

747 

719 

-3.7 

620 

-17.0 

6 

660 

701 

6.2 

735 

11.4 

A 

344 

367 

6.7 

344 

0.0 

20 

243 

237 

-2.5 

230 

-5.3 

13 

165 

173 

4.8 

230 

39.4 

18 

152 

130 

-14.5 

161 

5.9 

22 

77 

64 

-16.9 

23 

-70.1 

17 

69 

61 

-11.6 

23 

-66.7 

52 

54 

36 

-33.3 

46 

-14.8 

23 

50 

59 

18.0 

46 

-8.0 

9 

28 

0 

-100.0 

0 

-100.0 

19 

27 

23 

-14.8 

0 

-100.0 

12 

12 

13 

8.3 

23 

91.7 

42 

7 

3 

-57.1 

0 

-100.0 

3 

6 

0 

-100.0 

0 

-100.0 

14 

6 

8 

33.3 

0 

-100.0 

25 

4 

10 

150.0 

23 

475.  u 

8 

0 

8 

-- 

0 

-- 

UND** 

75 

-- 

-- 

-- 

-- 

TOTALS 

— 

10286 

10283 

10287 

*See  Section  VII 

**Undetermi ned  vegetation  types  (from  vegetation  map  sheets) 
Weighted  mean  deviation  (small  cells)  = 3.3 
Weighted  mean  deviation  (large  cells)  = 7.1 
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The  nonparametric  Spearman  rank  correlation  method  (Reference  7)  was 
used  to  compare  the  pooled  data  in  Table  15.  The  computing  formula  for 
the  Spearman  rank  correlation  coefficient  is: 

1 - 6zd^ 


The  value  of  r may  range  from  -1  to  +1.  The  interpretation  of  r is  that 
as  its  value  approaches  +1,  the  two  groups  have  increasingly  similar  rank 
correlation. 


Two  Spearman  rank  correlation  tests  were  performed  on  the  data  in  Table 
15.  In  the  first  test,  acres  (planimetered)  were  compared  to  acres  (small 
cells)  for  the  different  vegetation  types.  The  calculated  value  of  rs  is: 

r$  = 0.958 

In  the  second  test,  acres  (planimetered)  were  compared  to  acres  (large 
cells)  for  the  different  vegetation  types.  The  calculated  value  is: 

rs  = 0.920 

Both  tests  showed  highly  significant  rank  correlations.  The  critical 
value  of  the  Spearman  rank  correlation  coefficient  for  the  above  two  tests 
is: 


rs. 05(2)20  0-447  P < 0-001 

The  interpretation  which  can  be  made  from  the  Spearman  rank  correlation 
tests  above  is  that  both  methods  of  estimating  areas  of  vegetation  (small 
cell  and  large  cell ),  when  compared  to  planimetered  areas,  show  highly 
significant  correlation  when  their  areas  are  ranked.  The  percent  deviation 
of  the  planimetered  areas  from  the  small  cell  and  large  cell  areas  must  be 
considered  to  view  the  Spearman  rank  correlation  test  in  perspective. 

6.  COMPARISON  OF  THE  SMALL  CELL  COUNT  PROGRAM  AND  GRID 

This  study  compares  areas  of  vegetation  for  the  entire  base  as  estimated 
by  the  small  cell  (subcell)  count  program  and  GRID  (aggregated  to  23-acre 
cells).  Two  different  aggregation  methods  were  used.  In  the  first  method, 
the  computer  was  directed  to  determine  the  predominant  vegetation  type  of  a 
grid  cell  from  the  vegetation  types  of  the  nine  subcells.  In  the  second 
method,  the  computer  was  directed  to  determine  the  vegetation  type  for  a grid 
cell  by  using  the  vegetation  type  of  the  center  subcell.  Table  16  summarizes 
the  data.  The  grid  cell  counts  were  converted  to  acres  for  comparison. 

The  vegetation  types  in  Table  16  were  ranked  from  largest  to  smallest 
(in  numbers  of  acres)  using  estimated  acres  by  small  cells.  The  differences 
between  the  totals  for  both  GRID  methods  was  due  to  rounding  error. 


Percent  deviation  was  computed  by  the  formula: 
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TABLE  16.  AREAS  OF  VEGETATION  FOR  VANDENBERG  AFB. 


Vegetation  Small  Cells 
Type  Estimated 

(Code  No. )*  Acres 

GRID 

(Majori ty) 
Estimated 
Acres 

Percent 

Deviation 

GRID 

(Majori ty) 

GRID 

(Center) 

Estimated 

Acres 

Percent 

Deviation 

GRID 

(Center) 

16 

38610 

39422 

2.1 

38389 

-0.6 

6 

19337 

19746 

2.1 

191C3 

-1.2 

5 

14359 

14671 

2.2 

14350 

-0.1 

8 

9091 

9092 

0.0 

8931 

-1.8 

3 

4149 

4294 

3.5 

4248 

2.4 

23 

3815 

3697 

-3.1 

3903 

2.3 

7 

3302 

3421 

3.6 

3398 

2.9 

18 

2512 

2388 

-4.9 

2663 

6.0 

00 

2259 

1010 

-55.3 

2319 

2.7 

20 

1907 

1975 

3.6 

1929 

1.2 

12 

1658 

1676 

1.1 

1883 

13.6 

4 

1617 

1561 

-3.5 

1676 

3.6 

42 

436 

390 

-10.6 

528 

21.1 

24 

431 

436 

1.2 

321 

-25.5 

72 

357 

253 

-29.1 

298 

-16.5 

14 

352 

321 

-8.8 

321 

-8.8 

99 

321 

321 

0.0 

321 

0.0 

52 

283 

184 

-35.0 

321 

13.4 

13 

280 

321 

14.6 

321 

14.6 

19 

273 

298 

9.2 

184 

-32.6 

25 

255 

253 

-0.8 

230 

-9.8 

1 

227 

207 

-8.8 

321 

41 .4 

22 

194 

184 

-5.2 

161 

-17.0 

17 

166 

161 

-3.0 

115 

-30.7 

11 

158 

115 

-27.2 

207 

31.0 

9 

127 

161 

26.8 

92 

-27.6 

2 

61 

46 

-24.6 

69 

13.1 

26 

43 

46 

7.0 

46 

7.0 

62 

33 

0 

-100.0 

23 

-30.3 

10 

8 

0 

-100.0 

0 

-100.0 

TOTALS 

106621 

106650 

106671 

*See  Section  VII 

1 

Weighted  mean  deviation  GRID  (majority)  = 3.7 
Weighted  mean  deviation  GRID  (center)  = 2.1 
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percent  deviation 


acres(GRID)  - acres(small  cells) 
acres(sma11  cells) 


x 100  percent 


Acres  computed  by  the  small  cell  count  method  were  assumed  to  be  more  accurate 
than  acres  computed  by  the  aggregated  GRID  method.  Comparing  the  three 
methods,  there  is  less  than  4 percent  deviation  for  areas  decreasing  to  3500 
acres.  In  general,  the  percent  deviations  increase  as  the  areas  of  the 
vegetation  types  decrease. 


The  weighted  mean  deviation  was  used  to  compare  the  percent  deviation 
GRID  (majority)  to  the  percent  deviation  GRID  (center)  when  their  values  were 
weighted  by  using  the  small  cells  estimated  acres  for  each  vegetation  type 
in  Table  16.  See  Section  V #5  for  the  computing  formula  for  the  weighted 
mean  deviation. 


Using  the  data  in  Table  16,  the  following  values  were  calculated: 

weighted  mean  deviation  GRID  (majority)  = 3.7 
weighted  mean  deviation  GRID  (center)  = 2.1 

The  values  above  indicate  that  there  is  more  overall  variability  for  estimating 
acres  of  vegetation  using  estimated  acres  GRID  (majority)  compared  to 
estimated  acres  GRID  (center). 


The  nonparametric  Spearman  rank  correlation  method  was  used  to  compare 
areas  of  vegetation  for  the  enti re  base  as  estimated  by  small  cells  and  the 
aggregated  GRID  method.  Two  Spearman  rank  correlation  tests  were  performed 
on  the  data  in  Table  16.  In  the  first  test,  acres  (small  cells)  were 
compared  to  acres  (GRID,  aggregated  by  the  majority  of  subcells)  for  the 
different  vegetation  types.  The  calculated  value  of  r$  is: 

rs  = 0.983 

In  the  second  test,  acres  (small  cells)  were  compared  to  acres  (GRID, 
aggregated  by  the  center  subcell  vegetation  type).  The  calculated  value  is: 


r$  = 0.980 

The  calculated  values  of  r in  both  tests  showed  highly  significant 
rank  correlations.  The  critical  value  of  rg  for  the  above  two  tests  is: 


rs. 05(2)30 


0.362 


P < 0.001 


The  interpretation  from  the  Spearman  rank  correlation  tests  is  that  both 
GRID  methods  of  estimating  areas  of  vegetation  for  the  entire  base  when 
compared  to  areas  computed  by  the  small  cell  method  show  highly  significant 
correlation  when  their  areas  are  ranked.  The  calculated  Spearman  rank 
correlation  coefficients  are  almost  identical  for  the  two  comparisons. 

The  percent  deviations  for  the  two  aggregation  methods  show  close 
similarity,  especially  for  the  larger  areas.  As  the  areas  decrease, 
departures  from  the  areas  estimated  by  small  cells  increase.  The  aggregated 
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GRID  method  by  the  majority  of  subcells  seems  to  be  the  most  severely  affected. 
This  is  probably  due  to  the  patchiness  of  the  scarce  vegetation  types. 

The  best  aggregation  method  to  display  vegetation  for  GRID  studies, 
based  on  the  data  in  Table  16,  is  the  center  subcell  method.  This  method 
can  be  modified  to  more  accurately  display  areas  of  scarce  vegetation  (less 
than  200  acres)  by  having  GRID  display  areas  of  partial  occurrence  of  the 
scarce  vegetation  types.  More  research  is  needed  on  this  aspect  of 
aggregation. 
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SECTION  VI 
USERS  MANUAL 


1.  INTRODUCTION 

The  following  Users  Manual  is  written  specifically  for  the  use  of  the 
GRID  computer  graphics  program  and  the  Search/Count  program  with  the 
environmental  data  base  developed  for  VAFB.  The  GRID  Manual  (Reference  8) 
was  used  as  a key  reference  in  writing  sections  of  this  Users  Manual 
pertaining  to  GRID. 

The  GRID  program  was  developed  at  the  Harvard  Laboratory  for  Computer 
Graphics  and  Spatial  Analysis.  It  provides  an  efficient  means  for  the 
graphic  display  of  large  quantities  of  information  collected  on  the  basis 
of  a rectangular  coordinate  grid.  It  has  the  capacity  of  manipulating  up 
to  10,000  different  grid  cells. 

The  Search/Count  program  is  essentially  a frequency  counting  program. 

It  can  be  used  in  conjunction  with  GRID  or  independently  to  locate  and  print 
out  data  base  variables  for  selected  grid  cells. 

GRID  and  the  Search/Count  program  are  written  in  FORTRAN  IV.  They  are 
operational  on  the  IBM  360  computer  at  San  Diego  State  University  and  the 
Burroughs  3500  computer  at  VAFB.  The  computer  programs  in  the  EPS  are 
designed  for  computer  users  rather  than  computer  programmers.  They  require 
the  user  to  have  a basic  knowledge  of  FORTRAN  IV. 

Complete  listings  for  the  GRID  program  and  the  Search/Count  program  are 
given  in  Appendixes  A and  G.  Examples  in  the  Users  Manual  and  programs 
listings  in  the  appendixes  are  specifically  adapted  for  the  Burroughs  3500 
computer.  The  user  is  advised  to  consult  with  a programmer  to  adapt  the 
computer  programs  to  other  computer  systems. 

2.  BASIC  PRINCIPLES  OF  GRID 

The  GRID  program  associates  a data  value  with  a cell  on  a grid.  The 
program  processes  the  data  in  the  order  in  which  it  prints  the  map. 
Therefore,  the  data  must  be  in  correct  order.  The  program  begins  processing 
the  data  starting  with  the  grid  cell  in  the  northwest  corner  of  the  map  and 
continues  horizontally  row  by  row  and  from  left  to  right  in  each  row.  The 
numbers  below  represent  the  order  in  which  the  first  40  grid  cells  are 
processed  and  printed  for  GRID  maps  of  VAFB: 

1 2 

3 4 5 6 

7 89  10  11  12  13  14  15 

16  17  18  19  20  21  22  23  24  25  26  27 

28  29  30  31  32  33  34  35  36  37  38  39  40 

To  obtain  a map,  the  user  must  provide  the  following  packages  or  sets  of 

instructions:  Subroutine  Flexin,  Irregular  Outline  Package,  Map  Package, 

and  Data  Input  Package. 
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3.  SUBROUTINE  FLEXIN 

Flexin  is  the  FORTRAN  subroutine  in  GRID  which  is  modified  by  the 
user  to  produce  a map.  Flexin  is  called  by  the  main  program  once  for  each 
data  cell  that  is  mapped.  The  user  specifies  two  sets  of  instructions  in 
Flexi n : 

a.  Where  the  data  value  is  located  on  the  data  card  or  tape. 

b.  What  analysis,  if  any,  is  to  be  performed  on  the  data  variable 
or  variables  to  derive  the  value  that  is  mapped. 

The  subroutine  statement  for  Flexin  is: 

SUBROUTINE  FLEX  I N ( I FORM,  T,  FIRST  ) 

The  variables  IFORM  and  FIRST  are  carried  into  the  subroutine  as  control 
variables.  The  data  value  to  be  mapped  is  carried  back  to  the  main  program 
as  T (once  for  every  value  to  be  mapped).  T is  mapped  according  to  the 
instructions  specified  in  the  Map  Package.  The  following  examples  (1,  2 and 
3)  demonstrate  the  use  of  IFORM  and  T.  The  variable  FIRST  can  be  used  as  a 
logical  variable  in  Subroutine  Flexin.  It  is  true  on  the  first  entry  to 
Flexin  and  false  on  all  other  entries  from  the  main  program.  This  use  of 
FIRST  is  automatically  set  in  the  GRID  program  outside  of  Subroutine  Flexin. 

Example  1 shows  how  to  use  Flexin  to  map  the  elevation  data  coded  for 
VAFB.  This  is  a relatively  simple  use  of  Flexin. 

Example  1 : 

A Simple  Use  of  Subroutine  Flexin 

SUBROUTINE  FLEXIN  ( IFORM,  T,  FIRST  ) 

READ ( 11,1)  EL 

1 FORMAT ( 16X , F1.0) 

T = EL 

RETURN 

END 

This  subroutine  instructs  the  computer  to  read  the  variable  EL  from  the 
data  tape.  The  FORMAT  statement  specifies  where  the  variable  is  on  the 
data  tape.  The  fourth  statement  assigns  the  variable  T the  elevation  value 
assigned  to  variable  EL.  The  elevation  value  to  be  mapped  for  the  grid 
cell  is  returned  to  the  main  program  as  the  value  of  T in  the  subroutine 
statement.  Flexin  is  called  once  for  each  grid  cell  that  is  mapped.  A new 
value  for  T is  read  and  returned  to  the  main  program  for  the  value  to  be 
mapped  for  each  grid  cell.  T is  mapped  according  to  instructions  given  in 
the  Map  Package.  For  example,  the  Map  Package  could  be  set  up  to  map  the 
elevation  data  into  10  levels  --  for  the  10  possible  elevation  codes. 

Example  2 shows  how  to  use  Flexin  in  a more  detailed  analysis  to  produce 
a map  which  will  display  areas  of  suitable  habitat  for  the  California  Legless 
Lizard.  Both  the  soil  and  vegetation  data  were  used  in  making  the  map. 

Figure  4 shows  the  GRID  map  for  tin's  example. 
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Example  2: 

A Complex  Use  of  Subroutine  Flexin 


Explanatory  Notes 


SUBROUTINE  FLEX  I N ( IFORM,  T,  FIRST) 

REAL  SOIL (22)  a. 

DATA  SOIL/4. ,5. ,6. ,25. ,28. ,59. ,60. ,63. ,64. ,76. ,113. , 

11  14. ,133. ,134. ,138. ,168. ,170. ,206. ,207. ,208.  ,217.  ,71 ./  D- 

NTYPE=  22 

RE AD (11 ,11)  S,V  c. 

11  FORMAT ( 3X , F3 . 0 , 30X , F2 . 0 ) 

T=1 0 . d. 

DO  300  J=  1 ,NTYPE 
I F ( S . EQ. SOIL ( J ) ) GO  TO  450 
300  CONTINUE 

I F ( V . EQ .12.)  T=  1. 

IF( V . EQ. 8. ) T=2. 

IF ( V . EQ .11.)  T=3. 

IF( V . EQ. 3)  T=4. 

IF ( V . EQ. 31 . ) T=4. 

GO  TO  500 

450  IF ( V . NE . 12. ) T=5. 

IF (V . NE . 8. ) T=5. 

I F ( V . NE .11.)  T=5. 

IF(V.NE. 3. ) T=5. 

IF( V . NE . 31 . ) T=5. 

IF (V . EQ . 1 2 . ) T=6. 

IF (V . EQ. 8. ) T=7. 

IF (V . EQ. 1 1 . ) T=8. 

IF (V . EQ. 3. ) T=9. 

IF(V . EQ. 31 . ) T=9 . e. 

500  CONTINUE 
RETURN 
END 

Explanatory  Notes  for  Example  2: 

a.  This  statement  creates  storage  space  for  22  soil  codes. 

b.  Twenty-two  soil  types  were  identified  as  suitable  for  the  Legless  Lizard. 
These  soil  codes  are  entered  into  the  storage  spaces  created  in  a. 

c.  A soil  code  and  vegetation  code  are  read  for  a grid  cell  from  the  data 
tape. 

d.  The  remainder  of  the  subroutine  compares  criteria  for  a suitable  habitat 
for  the  Legless  Lizard  with  each  vegetation  and  soil  code  in  the  data 
base.  Ten  different  levels  were  selected  for  printing  the  map.  Each  grid 
cell  is  classified  into  one  of  the  levels  based  on  its  vegetation  and 
soil  codes.  The  level  is  returned  to  the  main  program  as  the  value  of 

T.  The  levels  for  suitable  habitat  are: 

Level  1 = Coastal  strand  vegetation 

Level  2 = Coastal  sage  scrub-stabilized  dune  phase  vegetation 
Level  3 = Coastal  bluff  vegetation 
Level  4 = Oak  woodland  vegetation 
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Level  5 = Suitable  soil  associations 

Level  6 = Coastal  strand/suitable  soil  overlap 

Level  7 = Coastal  sage  scrub-stabilized  dune  phase/suitable  soil  overlap 
Level  8 = Coastal  bluff/suitable  soil  overlap 
Level  9 = Oak  woodland/suitable  soil  overlap 

Level  10  = Areas  which  are  not  suitable  habitat  for  the  Legless  Lizard 

e.  Once  the  value  of  T is  determined  for  a grid  cell  in  Subroutine  Flexin, 
it  is  mapped  according  to  instructions  provided  in  the  Map  Package. 

Example  3 shows  a relatively  simple  use  of  the  control  variable  IFORM, 
if  more  than  one  map  is  desired  with  one  computer  submission.  In  this  example, 
the  computer  is  directed  to  produce  two  maps. 

Example  3: 

Use  of  IFORM  With  Subroutine  Flexin 


GO  TO  ( 1,  2 ),  IFORM 

1 CONTINUE 
READ (11  ,15)  EL 

15  FORMAT ( 16X.F1.0) 

T = 2. 

IF  EL.FQ.  0.  T=1 . 

IF  EL.EQ.  1.  T=1 . 

IF  EL.EQ.  2.  T=1 . 

RETURN 

2 CONTINUE 
READ ( 1 1,15)  EL 
T=2 . 

IF  EL.EQ. 6.  T=1 . 

IF  EL.EQ. 7.  T=1 . 

IF  EL.EQ. 8.  T=1 . 

IF  EL.EQ. 9.  T=1 . 

RETURN 
END 

Explanatory  Notes  for  Example  3: 

a.  The  variable  IFORM  is  carried  into  the  subroutine  as  a control  variable. 
The  value  to  be  mapped  is  carried  back  to  the  main  program  as  T (Flexin 
is  called  once  for  every  value  to  be  mapped). 

b.  Since  2 maps  are  to  be  made,  Flexin  is  split  into  2 segments.  The 
variable  IFORM  is  given  its  value  in  Elective  2,  field  1 of  the  Map 
Package  for  each  map.  Two  map  packages  are  required  to  produce  2 maps. 
After  reading  the  GO  TO  statement,  the  program  jumps  to  the  statement 
N CONTINUE,  where  N is  the  value  of  IFORM.  If  IFORM  equals  2,  as 
specified  in  the  Map  Package,  the  program  will  execute  statement  2 
CONTINUE  and  the  statements  following  it. 

In  Example  3 there  are  only  2 routines  in  Subroutine  Flexin.  By  extending 
the  GO  TO  statement,  as  many  routines  as  needed  may  be  used.  Each 
routine  is  sandwiched  between  a CONTINUE  statement  (which  indicates  the 
beginning  of  each  routine)  and  a RETURN  statement  (which  sends  the  value 
T back  to  the  main  program). 


Explanatory  Notes 

a. 

b. 

c. 


d. 
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c. 


d. 


This  routine 
in  2 levels, 
feet.  Level 
600  feet. 

This  routine 
in  2 levels. 
Level  2 displ 


instructs  the  computer  to  produce  a map  displaying  elevation 
Level  1 displays  all  grid  cells  coded  from  zero  to  600 
2 displays  all  grid  cells  coded  at  elevations  greater  than 

instructs  the  computer  to  produce  a map  displaying  elevation 
Level  1 displays  all  grid  cells  at  1200  feet  or  greater, 
ays  all  grid  cells  less  than  1200  feet. 


4.  IRREGULAR  OUTLINE  PACKAGE 

This  data  set  specifies  the  shape  of  the  outline  for  the  grid  map.  The 
irregular  outline  data  determine  the  number  of  cells  from  the  left  and  right 
sides  of  the  map,  in  each  row,  that  are  left  blank.  The  following  data  cards 
are  used  to  make  the  Irregular  Outline  Package: 


a.  On  the  first  card,  IRREGULAR  OUTLINE,  is  punched  in  columns  1-17. 

b.  On  the  last  card,  99999,  is  punched  in  columns  1-5. 

c.  Between  the  first  and  last  card,  a series  of  cards  with  the 

following  format  are  punched: 

(1)  In  columns  1-5,  the  number  of  successive  rows  for  which  the 
format  is  repeated. 

(2)  In  columns  6-10,  the  number  of  blank  cells  at  the  beginning 
of  the  row. 

(3)  In  columns  11-15,  the  number  of  blank  cells  at  the  end  of 
the  row. 


The  numbers  above  must  be  punched  as  integer  numbers.  They  must  be 
right  justified  and  contain  no  decimal  points.  The  program  processes  the 
cards  in  order.  The  first  card  refers  to  the  top  row  (or  rows,  as  specified 
in  columns  1-5).  The  second  card  refers  to  the  second  row  (or  rows)  etc. 

The  irregular  outline  is  the  first  package  to  prepare  in  doing  a GRID 
study.  The  irregular  outline  for  VAFB  was  prepared  from  the  Base  Master 
Plan  maps  (Reference  4)  and  U.S.  Geological  Survey  maps.  The  C-1.2  series 
of  the  Base  Master  Plan  maps  are  blocked  into  1000-foot  grid  cells.  This 
made  it  relatively  easy  to  determine  the  irregular  outline  to  enclose  the 
base  boundaries. 

Once  the  irregular  outline  was  determined,  it  was  coded  on  to  coding 
forms  and  the  data  cards  were  keypunched.  The  irregular  outline  data  deck 
contains  102  data  cards  to  specify  the  146  rows  of  grid  cells  used  for  the 
VAFB  GRID  map.  Example  4 shows  how  the  data  deck  is  arranged  for  the  VAFB 
irregular  outline.  See  Appendix  B for  the  complete  listing  of  the  VAFB 
irregular  outline. 

Example  4: 

Arrangement  of  the  Data  Deck  for  the  VAFB  Irregular  Outline 
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column  mum 

12345678901234567 

IRREGULAR  OUTLINE 
1 7 59 

1 7 57 

102  data  cards 

2 45  18 

99999 

Once  the  irregular  outline  was  determined  and  keypunched,  it  was  used 
for  all  successive  GRID  maps.  The  irregular  outline  is  not  usually  changed 
in  a study  area  once  it  is  set.  If  the  irregular  outline  is  changed,  the 
data  coded  for  the  grid  cells  must  be  readjusted  to  fit  in  the  new 
i rregular  outl i ne. 

5.  MAP  PACKAGE 

This  package  instructs  the  computer  on  how  to  make  the  GRID  map.  A set 
of  map  electives  are  used  to  specify  the  instructions  for  making  the  map. 

a.  Preparation  of  the  Package 

The  map  package  is  prepared  according  to  the  following  procedure: 

(1)  On  the  first  card  of  the  package,  punch  MAP  in  columns  1-3. 

(2)  On  the  last  card  of  the  package,  punch  99999  in  columns  1-5. 

(3)  On  the  second,  third,  and  fourth  cards,  punch  the  title  you 

wish  to  have  appear  below  the  map.  One  or  more  of  these  cards 
may  be  left  blank,  but  all  3 cards  must  be  included. 

(4)  On  the  cards  to  be  inserted  between  the  fourth  and  last  card, 

the  desired  electives  are  punched.  Whenever  a map  elective 
is  not  called  for,  the  standard  result  described  under  each 
elective  will  automatically  occur.  Electives  1 and  7 must  be 
included  because  there  is  no  standard  condition  created  by  the 
program  for  them. 

b.  Standard  Format  for  the  Electives 

A standard  format  is  used  for  the  electives  with  the  exception  of 
electives  7,  10,  and  13.  The  format  is: 

(1)  The  elective  number  is  punched  as  an  i nteger  in  columns  4 and 
5 (right  justified). 

(2)  Six  fields  of  10  columns  each  are  used  for  specifying  the 
electives: 

Field  Columns 

1 11-20 

2 21-30 

3 31-40 

4 41-50 

5 51-60 

6 61-70 
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The  numbers  punched  in  any  of  the  six  fields  are  real  numbers  and  should 
contain  a decimal  point.  The  numbers  may  be  punched  anywhere  within  the 
columns  assigned  to  that  field.  The  decimal  point  may  be  omitted  if  the 
number  is  right  justified. 

c.  Elective  1 : 

Grid  (1  card) 

The  elective  specifies  the  parameters  for  the  rectangular  grid  that 
is  to  be  mapped.  Infield  1 specify  the  number  of  rows  of  grid  cells  down 
the  map.  In  field  2 specify  the  number  of  columns  of  grid  cells  across  the 
map. 


Fields  3 and  4 specify  the  size  of  each  printed  grid  cell.  Field  3 
specifies  the  number  of  characters  down  and  field  4 specifies  the  number  of 
characters  across.  If  fields  3 and  4 are  left  blank,  the  printer  cell  size 
will  be  4 x 5 or  \ inch  square. 

For  the  VAFB  irregular  outline,  a grid  cell  of  1 character  (lxl) 
produces  a map  of  approximately  8!2  by  19^  inches  when  8 characters  per  inch 
are  printed  by  the  line  printer.  This  map  is  slightly  distorted  vertically 
since  the  line  printer  will  print  10  characters  per  inch  horizontally.  The 
1 x 1 map  is  suitable  for  most  studies. 

A larger  map  with  slightly  less  distortion  can  be  produced  when  the 
grid  cell  size  is  2 x 3 (2  characters  vertically  by  3 characters  horizontally). 
The  line  printer  is  set  at  6 characters  per  inch  vertically  for  this  map. 

The  map  is  produced  in  2 parts  by  the  GRID  program.  The  halves  must  be  taped 
together  to  form  a complete  map.  The  map  measures  approximately  22  by  49*2 
inches  when  measured  from  the  bordering  row  and  column  numbers.  The  2x3 
grid  cell  map  takes  approximately  twice  as  much  computer  time  to  produce  as 
the  1 x 1 grid  cell  map.  Example  5 shows  how  to  punch  Elective  1 for  a 2 x 
3 grid  cell  size  for  the  VAFB  irregular  outline. 

Example  5: 

Elective  1 

Column  5 11-20  21-30  31-40  41-50 

1 146.0  68.0  2.0  3.0 

d.  Elective  2: 

Data  (1  card) 

This  elective  controls  the  input  options  for  the  data.  Data  Option 
A (see  section  on  the  data  package)  was  used  for  the  VAFB  data  base.  To 
activate  Data  Option  A,  specify  a number  greater  than  zero  in  field  1.  The 
number  specified  in  field  1 is  transferred  to  Subroutine  Flexin  as  the  value 
of  IF0RM.  The  use  of  IF0RM  is  discussed  in  the  section  on  Subroutine  Flexin. 
Example  6 shows  how  Elective  2 was  used  for  the  Vandenberg  study. 

Example  6: 

Elective  2 

Col  umn 5 20 

2 2 
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e.  Elective  3: 

Number  of  Levels  (1  card) 

This  elective  is  used  to  specify  the  number  of  levels  or  class 
intervals  used  for  the  GRID  map  (from  2 to  10).  The  number  desired  is 
punched  as  a decimal  number  in  field  1.  Standard  is  10  levels.  Example 
7 shows  how  to  punch  Elective  3 for  an  8 level  map. 

Example  7: 

Elective  3 

Column  5 11-20 

3 8.0 

f.  Elective  4: 

Value  Range  Minimum  (1  card) 

This  elective  is  used  to  specify  a number  as  the  minimum  value  of 
the  total  value  range  of  the  data.  The  number  is  punched  as  a decimal 
number  in  field  1.  Standard  is  to  use  the  minimum  value  of  the  data. 

g.  Elective  5: 

Value  Range  Maximum  (l  card) 

This  elective  is  used  to  specify  a number  as  the  maximum  value  of 
the  total  value  range  of  the  data.  The  number  is  punched  as  a decimal  number 

in  field  1.  Standard  is  to  use  the  maximum  value  of  the  data. 

h.  Elective  6: 

Value  Range  Intervals  (1  to  2 cards) 

This  elective  controls  the  value  range  for  each  level  or  interval. 

The  total  value  range  of  the  data,  as  modified  by  the  minimum  and  maximum 

of  the  data,  (Electives  4 and  5),  will  be  divided  up  into  the  number  of 
levels  specified  in  Elective  3.  Standard  is  to  have  each  level  or  interval 
assigned  an  equal  range. 

To  specify  the  desired  range  for  each  level,  use  values  proportionate 
to  the  size  of  the  desired  ranges.  These  should  be  punched  as  decimal 
numbers.  See  Example  8.  If  there  are  more  than  6 levels,  continue  on  a 
second  card,  punching  the  number  for  the  seventh  level  in  field  1,  for  the 
eighth  level  in  field  2,  etc. 

Example  8: 

Unequal  Value  Range  Levels 

The  data  are  to  be  divided  in  4 levels:  the  lowest  10  percent,  the 

next  25  percent,  the  next  35  percent,  and  the  remainder.  One  card  for 
Elective  6 would  be  punched  as  follows: 

Column  5 11-20  21-30  31-40  41-50 

6 10.  25.  35.  30. 

i . Electi ve  7 : 

Symbolism  (5  cards) 

This  elective  specifies  the  symbolism  that  will  be  printed  on  the 
map.  No  standard  symbolism  is  stored  in  the  program.  Therefore,  this 
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elective  must  be  included  in  the  map  package.  All  5 cards  must  be  included 
each  time  the  program  is  run. 

On  card  1 , punch  the  elective  number  7 in  column  5. 

On  card  2,  punch  as  follows  (any  print  characters  may  be  used  for  symbolism): 

Columns  1-10  are  used  to  specify  the  general  symbolism  for  each  level 
(column  1 for  the  symbol  to  designate  the  first  level,  etc.  --  for  as 
many  levels  as  are  to  be  used). 

Columns  11-20  are  used  to  specify  the  special  symbolism  for  the  flag 
points  (column  11  for  the  symbol  to  designate  flag  points  in  the 
first  level,  etc.).  The  flag  point  is  the  central  character  of  a 
gri d cell. 

Column  21  is  used  to  specify  the  symbolism  for  a value  less  than  the 
minimum  specified  in  Elective  4. 

Column  22  is  used  to  specify  the  flag  point  symbolism  for  a low  va'ue. 

Column  23  is  used  to  specify  the  symbolism  for  a value  greater  than  the 

maximum  specified  in  Elective  5. 

Column  24  is  used  to  specify  the  flag  point  symbolism  for  a high  value. 

Column  25  is  used  to  specify  background  symbolism.  Background  symbolism 

appears  outside  the  outline  of  the  study  area.  This  column  is 
normally  left  blank  to  indicate  that  no  symbolism  is  to  appear 
outside  of  the  study  area. 

On  cards  3,  4 and  5,  punch  in  the  columns  given  above  any  overprint  desired. 
Any  printer  characters  may  be  used.  Example  9 shows  a gray  scale  symbolism 
for  10  levels. 


Example  9: 

Gray  Scale  Symbolism  for  10  levels 
This  example  shows: 

(1)  A gray  scale  for  10  levels  of  symbolism  (columns  1-10). 

(?)  Flag  point  symbolism  for  the  10  levels  (columns  11-20). 

(3)  Blank  low  value  symbolism  and  flag  point  (columns  21-22). 

(4)  Blank  high  value  symbolism  and  flag  point  (columns  23-24). 

(5)  Blank  background  symbolism  (column  25). 

Column  111111111122222222 

123456789012345678901234567 

1 7 

2 . ,,+XOOOOO01 23456789 

3 . -X0000A 

4 -=XX 

5 V 

A 2 x 3 grid  cell  for  level  7 would  be  printed  on  the  GRID  map  as 
follows : 


card 

card 

card 

card 

card 


999 

979 

To  make  a gray  scale  for  less  than  10  levels,  the  level  symbols  can  be 
eliminated  in  the  following  order: 


for:  9 levels 

8 levels 
7 levels 
6 levels 
5 levels 


eliminate:  2 

2,  9 
2,  9,  8 
2,  9,  8,  3 
2,  9,  8,  3,  6 


The  flag  point  symbolism  should  be  adjusted  accordingly. 


j . Elective  8: 

Flag  Point  (1  card) 

The  flag  point  is  the  central  character  of  a grid  cell.  The  flag 
point  symbolism  specified  in  Elective  7 is  printed  at  the  flag  point.  If 
flag  point  symbolism  is  not  desired,  specify  1.0  in  field  1.  Standard  is 
special  symbolism  at  the  flag  point.  When  a mao  is  made  with  a 1 x 1 
character  grid,  the  flag  point  symbolism  is  automatically  suppressed. 

k.  Elective  9: 

Histogram  (1  card) 

This  elective  controls  the  printing  at  the  bottom  of  the  map. 
Specify  1.0  in  field  1 to  generate  a histogram  bar  chart  at  the  bottom  of 
the  map.  This  bar  chart  shows  the  frequency  of  grid  cells  in  each  level. 
Specify  1.0  in  field  2 to  suppress  numeric  information  which  is  included 
with  the  levels.  Standard  is  no  bar  chart  and  inclusion  of  numeric 
information. 


l . Electi ve  10 : 

Text  (3  to  32  cards) 

This  elective  is  used  to  print  up  to  30  lines  of  text  below  the 
map.  Standard  is  to  have  no  text.  However,  some  explanatory  text  is 
usually  desirable  for  interpreting  the  map. 

To  use  this  elective: 

(1)  On  card  1,  punch  the  identifying  elective  number  10  in  columns  4-5. 

(2)  On  not  more  than  30  other  cards,  to  be  inserted  between  the  first 
and  last,  punch  in  columns  1-72  any  supplementary  information  to 
be  used  as  the  text  for  the  map. 

(3)  On  the  last  card,  punch  ENDTEXT  in  columns  1-7. 

m.  Elective  13: 

Grid  Numbering  (1  card) 

This  elective  generates  row  and  column  numbers  on  all  four  sides  of 
the  grid  map  to  assist  the  user  in  locating  individual  cells  on  the  map. 

The  top  left  hand  cell  of  the  grid  is  called  the  Reference  Grid  Cell 
(RGC).  It  provides  the  coordinates  from  which  all  the  rows  and  columns  are 
numbered.  If  the  coordinates  of  the  RGC  are  not  specified,  the  program 
assumes  them  to  be: 


Column  = 1 

Row  = N (N  is  the  number  of  rows  specified  in  Elective  1) 

In  field  1 specify  1.0  for  grid  numbering.  In  field  2 specify  the  column 
number  of  the  RGC.  The  standard  for  Elective  13  is  no  numbering.  Example 
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10  shows  the  use  of  Elective  13  for  the  VAFB  irregular  outline. 

Example  10: 

Elective  13  for  the  VAFB  irregular  outline 

Column  4-5  18-20  28-30  36-40 

13  1.0  1.0  147.0 

When  Elective  13  is  specified  using  the  numbers  given  in  Example  10, 
the  row  numbers  will  correspond  to  the  row  numbers  given  in  the  C-1.2 
series  of  maps  in  the  Base  Master  Plan.  The  column  designators  have  been 
changed  to  correspond  to  the  alpha  designators  given  in  the  C-1.2  series  of 
maps  in  the  Base  Master  Plan.  Subroutine  FLAT0N  on  the  GRID  program  was 
modified  to  print  the  alpha  designators  when  Elective  13  is  specified. 
Therefore,  when  Elective  13  is  used,  the  alphanumeric  coordinates  of  a grid 
cell  correspond  to  the  grid  system  and  registration  described  in  Section  VII. 

Example  1 1 : 

An  Example  of  a Map  Package 

Example  2 showed  how  to  use  Flexin  to  produce  a map  which  displays  areas 
of  suitable  habitat  for  the  California  Legelss  Lizard  (Anniella  pulchra). 

This  example  shows  how  the  map  package  is  set  up  to  produce  the  map.  The 
following  cards  are  used  to  set  up  the  map  package: 

Card  1 : 

Column  1-3 
MAP 


Card  2: 

Col  umn 2-63 

ANNIELLA  PULCHRA  MICROHABITAT  / DISTRIBUTION  ON  VANDENBERG  AFB 


Card  3: 

B1  ank 

Card  4: 

Blank 

Card  5: 

Elective 

1 

Col umn 

5 

16-20 

27-30 

38-40 

48-50 

1 

146.0 

68.0 

1 .0 

1.0 

Card  6: 

Electi ve 

2 

Col umn 

5 

20 

2 2 


Card  7:  Elective  4 


I 


) 


Card  8:  Elective  5 

Column  5 11-14 

5 10.5 

Card  9:  Elective  6 

Column  5 11-13  21-23  31-33  41-43  51-53  61-63 

6 1.0  1.0  1.0  1.0  1.0  1.0 

Card  10:  Continuation  of  Elective  6 

Column  11-13  21-23  31-33  41-43 

TTo  ~77o  TTo  ' T7o 

Card  11:  Electi ve  7 

Col  umn 5 

7 

Card  12:  Continuation  of  Elective  7 

Column 1-20 

' „+XOOOOO.  123456789. 

Card  13:  Continuation  of  Elective  7 

Column  2-9 

. -X0000A 

Card  14:  Continuation  of  Elective  7 

Column 6-9 

-=XX 

Card  15:  Continuation  of  Elective  7 

Col umn 9 

V 

Card  16:  Elective  9 

Column  5 18-20 

9 1.0 

Card  17:  Elective  10 

Column  4-5^ 

' 10 
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Card  18: 


Continuation  of  Elective  10 


Col  umn 7-24 

THE  MAP  LEGEND  IS: 

Card  19:  Continuation  of  Elective  10 

Column  7-23 

LEVEL  1=  COASTAL  STRANDTT 

Card  20:  Continuation  of  Elective  10 

Column  7-60 

LEVEL  2=  COASTAL  SAGE  SCRUB-STABILIZED  DUNE  PHASE  ( 7 

Card  21:  Continuation  of  Elective  10 

Column  7-33 

LEVEL  3=  COASTAL  BLUFF  ( J 

Card  22:  Continuation  of  Elective  10 

Col  umn 7-31 _ 

LEVEL  4=  OAK  WOODLAND  ( ) 

Card  23:  Continuation  of  Elective  10 

Column  7-46 

LEVEL  5=  SUITABLE  SOIL  ASSOCIATIONS  fT 
Card  24:  Continuation  of  Elective  10 

Column  7-58 

LEVEL  6=  COASTAL  STRAND  / SUITABLE  SOIL'OVERLAP'TT' 

Card  25:  Continuation  of  Elective  10 

Column  7-71 

LEVEL  7=  COASTAL  SAGE  SCRUB- STABILIZED  DUNE  PHASE  / SUITABLE  SOIL 

Card  26:  Continuation  of  Elective  10 

Col umn  16-27 

OVERLAP  ( ) 

Card  27:  Continuation  of  Elective  10 

Column  7-57 

LEVEL  8=  COASTAL  BLUFF  / SUITABLE  SOIL  OVERLAP'd 
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Card  28:  Continuation  of  Elective  10 


Col  umn 7-56 

LEVEL  9=  OAK  WOODLAND  / SUITABLE  SOIL  OVERLAP  ("7 

Card  29:  Continuation  of  Elective  10 

Column 7-22 

BACKGR0UND=  

Card  30:  ENDTEXT 

Column  1-7 

ENDTEXT 

Card  31 : Electi ve  1 3 

Column  4-5  18-20  28-30  36-40 

13  1.0  1.0  147.0 

Card  32:  Last  Card  of  Map  Package 

Column  1-5 

99999 

6.  DATA  INPUT  PACKAGE 

The  GRID  program  has  two  input  options  for  the  data.  These  input 
options  optimize  the  use  of  the  program  under  two  different  types  of  opera- 
tions . 


a.  Data  Option  A. 

Option  A uses  GRID  as  an  independent  program  in  which  Subroutine 
Flexin  is  used  to: 

(1 ) Read  the  data. 

(2)  Perform  calculations  on  the  data  so  as  to  generate  the  value 
to  be  mapped. 

In  Option  A,  the  data  are  processed  one  cell  at  a time.  This  allows 
flexibility  in  the  organization  of  the  data.  The  Vandenberg  data  base  is 
designed  to  operate  using  Option  A. 

Each  time  a map  is  made  with  GRID  the  basic  data  file  must  be  read. 

The  Vandenberg  data  file  is  available  on  either  cards  or  magnetic  tape. 

Since  the  data  file  is  rather  large  (4646  cards),  it  is  more  efficient  to 
submit  the  data  file  on  magnetic  tape.  This  permits  the  user  to  rewind  the 
fi le  between  maps. 

Data  Option  A is  activated  by  specifying  a number  greater  than  zero  in 
field  1 of  Elective  2 (see  Map  Package).  This  number  is  transferred  to 
Subroutine  Flexin  as  the  value  of  IFORM.  I FORM  is  discussed  in  Section  VI3. 
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b.  Data  Option  8 

This  option  was  not  used  with  the  Vandenberg  data  base.  In  Option  B, 
the  data  used  to  create  the  map  are  transferred  to  the  GRID  program  as  a 
series  of  binary  arrays.  One  array  specifies  each  row  of  the  map.  The  program 
will  expect  one  real  value  in  the  array  for  each  cell  in  a row. 

Option  B is  different  from  Option  A in  that  it  uses  GRID  as  a final  job 
step  in  a series  of  job  steps  to  produce  the  map. 

Option  B is  used  automatically  if  Elective  2 is  not  specified.  It  is 
also  used  if  zero  is  specified  in  field  1 of  Elective  2. 

7.  COMPUTER  SUBMISSIONS 

After  the  packages  have  been  prepared  they  must  be  placed  in  the  correct 
order, together  with  the  control  cards  needed  for  submission  to  the  computer. 

The  following  procedure  is  used  to  operate  the  GRID  program  with  the  VAFB  data 
base  on  the  Burroughs  3500  computer. 

a.  To  compile  the  GRID  program  and  load  to  the  program  library  use  the 
following  control  cards: 

(1 ) Column 1-37 

1 

2C0MPILE  SNFT00  XFORTN  LIB  DATA  CARDS 
3 

(2)  Col  umn 1-5 

11 

3 LST1 
8 

GRID  program  (including  Subroutine  Flexin) 

(3)  Column  1-4 

1 

2END 

3 

b.  To  execute  the  GRID  program  and  produce  the  map  use  the  following  control 
cards : 

(1)  Column  1-21 

1 

2EX  SNFT00  DATA  FILE5 
3 

Data  on  which  the  program  operates 

(2)  Column  1-4 

1 

2END 

3 
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c.  Data  on  which  the  program  operates: 

The  data  on  which  the  program  operates  (see  above)  consist  of  the 
Irregular  Outline  Package,  the  Map  Package,  and  the  Data  Input.  These  packages 
must  be  in  the  correct  order  in  the  input  card  deck. 

(1)  The  order  for  submitting  the  deck  when  the  Data  Input  is  read  from  cards 
is : 

Control  Card 

Irregular  Outline  Package 

Map  Package(s) 

Data  Input 

Control  Card 

The  end  of  the  Data  Input  is  signaled  by  a card  with  END  punched  in 
columns  1-3.  This  card  immediately  follows  the  last  data  card. 

The  Irregular  Outline  Package  must  precede  the  Map  Package  to  which  it 
refers.  Once  the  irregular  outline  has  been  specified  it  will  be  used  for 
every  Map  Package  in  the  computer  submission. 

Each  time  the  program  reads  a Map  Package  it  will  attempt  to  make  a 
map.  There  is  no  limit  to  the  number  of  Map  Packages  in  any  one  submission. 

If  the  Data  Input  is  read  from  cards,  it  must  immediately  follow  the 
Map  Package  to  which  it  refers. 

(2)  The  order  for  submitting  the  deck  when  the  Data  Input  is  read  from 
magnetic  tape  is: 

Control  Card 

Irregular  Outline  Package 

Map  Package(s) 

Control  Card 

The  end  of  the  Data  Input  is  signaled  by  a card  with  END  punched  in 
columns  1-3  immediately  following  the  last  card  in  the  Map  Package.  The  READ 
statement  in  Subroutine  Flexin  must  be  changed  to  indicate  that  the  data  is 
being  read  from  tape  instead  of  cards. 

8.  SEARCH/COUNT  PROGRAM 

This  program  can  be  used  in  conjunction  with  GRID  or  independently  as  a 
feature  of  the  EPS  developed  for  VAFB.  The  Search/Count  program  will  locate 
and  print  out  data  base  variables  for  selected  grid  cells.  By  entering  the 
alphanumeric  coordinates  of  selected  grid  cells,  the  program  will  print  out 
the  data  base  variables  for  each  grid  cell,  count  the  frequency  of  each  subvari- 
able, and  convert  the  frequency  to  estimated  acres.  The  Search/Count  program 
is  useful  for  studies  such  as  computing  the  areas  of  various  vegetation  types 
that  would  be  affected  by  alternate  placements  of  a runway-launch  complex. 

The  Search/Count  program,  like  GRID,  is  written  in  FORTRAN  IV  and  requires 
the  user  to  have  a basic  knowledge  of  computer  programming.  A listing  of  the 


Search/Count  program,  as  used  in  a case  study,  is  given  in  Appendix  G.  The 
user  should  refer  to  this  listing  as  a guide  in  setting  up  the  program. 

Certain  statements  within  the  Search/Count  program  must  be  modified  to  fit 
the  user's  needs.  These  modifications  are  given  in  the  following  sections. 

a.  Title 

A title  identifying  the  particular  use  of  the  Search/Count  program  should 
be  included  at  the  beginning  of  the  program  using  comment  cards. 

b.  REAL*8  Statement 

This  statement  must  be  modified  to  reflect  the  exact  number  of  23-acre 
grid  cells  to  be  counted  in  the  program.  To  count  43  grid  cells,  the  statement 
would  be: 

REAL*8  C(43) , CX 

c.  Size  of  Counting  Arrays 

The  size  of  arrays  used  in  the  counting  program  need  to  be  adjusted  for 
the  number  of  cells  counted.  For  counting  up  to  100  grid  cells,  the  array 
sizes  used  for  the  program  in  Appendix  G can  be  used. 

d.  Data  Statement 

This  statement  is  used  to  enter  the  alphanumeric  coordinates  of  the 
selected  grid  cells.  The  general  form  is: 

DATA  C/ ' xxxxx ' , 1 xxxxx 1 , ' xxxxx ' / 

Columns  7-72  are  used  for  the  Data  Statement.  Since  only  7 grid  cells  can  be 
entered  on  the  first  card,  continuation  cards  are  required.  The  general  form 
for  the  alphanumeric  coordinate  of  a grid  cell  is:  'xxxxx'. 
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SECTION  VII 

QUANTITATIVE  ECOLOGICAL  DATA  BASE 
1.  GRID  SYSTEM  AND  REGISTRATION 

The  grid  system  and  registration  used  for  the  study  of  VAFB  is  based  on 
the  California  Coordinate  System.  The  basic  units  of  this  system  are  feet. 

The  California  Coordinate  System  was  chosen  for  the  Vandenberg  study  because 
of  its  flexibility.  Most  maps  used  in  relation  to  the  project  were  registered 
in  California  coordinate  units.  This  greatly  facilitated  the  transfer  of 
information  into  the  data  base. 

A basic  grid  cell  size  of  1000  by  1000  feet  was  adopted  to  register  all 
information  collected  for  the  computerized  data  base.  This  represents  a 
subdivision  into  24  equal  parts  of  the  standard  4000  by  6000-foot  sections  of 
the  California  Coordinate  System.  The  problem  in  determining  grid  cell  size 
is  that  it  should  not  be  too  small  to  be  unfeasible  for  storage  and  economy 
and  it  should  be  small  enough  for  the  required  analysis  purposes.  The  1000- 
foot  grid  cell  has  been  used  effectively  in  previous  land  use  studies  (Refer- 
ence 3).  It  was  determined  to  be  suitable  for  this  study  only  after  thorough 
evaluation  and  discussion  of  the  scope  of  the  project. 

The  1000-foot  grid  cell  encompasses  an  area  of  approximately  23  acres. 
There  are  4646  of  the  1000-foot  grid  cells  enclosing  the  approximate  boundaries 
of  VAFB.  If  a smaller  grid  cell  size  is  required  for  registering  data,  the 
1000-foot  cell  can  be  subdivided  into  nine  2.55-acre  grid  cells  measuring 
333.33  feet  on  each  side  (see  Figure  7).  For  example,  sensitive  areas  that 
need  more  detailed  analysis  could  be  redefined  from  the  larger  23-acre  units 
into  the  smaller  2.55-acre  units. 
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Figure  7.  Grid  Cell  Scales 

An  alphanumeric  registration  system  was  adopted  to  locate  a particular 
grid  cell  in  the  Vandenberg  area.  Each  1000-foot  vertical  line  is  coded  with 
a two-letter  designator.  The  1000-foot  horizontal  lines  are  coded  with  number 
designators.  See  C-1.2  series  of  maps  in  the  Base  Master  Plan  (Reference  4). 

A series  of  letters  and  numbers  will  then  specify  a particular  grid  cell  in 
the  Vandenberg  area.  The  cell  is  identified  by  a point  in  the  southwest 
corner  of  the  grid  cell.  For  example,  RA055  represents  the  cell  where  the 
vertical  line  RA  intersects  the  horizontal  line  055.  Figure  8 illustrates  this. 


057  j- 

056  ^ 

055  : 

054  ! i 

RA  SA  TA 

Figure  8.  Grid  Ceil  Coordinates 

If  the  coordinate  of  the  grid  cell  specified  in  Figure  7 is  RA055,  the 
coordinate  of  the  fourth  subcell  is  specified  as  RA055.4. 

2.  SOIL  DATA  BASE 

The  U.S.  Department  of  Agriculture  soil  surveys  of  Northern  Santa  Barbara 
area  and  Santa  Barbara  area  (References  9,  10)  were  used  to  classify  the  soils 
occurring  within  VAFB.  The  soils  were  classified  and  coded  according  to: 

a.  soil  series  and  phase  or 

b.  land  type 

c.  slope  of  the  land 

Soils  with  similar  profiles  make  up  a soil  series.  A soil  series  is 
divided  into  phases  based  on  the  texture  of  the  surface  soil  and  slope.  A land 
type  is  soil  material  that  cannot  be  classified  into  a soil  series  because  it 
is  too  rocky,  shallow,  or  severely  eroded.  A land  type  is  given  a descriptive 
name  in  contrast  to  a soil  series  which  is  named  for  the  town  or  geographic 
feature  near  the  place  where  the  soil  of  that  series  was  first  observed  and 
mapped. 

A transparent  grid  overlay  was  used  with  the  soil  survey  maps  to  code 
the  soils.  Each  grid  cell  on  the  overlay  represented  the  standard  23-acre 
grid  cell.  The  predominant  soil  series  or  land  type  occurring  on  each  23-acre 
grid  cell  was  used  for  the  overall  interpretation  of  that  cell.  Forty-three 
different  soil  series  and  13  land  types  were  classified  on  VAFB  using  this 
coding  method. 

The  combination  of  the  soil  series/phase  or  land  type  with  the  slope 
characteristic  can  be  used  to  determine  the  erosion  potential  on  the  soil. 
Similarly,  it  can  be  used  to  predict  the  potential  occurrence  of  a plant  or 
animal  according  to  its  substrate  preference.  One  hundred  eighty-three 
different  soil/slope  categories  were  coded  on  VAFB  using  the  above  system. 

Table  17  lists  the  soil/slope  categories  in  the  data  base. 

The  soil/slope  type  for  each  grid  cell  was  coded  on  computer  data  cards 
in  integer  form  using  the  code  numbers  listed  in  Table  17.  Twenty-five  fields 
of  three  columns  each  were  defined  in  columns  1-75  of  each  data  card. 

Columns  76-80  were  assigned  an  identification  number  for  ordering  purposes. 
Therefore,  the  soil/slope  types  for  25  grid  cells  were  coded  on  each  computer 
data  card. 


1 


TABLE  17.  SOIL/SLOPE  CODES  AND  CATEGORIES  FOR  VANDENBERG  AFB 


' Code 
(Number 

Symbol 

Name 

; 0 

Ocean 

2 

AgC 

Agueda  silty  clay  loam,  0 to  2%  slopes 

| 4 

ArD 

Arnold  sand,  5 to  15%  slopes 

5 

ArF 

Arnold  sand,  15  to  45%  slopes 

6 

ArF3 

Arnold  sand,  9 to  45%  slopes,  severely  eroded 

16 

Bd 

Bayshore  loam,  drained 

! 18 

Bg 

Bayshore  silty  clay  loam 

, 23 

BnB2 

Betteravia  loamy  sand,  dark  variant,  0 to  5%  slopes,  eroded 

| 24 

BnD2 

Betteravia  loamy  sand,  dark  variant,  5 to  15%  slopes,  eroded 

! 25 

BoA 

Botella  loam,  0 to  2%  slopes 

1 27 

BoC 

Botella  loam,  2 to  9%  slopes 

28 

BoD2 

Botella  loam,  2 to  15%  slopes,  eroded 

29 

BsA 

Botella  loam,  slightly  wet,  0 to  2%  slopes 

30 

BtA 

Botella  clay  loam,  0 to  2%  slopes 

31 

BtA2 

Botella  clay  loam,  0 to  2%  slopes,  eroded 

32 

BtC 

Botella  clay  loam,  2 to  9%  slopes 

33 

BtD2 

Botella  clay  loam,  2 to  15%  slopes,  eroded 

j 34 

BwA 

Botella  clay  loam,  wet,  0 to  2%  slopes 

1 35 

Ca 

Camarillo  sandy  loam 

j 37 

Cc 

Camarillo  very  fine  sandy  loam 

38 

Cd 

Camarillo  silty  clay  loam 

i 43 

ChD 

Chamise  shaly  loam,  9 to  15%  slopes 

44 

ChF 

Chamise  shaly  loam,  15  to  45%  slopes 

46 

ChG2 

Chamise  shaly  loam,  30  to  75%  slopes,  eroded 

48 

CmF 

Cl imara-Toomes  complex,  15  to  45%  slopes 

I 49 

CnB 

Coastal  beaches 

52 

CrF 

Contra-Costa-Lodo  loams,  30  to  45%  slopes 

55 

CtA 

Corral itos  sand,  0 to  2%  slopes 

56 

CtD 

Corralitos  sand,  2 to  15%  slopes 

57 

CtD2 

Corralitos  sand,  9 to  15%  slopes,  eroded 

58 

CuA 

Corralitos  loamy  sand,  0 to  2%  slopes 

59 

CuC 

Corralitos  loamy  sand,  2 to  9%  slopes 

60 

CuD 

Corralitos  loamy  sand,  9 to  15%  slopes 

62 

CwE 

Crow  Hill  loam,  15  to  30%  slopes 

63 

CwF 

Crow  Hill  loam,  30  to  45%  slopes 

64 

CwG 

Crow  Hill  loam,  45  to  75%  slopes 

65 

CwG3 

Crow  Hill  loam,  15  to  75%  slopes,  severely  eroded 

66 

DaD 

Diablo  silty  clay,  9 to  15%  slopes 

67 

DaE 

Diablo  silty  clay,  15  to  30.  slopes 

68 

DaF 

Diablo  silty  clay,  30  to  45%  slopes 

71 

DuE 

Dune  land 

74 

EdA2 

Elder  sandy  loam,  0 to  2%  slopes,  eroded 

75 

EdC2 

Elder  sandy  loam,  2 to  9%  slopes,  eroded 

76 

EdD2 

Elder  sandy  loam,  9 to  15%  slopes,  eroded 

78 

EmC 

Elder  loam,  2 to  9%  slopes 

79 

EnA2 

Elder  shaly  loam,  0 to  2%  slopes,  eroded 

80 

EnC2 

Elder  shaly  loam,  2 to  9%  slopes,  eroded  | 
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TABLE  17.  SOIL/SLOPE  CODES  AND  CATEGORIES  FOR  VANDENBERG  AFB  (Continued) 


I Code 
Number 

Symbol 

Name 

81 

EnD2 

Elder  shaly  loam,  9 to  15%  slopes,  eroded 

89 

GmG 

Gaviota  sandy  loam,  30  to  75%  slopes 

91 

GsE 

Gazos  clay  loam,  15  to  30%  slopes 

92 

GsF 

Gazos  clay  loam,  30  to  45%  slopes 

93 

GsG 

Gazos  clay  loam,  45  to  75%  slopes 

94 

GuE 

Gullied  land 

99 

LaF 

Landsl ides 

101 

LcE 

Linne  clay  loam,  15  to  30%  slopes 

102 

LcF 

Linne  clay  loam,  30  to  45%  slopes 

103 

LcG 

Linne  clay  loam,  45  to  75%  slopes 

106 

LmG 

Lopez  shaly  clay  loam,  15  to  75%  slopes 

107 

LoE 

Los  Osos  clay  loam,  15  to  30%  slopes 

108 

LoG 

Los  Osos  clay  loam,  30  to  75%  slopes 

109 

LsE 

Los  Osos -San  Benito  clay  loams,  15  to  30%  slopes 

110 

LsF 

Los  Osos-S an  Benito  clay  loams,  30  to  45%  slopes 

112 

MaA 

Marina  sand,  0 to  2%  slopes 

113 

MaC 

Marina  sand,  2 to  9%  slopes 

114 

MaE 

Marina  sand,  9 to  30%  slopes 

115 

MaE3 

Marina  sand,  9 to  30%  slopes,  severely  eroded 

116 

Mh 

Marsh 

118 

MnA 

Metz  loamy  sand,  0 to  2%  slopes 

119 

MnC 

Metz  loamy  sand,  2 to  9%  slopes 

123 

Mr 

Mocho  sandy  loam,  overflow 

126 

Mu 

Mocho  fine  sandy  loam 

127 

Mv 

Mocho  loam 

129 

Mx 

Mocho  silty  clay  loam 

130 

MyG 

Montara  rocky  clay  loan,  30  to  75%  slopes 

131 

NrB 

Narlon  sand,  0 to  5%  slopes 

132 

NsA 

Narlon  loamy  sand,  0 to  2%  slopes 

133 

NsC 

Narlon  loamy  sand,  2 to  9%  slopes 

134 

NsD 

Narlon  loamy  sand,  9 to  15%  slopes 

136 

NvC 

Narlon  sand,  hardpan  variant,  2 to  9 l slopes 

137 

OcA 

Oceano  sand,  0 to  2%  slopes 

138 

OcD 

Oceano  sand,  2 to  15%  slopes 

159 

Rs 

Ri verwash 

160 

RuG 

Rough  broken  land 

161 

SaA 

Salinas  loam,  0 to  2%  slopes 

162 

SaC 

Salinas  loam,  2 to  9;  slopes 

166 

SeD 

Salinas  and  Sorrento  loams,  9 to  15  slopes 

167 

SfD 

San  Andreas-Tierra  complex,  5 to  15%  slopes 

168 

SfE 

San  Andreas-Tierra  complex,  15  to  30%  slopes 

169 

SfF3 

San  Andreas-Tierra  complex,  9 to  45  slopes,  severely  eroded 

170 

SfG 

San  Andreas-Tierra  complex,  30  to  75  slopes 

172 

SgG 

San  Beni to-Di ablo  complex,  45  to  75  slopes 

173 

Sh 

Sandy  al luvial  land 

174 

Sk 

Sandy  alluvial  land,  wet 

175 

SmD 

Santa  Lucia  shaly  clay  loam,  9 to  15%  slopes 

176 

SmE 

Santa  Lucia  shaly  clay  loam,  15  to  30%  slopes 
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TABLE  17.  SOIL/SLOPE  CODES  AND  CATEGORIES  FOR  VANDENBERG  AFB  (Continued) 


Code 
IN  umber 

Symbol 

Name 

177 

SmF 

Santa  Lucia  shaly  clay  loam,  30  to  45%  slopes 

179 

SmG 

Santa  Lucia  shaly  clay  loam,  45  to  75%  slopes 

184 

SpG 

Sedimentary  rock  land 

185 

SrE 

Shedd  silty  clay  loam,  15  to  30%  slopes 

186 

SrF 

Shedd  silty  clay  loam,  30  to  45%  slopes 

187 

SrG 

Shedd  silty  clay  loam,  45  to  75%  slopes 

188 

SrG3 

Shedd  silty  clay  loam,  30  to  75%  slopes,  severely  eroded 

190 

SsF 

Shedd  silty  clay  loam,  diatomaceous  variant,  30  to  45%  slopes| 

196 

SvC 

Sorrento  loam,  2 to  9%  slopes 

204 

Szw 

Swamp 

205 

TaA 

Tangair  sand,  0 to  2%  slopes 

206 

TaC 

Tangair  sand,  2 to  9%  slopes 

207 

TcG 

Terrace  escarpments,  sandy 

208 

TdF 

Terrace  escarpments,  loamy 

211 

TmE 

Tierra  loamy  sand,  9 to  30%  slopes 

212 

TnC 

Tierra  sandy  loam,  2 to  9%  slopes 

213 

TnD2 

Tierra  sandy  loam,  9 to  15%  slopes,  eroded 

214 

TnE2 

Tierra  sandy  loam,  15  to  30%  slopes,  eroded 

215 

TrC 

Tierra  loam,  2 to  9%  slopes 

216 

TrD 

Tierra  loam,  9 to  15%  slopes 

217 

TrE2 

Tierra  loam,  15  to  30%  slopes,  eroded 

218 

TrE3 

Tierra  loam,  5 to  30%  slopes,  severely  eroded 

219 

TsF 

Tierra  clay  loam,  15  to  45%  slopes 

220 

TxG 

Toomes-Cl imara  complex,  30  to  75  slopes 

226 

Td 

Tangair  sand,  16  to  30%  slopes 

227 

Th 

Tangair  sand,  9 to  15%  slopes,  severely  eroded 

228 

Tf 

Tangair  sand,  9 to  15%  slopes 

229 

Su 

Santa  Lucia,  stony  soils,  undifferentiated,  31  + slopes 

230 

Lp 

Los  Osos , stony  soils,  undifferentiated,  31  + slopes 

231 

Te 

Tangair  sand,  16  to  30.  slopes,  moderately  eroded 

232 

Ta 

Tangair  loamy  sand,  16  to  30  slopes 

233 

Tb 

Tangair  loamy  sand,  9 to  15  slopes 

234 

Tc 

Tangair  loamy  sand,  9 to  15  slopes,  moderate’/  -ended 

235 

At 

Arguello  shaly  loam,  9 to  15  slopes 

236 

B1 

Baywood  loamy  sand,  gently  sloping,  3 to  8 

237 

Bn 

Baywood  loamy  sand,  over  Watsonville  soils,  lently  sloping. 

3 to  8%  slopes 

238 

Sq 

Santa  Lucia  shaly  loam,  16  to  30  slopes 

239 

Sr 

Santa  Lucia  shaly  loam,  31  to  45  slopes 

240 

Lr 

Los  Trancos  stony  loam,  16  to  45  slopes 

241 

Lf 

Los  Osos  clay,  31  to  45%  slopes 

242 

Yd 

Yolo  loam,  9 to  15%  slopes 

243 

Le 

Los  Osos  clay,  16  to  30%  slopes 

244 

Ja 

Jalama  shaly  sandy  loam,  3 to  15  slopes 

245 

Wh 

Watsonville  loam,  3 to  8 slopes 

246 

As 

Arguello  shaly  loam,  3 to  8%  slopes 

247 

Bp 

Baywood  loamy  sand,  9 to  15%  slopes 

248 

Ec 

Elder  shaly  clay  loam,  3 to  8%  slopes 
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TABLE  17.  SOIL/SLOPE  CODES  AND  CATEGORIES  FOR  VANDENBERG  AFB  (Continued) 


Code 

Number 

Symbol 

Name 

249 

St 

Santa  Lucia  stony  clay  loam,  16  to  30%  slopes 

250 

Bk 

Baywood  loamy  fine  sand,  over  Watsonville  soils,  3 to  8% 
slopes 

251 

Ss 

Santa  Lucia  shaly  loam,  46%+  slopes 

252 

Nd 

Nacimiento  clay,  31  to  45%  slopes 

253 

Cn 

Climax  clay  (adobe),  31  to  45%  slopes 

254 

Cm 

Climax  clay  (adobe),  16  to  30%  slopes 

255 

Gk 

Gaviota  stony  soils,  undi fferenti ated , 31%+  slopes 

256 

Ne 

Nacimiento  clay,  31  to  45%  slopes,  moderately  eroded 

257 

Rh 

Rough  gullied  land,  Los  Osos  soil  material 

258 

Yg 

Yolo  loam,  0 to  2%  slopes 

259 

Ye 

Yolo  loam,  3 to  8%  slopes 

260 

Na 

Nacimiento  clay,  16  to  30%  slopes 

261 

Ct 

Crow  Hill  loam,  9 to  15%  slopes 

263 

Sa 

San  Andreas  fine  sandy  loam,  moderately  eroded,  9 to  15% 
slopes 

264 

Tp 

Tierra  fine  sandy  loam,  moderately  eroded,  31  to  45%  slopes 

265 

(Sh) 

San  Andreas  stony  soils,  undifferentiated,  46%+  slopes 

267 

Zv 

Zaca  stony  soils,  undifferentiated,  31%+  slopes 

270 

Zm 

Zaca  shaly  clay  loam,  16  to  30%  slopes 

271 

Zk 

Zaca  shaly  clay  loam,  31  to  45%  slopes 

272 

Zg 

Zaca  shaly  clay  loam,  16  to  30%  slopes 

273 

Za 

Zaca  clay,  16  to  30%  slopes 

274 

Jc 

Jalama  stony  soils,  undifferentiated,  16  to  45%  slopes 

275 

Jb 

Jalama  shaly,  sandy  loam,  16  to  30%  slopes 

276 

Zs 

Zaca,  shaly  clay  loam,  31  to  45%  slopes 

277 • 

Kb 

Kitchen  middens,  over  impermeable  soil  material 

278 

Zc 

Zaca  clay,  9 to  15%  slopes 

279 

Zd 

Zaca  clay,  31  to  45%  slopes 

280 

N1 

Nacimiento  stony  soils,  undifferentiated,  46%+  slopes 

281 

S 

San  Andreas  fine  sandy  loam,  16  to  30%  slopes 

282 

Ad 

Agueda,  gravelly  clay  loam,  3 to  8%  slopes 

283 

Ae 

Agueda,  gravelly  clay  loam,  9 to  15%  slopes 

284 

MN 

Montezuma  clay  (adobe),  3 to  8%  slopes 

285 

Nf 

Nacimiento  clay,  46%+  slopes 

286 

Zp 

Zaca  shaly  clay  loam,  9 to  15%  slopes 

287 

Ac 

Agueda  clay  loam,  9 to  15%  slopes 

288 

Zn 

Zaca  shaly  clay  loam,  moderately  eroded,  16  to  30.  slopes 

289 

Z1 

Zaca  nonstony  soils,  undifferentiated,  46.+  slopes 

290 

Wn 

Watsonville  loam,  9 to  15%  slopes 

291 

MV 

Montezuma  clay  loam,  3 to  8%  slopes 

292 

Tg 

Tangair  sand,  moderately  eroded,  9 to  15  slopes 
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3.  EXPOSURE  DATA  BASE 

The  exposure  or  aspect  of  the  land  was  coded  for  each  23-acre  grid  cell. 

A transparent  grid  overlay  was  used  with  the  C series  of  maps  in  the  Base 
Master  Plan  (Reference  4)  to  code  the  exposure  cell  by  cell.  If  multiple 
exposures  occurred  within  a grid  cell,  the  predominant  exposure  was  used  to 
describe  the  whole  cell.  Each  cell  was  assigned  an  integer  exposure  code  for 
processing  exposure  into  the  data  base  for  the  GRID  computer  program. 

A series  of  25  different  exposure  codes  were  used  to  evaluate  the  terrain. 
Codes  1-9  were  used  if  a cell  contained  only  a predominant  single  exposure 
(example:  south  facing).  Codes  10-17  were  used  for  multiple  exposures  in  a 

cell  caused  by  drainage  areas.  For  example,  a north  sloping  drainage  provides 
both  east  and  west  exposure.  Codes  18-25  were  used  for  multiple  exposures  in 
a cell  caused  by  a ridge  line.  For  example,  a ridge  sloping  to  the  northeast 
has  northwest  and  southeast  exposure.  Table  18  lists  the  codes  and  categories 
used  for  describing  exposure. 

The  exposure  data  were  coded  on  computer  data  cards  in  two-column  fields. 
Thirty-seven  fields  of  two  columns  each  were  defined  in  columns  1-74.  Columns 

75- 80  were  assigned  a card  identification  number.  Using  this  coding  system, 

37  grid  cells  were  contained  on  each  data  card. 

4.  ELEVATION  DATA  BASE 

The  elevation  of  the  land  was  coded  for  each  grid  cell.  Elevation  was 
coded  using  the  same  basic  procedure  described  for  coding  exposure.  The 
highest  elevation  in  a grid  cell  determined  the  overall  elevation  assigned  to 
that  cell.  Elevation  was  coded  in  200-foot  intervals.  Using  this  system, 

10  different  codes  describe  the  elevation  from  sea  level  to  the  highest  point 
on  VAFB.  All  areas  1800  feet  and  above  were  combined  into  one  elevation  class. 
Table  19  lists  the  elevation  codes  and  categories  used  in  making  the  data  base. 

The  elevation  data  were  coded  on  computer  data  cards  in  one-column  fields. 
Seventy-five  fields  of  one-column  each  were  defined  in  columns  1-75.  Columns 

76- 80  were  assigned  a card  identification  number. 

The  exposure  and  elevation  data  can  be  used  for  a variety  of  purposes. 
Examples  include:  hydrologic  studies,  meteorological  studies,  and  studies 

such  as  potential  effects  of  exhaust  clouds  from  rocket  launches  on  land  areas 
below  the  inversion  layer. 

5.  VEGETATION  DATA  BASE 

The  vegetation  data  base  was  developed  from  the  vegetation  maps  prepared 
for  VAFB.  The  maps  are  transparent  overlays  to  the  Base  Master  Plan,  map 
series  C-l  (66  sheets,  1 inch:800  feet  scale).  After  the  maps  were  drawn  the 
vegetation  was  coded  for  the  computerized  data  base.  A transparent  grid 
overlay  was  used  with  the  vegetation  maps  to  code  the  data.  The  vegetation 
data  were  coded  at  the  2.55-acre  subcell  level.  Each  23-acre  grid  ce1!  was 
divided  into  nine  subcells  using  a transparent  grid  overlay  (see  Figure  7). 

The  dominant  type  of  vegetation  of  each  of  the  subcells  was  coded  onto 
data  forms.  In  many  cases,  the  determination  of  the  dominant  type  was  a 
complicated  matter  because  the  scale  of  the  maps  was  such  that  one  2.55-acre 
cell  could  contain  several  different  vegetation  types.  If  one  vegetation  type 
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TABLE  18.  EXPOSURE  CODES  AND  CATEGORIES  FOR  VANDENBERG  AFB. 


Code 
IN  umber 

Exposure 

1 Undefinable  direction,  less  than  10  degrees 

2 North 

3 Northwest 

! 4 West 

I 5 Southwest 

6 South 

7 Southeast 

8 East 

9 Northeast 

Codes  10-17  describe  multiple  exposures  in  a cell  due 
Code 

Number  Exposure  Drainage 

of  slope 

to  drainage  areas. 

1 

! 10 

East  and  West 

North 

n 

East  and  West 

South 

: 12 

Northwest  and  Southeast 

Northeast 

i 13 

Northwest  and  Southeast 

Southwest 

1 14 

North  and  South 

East 

15 

North  and  South 

West 

! 16 

Northeast  and  Southwest 

Northwest 

17 

Northeast  and  Southwest 

Southeast 

Codes  18-25  describe  multiple  exposures 

in  a cell 

due 

to 

ridge  lines. 

Code 

iNumber 

Exposure 

Ri  dqe 

18 

East  and  West 

North 

19 

East  and  West 

South 

20 

Northwest  and  Southeast 

Northeast 

21 

Northwest  and  Southeast 

Southwest 

22 

North  and  South 

East 

1 

23 

North  and  South 

West 

24 

Northeast  and  Southwest 

Northwest 

25 



Northeast  and  Southwest 

Southeast 

1 
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TABLE  19.  ELEVATION  CODES  AND  CATEGORIES 
FOR  VANDENBERG  AFB . 


Code 

Number 

Elevation  (Feet) 

0 

0 - 200 

1 

200  - 400 

! 

i 2 

400  - 600 

3 

600  - 800 

4 

800  - 1000 

5 

1000  - 1200 

6 

1200  - 1400 

7 

1400  - 1600 

8 

1600  - 1800 

9 

1800  and  above 

It 


* 
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clearly  covered  the  subcell  , it  was  designated  as  the  type  for  the  entire 
subcell . If  one  vegetation  type  was  not  dominant  in  area  within  the  subcell, 
relations  of  vegetation  types  to  other  types  within  and  outside  the  subcell 
were  taken  into  account. 

The  vegetation  types  were  given  an  integer  code  for  processing  into  the 
data  base.  Table  20  lists  the  codes  and  categories. 

One  computer  data  card  was  used  for  each  grid  cell  of  the  vegetation  data 
base.  Table  21  lists  the  fields  and  columns  used  for  each  data  card.  The 
field  used  for  describing  optional  information  for  each  subcell  was  used  to 
code  subcells  containing  fire  breaks.  The  code  number  1 in  the  optional  field 
indicates  a fire  break  occurs  in  that  subcell.  Columns  45-68  on  each  grid 
cell  data  card  are  blank  and  can  be  used  to  code  additional  information  about 
each  gri d cell. 

The  GRID  program  processes  the  vegetation  data  for  a GRID  map  on  a 23-acre 
cell  size.  Therefore,  Subroutine  Flexin  must  be  modified  to  assign  an  overall 
vegetation  type  to  a grid  cell  on  the  basis  of  the  vegetation  types  assigned 
to  the  9 subcells.  Examples  12  and  13  show  two  different  modifications  of 
Flexin  for  assigning  an  overall  vegetation  type  to  a grid  cell.  Other  modifi- 
cations of  Flexin  are  possible  depending  upon  the  user's  needs  and  knowledge 
of  FORTRAN. 

Example  12: 

Assigning  Vegetation  Type  by  the  Center  Subcell 

The  FORTRAN  statements  given  in  this  example  of  Subroutine  Flexin  direct 

the  computer  to  assign  an  overall  vegetation  type  to  a grid  cell  based  on  the 

vegetation  type  of  the  center  subcell  (subcell  5).  It  can  be  demonstrated 
statistically  that  this  method  will  give  an  unbiased  estimate  of  vegetation 
coverage  for  all  the  grid  cells.  The  FORTRAN  statements  are: 

Explanatory  Notes: 

SUBROUTINE  FLEX  IN ( IFORM,  T,  FIRST) 

READ ( 1 1 , 11)  V 

11  FORMAT (36X  ,F2.0)  a. 

T=5.  b. 

I F ( V . EQ .12.  ) T=1 . 

I F ( V . EQ . 8 . ) T=2. 

I F ( V . EQ .11.)  T=3. 

IF (V . EQ . 3. ) T=4. 

I F ( V . EQ . 31 . ) T=4. 

RETURN 

END 

Explanatory  Notes  for  Example  12: 

a.  This  FORMAT  statement  and  the  READ  statement  direct  the  computer  to  read 
the  vegetation  type  for  subcell  5. 

b.  The  remainder  of  the  FORTRAN  statements  in  the  subroutine  direct  the  computer 
to  produce  a vegetation  map  of  5 levels,  based  on  the  instructions  given 

in  the  Map  Package.  In  this  case  the  Map  Package  would  be  set  up  to  produce 
a GRID  map  showing  4 types  of  vegetation  and  the  remainder  as  background. 
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TABLE  20.  VEGETATION  CODES  AND  CATEGORIES  FOR  VANDENBERG  AFB . 


Code 

Number  Category 

1 

Bishop  pine  forest 

72 

Bishop  pine  forest 

- sparse  phase 

2 

Tanbark  oak  forest 

3 

Foothill  woodland 

31 

Foothill  woodland  - 

dense  phase 

4 

Riparian  woodland 

42 

Riparian  woodland  - 

sparse  phase 

5 

Chaparral 

52 

Chaparral  - sparse 

phase 

6 

Coastal  sage  scrub 

- normal  phase 

62 

Coastal  sage  scrub 

- sparse  phase 

7 

Coastal  sage  scrub 

- Salvia  leucophylla  phase 

8 

Coastal  sage  scrub 

- stabilized  dune  phase 

9 

Wet  soil  scrub 

10 

Fluckleberry  scrub 

11 

Coastal  bluff 

12 

Coastal  strand 

13 

Coastal  salt  marsh 

14 

Freshwater  marsh 

16 

Grassland  - annual 

17 

Miscellaneous  native  herb  communities 

18 

Ruderal  vegetation 

19 

Planted  trees 

20 

Agricultural  plantings 

21 

Non-agricul tural  plantings 

22 

Freshwater 

23 

Man-made  facilities 

and  cantonement 

24 

Disked  areas 

25 

Naturally  bare  soil 

26 

Acer  negundo  stands 

99 

Land  not  within  the 

base  boundary 

00 

Ocean 
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TABLE  21.  VEGETATION  DATA  CARD  FIELDS. 


Col umn 

Field 

2-6 

Alphanumeric  coordinates  for  the  grid  cell 

10-11 

Subcel 1 

1 

vegetation  code 

12 

Subcell 

1 

optional  field 

14-15 

Subcel 1 

2 

vegetation  code 

16 

Subcel 1 

2 

optional  field 

18-19 

Subcel 1 

3 

vegetation  code 

20 

Subcel 1 

3 

optional  field 

22-23 

Subcel 1 

4 

vegetation  code 

24 

Subcel 1 

4 

optional  field 

26-27 

Subcel 1 

5 

vegetation  code 

28 

Subcel  1 

5 

optional  field 

30-31 

Subcel 1 

6 

vegetation  code 

32 

Subcell 

6 

optional  field 

34-35 

Subcel 1 

7 

vegetation  code 

36 

Subcel 1 

7 

optional  field 

38-39 

Subcel 1 

8 

vegetation  code 

40 

Subcel 1 

8 

optional  field 

42-43 

Subcel 1 

9 

vegetation  code 

44 

Subcel 1 

9 

optional  field 

69-72 

Date  the  card  was  prepared 

77-80 

Identification  number  for  each  data  card 

I - 
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The  mapping  levels  are: 

Level  1 = Coastal  strand 

Level  2 = Coastal  sage  scrub-stabilized  dune  phase 
Level  3 = Coastal  bluff 
Level  4 = Oak  woodland 

Level  5 = Background  symbolism  (all  other  vegetation  types) 

Example  13: 

Assigning  Vegetation  Type  by  the  Majority  of  Subcells 

The  FORTRAN  statements  given  in  this  example  of  Subroutine  Flexin  direct 
the  computer  to  assign  an  overall  vegetation  type  to  a grid  cell  based  on  the 
vegetation  types  of  the  majority  of  subcells.  The  FORTRAN  statements  are: 


Explanatory  Notes : 

SUBROUTINE  FLEX  IN ( IFORM,  T,  FIRST) 

REAL  A(9) , B(9) 

\ DO  100  1=1,9 

‘ 100  B ( I ) =0 . 

READ ( 1 1,11)  (A(I), 1=1,9) 

11  FORMAT ( 1 8X  ,9 (2X  ,F2 . 0) ,26X ) 

DO  200  IN=1 ,9 
DO  300  1 = IN ,9 

300  IF(A(IN) .EQ.A(I))  B(IN)=B(IN)+1 
200  CONTINUE 
T=A(1) 

DO  400  1=1,8 

400  IF ( B( I ) . LT . B( I +1 ) ) T=A(I+1)  a. 

• IF(T . EQ. 3. ) GO  TO  590  b. 

IF (T. EQ.31 . ) GO  TO  590 
IF(T. EQ. 1 . ) T=5. 

IF (T . EQ. 5 . ) GO  TO  600 
IF (T. EQ. 12 . ) T=l. 

IF (T. EQ. 1 . ) GO  TO  600 
. IF (T . EQ . 2 . ) T=5. 

I F ( T . EQ . 5 . ) GO  TO  600 
I F ( T . EQ . 8 . ) T=2. 

IF (T. EQ. 2. ) GO  TO  600 
• I F (T . EQ .11.)  T=3. 

IF (T. EQ. 3. ) GO  TO  600 
T=5 . 

IF (T. EQ.5. ) GO  TO  600 
590  T=4. 

600  CONTINUE 
RETURN 
END 


Explanatory  Notes  for  Example  13: 

a.  The  value  of  T is  the  vegetation  type  of  the  majority  of  subcells. 

b.  The  remainder  of  the  FORTRAN  statements  direct  the  computer  to  produce  a 
vegetation  map  of  5 levels,  based  on  the  instructions  provided  to  the  Map 
Package.  The  mapping  levels  are: 
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Level  1 = Coastal  strand 

Level  2 = Coastal  sage  scrub-stabilized  dune  phase 
Level  3 = Coastal  bluff 
Level  4 = Oak  woodland 

Level  5 = Background  symbolism  (all  other  vegetation  types) 

Examples  of  different  uses  for  the  vegetation  data  base  include: 

a.  Mapping  and  calculating  areas  sensitive  to  erosion  due  to  unstable 
vegetation. 

b.  Predicting  occurrences  of  rare,  endangered,  or  depleted  species  of 
plants  and  animals  based  on  their  vegetation  preferences. 

c.  Mapping  and  calculating  areas  of  different  vegetation  types  that  would 
be  affected  by  alternative  placements  of  construction  (runways,  launch 
pads,  etc.). 

6.  MERGED  DATA  BASE 

In  order  to  do  analyses  with  GRID  using  more  than  one  variable,  a merged 
data  base  was  produced.  The  data  base  variables  for  each  grid  cell  were 
merged,  resulting  in  a common  data  set.  The  merged  data  base  is  similar  to  the 
vegetation  data  base  in  that  one  computer  data  card  is  used  to  describe  the 
data  variables  for  each  grid  cell.  Therefore,  the  merged  data  base  contains 
4646  data  cards. 

Table  22  lists  the  fields  and  columns  describing  a grid  cell  for  each  data 
card.  The  merged  data  base  is  available  for  GRID  and  other  computer  programs 
on  either  data  cards  or  magentic  tape  at  VAFB. 

7.  UPDATE  PROCEDURES 

The  merged  data  base,  because  of  its  design,  can  easily  be  updated  on  a 
cell  by  cell  basis.  Example  14  illustrates  the  procedure  for  updating  data 
cards. 

Example  14: 

Update  Procedures  for  Data  Cards 

New  firebreaks  were  made  in  a number  of  areas  on  the  base.  It  is 
desireable  to  update  the  data  base  to  include  the  new  firebreaks.  The  following 
procedure  can  be  used: 

a.  Locate  the  new  firebreaks  on  the  Base  Master  Plan  Maps  (C-1.2  series). 

b.  Identify  the  coordinates  of  the  grid  cells  containing  the  new  fire  breaks  -- 
for  example,  0A100,  PA100,  and  QA100. 

c.  Locate  the  data  cards  in  the  data  set  which  display  the  coordinates. 

d.  Keypunch  new  data  cards  to  replace  and  update  the  affected  data  cards. 

e.  Return  the  updated  data  cards  to  their  proper  location  in  the  data  deck. 

Data  values  can  be  updated  on  the  magnetic  tape  by  having  the  computer 
search  the  data  tape  for  the  desired  cell  coordinates  and  writing  the  new  data 
values  into  their  proper  location  on  the  tape.  Example  15  illustrates  the 
procedure  with  a FORTRAN  program  to  update  firebreaks  in  subcell  5 for  the  grid 
cells  0A100,  PA100,  and  QA100. 
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TABLE  22.  MERGED  DATA  DECK  FIELDS. 


Col  umri 

Field  Description 

2-3 

S= 

4-6 

Soil/slope  code 

8-10 

EX= 

11-12 

Exposure  code 

14-16 

EL= 

17 

Elevation 

code 

19-20 

V- 

21-22 

Subcell 

1 

vegetation  code 

23 

Subcel 1 

1 

optional  field 

25-26 

Subcel 1 

2 

vegetation  code 

27 

Subcell 

2 

optional  field 

29-30 

Subcel 1 

3 

vegetation  code 

31 

Subcel 1 

3 

optional  field 

33-34 

Subcel 1 

4 

vegetation  code 

35 

Subcel 1 

4 

optional  field 

37-38 

Subcel 1 

5 

vegetation  code 

39 

Subcell 

5 

optional  field 

41-42 

Subcel 1 

6 

vegetation  code 

43 

Subcel 1 

6 

optional  field 

45-46 

Subcell 

7 

vegetation  code 

47 

Subcel  1 

7 

optional  field 

49-50 

Subcell 

8 

vegetation  code 

51 

Subcell 

8 

optional  field 

53-54 

Subcel 1 

9 

vegetation  code 

55 

Subcel 1 

9 

optional  field 

66-69 

Date  the  card  was  prepared 

71-75 

Alphanumeric  coordinates  for  the  grid  cell 

77-80 

Identification  number  for  the  data  card 
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Example  15: 

A FORTRAN  Program  to  Update  the  Data  Tape  for  Firebreaks 
1 

2C0MPILE  SNFT10  XFORTN  LIB  DATA  CARDS 
3 

11 

3LST1 

8 

IDENT  MAIN 

FILE  1 1=ANFT0T , UN  I T=TAPE ,REC0RD=80 , BLOCKINGS .UNLABELED  ,F  IXED 
C UPDATE  PROGRAM 
INTEGER  V,  ID 
REAL*8  C(3) , CX 
V=1 

NTYPE  = 3 

DATA  C/'OAIOO',  * PAIOO* , ‘ QA1 00  * / 

1 READ (11  ,2)  CX,  ID 

2 FORMAT (70X.A5,  IX,  14) 

DO  3 J=1  , NTYPE 

IF (CX.EQ.C(J) ) GO  TO  4 

3 CONTINUE 

IF ( ID. LT. 4646)  GO  TO  1 

4 WRITE(4 ,5 ) V 

5 FORMAT (38X,  II) 

IF ( ID. LT.4646)  GO  TO  1 
END 


1 

2END 


J 

APPENDIX  A 

GRID  PROGRAM  LISTING  AS  ADAPTED 

I 

TO  THE  BURROUGHS  3500  COMPUTER 


¥ 


IDeNT  MAIN 

i-  ILL  11  = ANF  T'TtUM  T = TAPt,WtC /)kO  = U.QL  JCK1NG  = 1 « L ABt  LE U . E I Xr  O 

L ARGt  C.  A,  L#  L IM 

SEGMENT  MAPI N,  IN OAT  A , FL  A TON, HmLAD , FLEX l N 


OHIO  PR 

GRID  MAP  PROGRAM  VERSION  3 


PREPAkF.O  O Y DAVID  h»  SlNTON 

LABOf-ATOKV  FOR  COMPUTER  GRAPHICS 
marvapc  university 
AUGUST  1969 


O G R A M GROOwOi.) 

G4D0  •.  c ao 
G IDO(.  «J  JO 
G-luQO  04  0 
o ID-  iiv5u 
G IDO  0 0 60 
G 4D0  C 0 7 0 
G IDvOOeO 
G IOC  w'O  90 
w ICO 
JCGRDOOl  I 0 
G IDO  w 1 20 
I T GiTDO  ' l 30 
UPT I UN  .GROOO 140 
GROOw ISO 

MODIFIED  FUR  VANDENBFRG  AFT  STUDY  OCTOBER.  1974  F • STUT L • R . RE  I L L Y ol0Gul60 


MODIFIED  BY  T “'M  HAkTSuA 


v.E  1G.  DLPT*  S A*  j JltGG  STATE  UNIV. 


THIS  VERSION  OF  GRID  IS  MODIFIED  TO  RUN  UN  IBM  360/4'  (DOS) 
ACCEPT  10.00C  DATA  CELLS  WITHOUT  USING  THE  MULTIPLE  DATA  SET 


C 

CUM40N/C/  MlCCOO) 

COMMON/ A / N CD  * NCAf  NGO*  NGA. NSHEET  *NCS»NCF.NCL  *NUMCH, NC  ST • NRS  T 
CJMMlN/B/  NLtVtL.  I FORM  ,NUD.  IHI  ST.NTX#0PT(25> 

COMMON  /L  IM/  VAi.max  , V 4 l M 1 n , PANGb  (20  I , VAL(  2 1 ) ,PRANGt(  20  ) . NF  RE  J ( 25  ) 
COMMON  /SYMFLS/  SYMO  it  <25,4  t «NS(  1 0 ) • T l TLE  ( 60  , TEA  T ( £>Jd) 

INTtG^R  SYMBOL 
DIMENSION  IAC2..0) 

LOGICAL  OPT 
DIMENSION  *FyS<5) 
l INTEGER  T I * l , T IM2,  T IM  3 , T I M4 

C INTEGER  l 

r l-  11 

DATA  K f YS/ JHM AP, 3HENU » 3H 1R R . JHDA  T . 3H9  99/ 

OPT (24 )=*FALSE • 

JPT  < ?5  ) - • F At  Sfc  . 

mapn=*. 

ILKRs  ' 

r -cau  control  vord 

1m  I F ( OPT  C 25) ) IEPR=IERRF1 
RtAU(S.IODu)  KEY 
OPT  ( 2r>)  = .F  ALSE  • 

DO  l 1=1,3 

IF (KEY  .EO.KE YS(  I ) ) G.  TJ  12 

l continue 


ERROR  M £ S S A G ?-  S 

4R I TE( 6*  20  . u ) KEY 
*•?!  TF  (M,9Cw'C  ) 

UFlLt  4 
it  4 IND  f; 
i T UP 

12  *MI TEC  ft# 1GI"  ) 
l r c(d»9Ct  r 1 
r ndf ilf  h 

if  « I M)  tt 

5 tup 

11  vjLi  T ) (1,2,3,4,51,  I 
/ m APN-M APN-  I GftR 

4PITF.(6#40K’)  MAPN 


U 4D*.  0 170 
G IDO  s,  ISO 
oRDO  0 190 
G I DO  0 2 1 C 
G ■!  DO  0 2 l 0 

G ROD 0221 
6R&0 .230 
G I JO  0 24  0 
G IDO C 250 
GRDO .?30 
GRDO  0270 
o ID  ) . 25  C 
GRD00290 
G -IDs  ) 30  0 
G 4D0  ) 3 I 0 
GRDJ  320 
GRDOO 33C 
G -ID  ) v J4  0 
GRDOO ISO 
G ROC . JoO 
G IDO  ' 37  w 
GIDO'-  380 
G ROD  39  0 
G fiO  J w 4 U 0 
G4D0r  4 1 0 
GRD-»  420 
G IDO  0 4 JC 
G IDO  H40 
G IDO  / 4 5 C 
G IDGj46j 
G 4000470 
GRDo  '480 
G IDO  . 4 90 
G 400  O 500 
G RO  'J  ' 510 
G IDDL52Q 
<j  IDC.  C S 30 
G4D0  JS4C 
URDOOSSO 
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O A 


P 


/ 

i M t*  fc  2 1 

C*  'w'  «1 

.*KlTE(Z,90wC  > 
fc.MU  FILE  Z 
^r.fclND  Z 
STOP 

4 WRl TE( 6. 2C  1C  ) 

WRITE! 8.90CC  ) 

LNDr  ILL  8 
REWIND  8 
STOP 

5 (j'J  TQ1CO 

r 

C RFAD  IN  IRREGULAR  OUTLINE 

a I P=0 

JPT (24  ) = .TRUE • 

REWIND  9 
WRITE ( 6. 1 050 ) 

106  READ (5 § 10 3C  I lK.NC8.NCE 

IF  (IK  .NE.99999)  GG  TU  107 
fJPT  ( 24  ) = IP  . GT  .0 
GU  TU  100 

1 w 7 IF ( IK.LT  . 1 ) l K — 1 
DC  105  1=1. IK 

IPs IPf 1 

WR I TF( 6.  1 C3C  ) IP.NCB.NCE 
1 r>  W»ITE(9)  NCB.NCE 
GO  TO  1 C 6 

i" 

C READ  IN  THE  MAP  CONTROLS 

1 MAPN=MAPN+1 
C CALL  CTIME(TIMl) 

ND=  0 

WFI  TE  ( 6, 1 0<*C  1 
CALL  MAP  IN(MAPN) 

I F ( OPT (25))  GO  TO  IOC 


RE  AC  THE  DATA 
2 . : ND=  ND 1 

I F ( OPT ( 2R ) .AND.  IP.NE.NCD)  GU  TO  120 
CALL  1 ND A TA(NCD*NCA.ND.  I A) 

I F ( OPT ( 25 ) ) GO  TU  ICC 

C loop  through  flatqn  for  each  sheet  of  the  map 

c 

r call  c t ime ( t I M2  > 

I f.H*  2 

3 i SH=  1 SH'f  1 

WRITE  (8.3oC0> 
w - ITE( 8.40  ,c  » MAPN,  I SHf ND 
WRITE  (8.3043) 

CALL  PL ATON( ISH.NO) 

•RITE  (8, 3v4C ) 

IF  ( ISH.L T .NSHEET ) GO  TU  3y ) 

IKNO.LT.NUD)  GO  TO  20  0 

C CALL  CTIME(TIM3) 


I rr  the  data  beljw  the  map 

white  ( R. 3 6C  ) T I TLE 

IP  ( NT* • G T ) WRITE  (6* s 80)  (TEAT!  I ) 
IF(PPT(12)  .AND.  PT ( 2 2 ) ) GJ  TJ  210 
CALL  HHEAD 
CAi  l CT  IME( T IMO ) 


I = 1 .NT  A) 


ij  iJO  C 560 
GRD0C570 
GRD0058C 
GRDOCbRO 
GRDO 06 CO 
GRDO 06 1 0 
GRUO  Gc  20 
GRDO 06 30 
G *00^640 
GRD0C650 
G<  DO  0 660 
G *000670 
GRDO  j 680 
GRDO  069  0 
GRD0C7C0 
GRD007 10 
GRD00720 
GRDO  0 7 3 0 
GRDOo  740 
GRD00750 
GRDO  0 760 
GROOC  770 
GRD0078C 
GRDO  0 790 
GRDC- 800 
GRD008 1 0 
GRDO  0 820 
G *000830 
GRDO  0 84  0 
GRDO  0850 
GRD0086C 
GRD0087C 
GRDOCQ0O 
GRDO  0890 
G *000900 
GRDO  09  i 0 
GRDO  092C 
G *D0  09  30 
GRDO  094  0 
GRDO  09o0 
GRDO  0960 
GRD00970 
G *000980 
GRJOs.990 
GRDOICCO 
GROG  10  10 
GRDO 1020 
G RDO  1 O 3 0 
oRDC 1 0 4 C 
6R6010  0 
GRDO 1060 
GRDO 10  70 
GRDO 1 080 
GRDO 1090 
GRDO 1100 
GRDO 1110 
GRDO 1 120 
GRDO  1 1 3C 
GRDO 1140 
GRDO 1 1 50 
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2 1 Q IF  ( .NOT  .OPT  < 1 5 I)  ,C)  TO  i 00 
C TI*(T IM2-TIM1 )/300. 

C T2=( T IM3-T I M2 )/300 * 

C T3=(T IM4-T IM1 )/300 . 

C T4*  T I M4/ 300  • 

C JEl 6. 9C50 > TI#7?*T 3,^* 

Gj  tu  ioo 
1000  F URM  AT  (All 

20  Ovl  FORMAT  (2H  -*A3.34M-  IS  AM  INVAI. 

1C1J  FURMAT  i//2X,  * 2H  IKQFC.O t..ip  0»T1 

1 EH  OF  P OWS,  ) 

4C10  FORMAT  < 1H  l , I Sf  l 1 H MA«>  MA0c~/l‘5M  PUN 
<01«  Ff.'RMAT  DATA  IS  NO  LONGER  * V*i.lO 

10-jC  FORMAT  < 1HI  ,«(/ I,  ! x#  I 2,Hir<Ne<iO«.  A . j u’,  !\ 

I 0 30  FORMAT  (JIM 

ICAO  F JHMAT  ( lnt  t B(  / ) • 1 X • ».2-‘AP 

1C00  FORMAT  I 1M1  • 8 ( ✓ I I 

aCCO  FORMAT  ( I X * 4HM  AP  v I 2*fH>SHC£  - , ,»  CM, OAT  4 

JCaO  FJPMAT  ( 2H  »,  I2(  l.'H ~ » I O'  1 

3t  60  FORMAT  I ✓✓// (1 X,?iA4/ I X )✓/ | 

33  80  FORMAT  (/ ( 4.CX  t 20AA  ) | 

9000 

L 9C  s: 


i ••  > 

FORMA'  (//*  time  In  seconds 

^ AP  • ,FH| i/I 7A» •-  TOTAL  * U- 
cAO 


G3O0 1 loO 
GRDO 1170 
GROO 1180 
GROQ l l 9C 

G»r>0  1 200 

GROO 1210 
GROO 1 220 
GROJ I 230 
GRDO 1 240 
CORRECT  MUMBGRD01250 
GR. >01260 
GROO 1 ? 70 
GRDO I 200 

. ) o^00 1290 

GRL)0  1 300 
,/)  GRD01310 

GROO  1 3 '0 

1 G9D01J30 

GRl)0  I 340 
GRDO 1350 
GROO I 360 

>*<***44  £ NO  JF  DATA  ► 1 1 £ A 4 **  4 * * A m GROO  1 3 70 

GRDO I 380 

k MPUT  afFa«3/|7\  * - PR  | NT  INvi  -0900  J3«0 

3.3/1 fX,»-  TQT*L  l t APS6D* ,F8. 3)  GR001«00 

6 903  I 4 '' 


to  input  cxn.pgl  wjroi 
*•  h S DO  vor  iPt.  C 2 F Y »'  he 


C Mpi.  F TF  I 
C N f POL  I 
• / i « # 1 7 H~  - 


. / 1 H . I ? H 


56  T 

-4) 
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I 


IDtNT  MAPIN 

SUBROUTINE  MAPIN(MAPN) 

COMMON/ A/  NCD,NCAf  NGD  f NO  A , NSHE  Fi  T , NCS  , NCF  , NCL  , NUMCH • NC  ST • NRST 
COMMON/Ei/  NLEVEL , I FORM  ,NUO  * lHlST,NTX,CiPT(2S) 

COMMON /L IM/  VALMAX .VALM1N, RANGE C 20  >, VAL(2l  ) ,PRANGE(20  ) , NFREQC  2 5) 
COMMON  /SYM8LS/  SYMBOL ( 25, 4 > fNS( I J ) f T I TLE ( 60 ) , TEX T ( 800 ) 

LOGICAL  OPT 
INTEGER  SYMBOL 
DIMENSION  A ( 6 I 

DATA  TX TEND/4 Hf  NDT / , P ACEND/4H9999/ 

IF (MAPN.GT. I ) GO  TO  200 

NT  X = 0 

N0D=0 

NGD=  4 

NbA«b 

NLEVEL  = 10 

NUMCH= 1 

I H I ST  =1  _ 

00  201  1=1,23 

20 1 OPT  ( I ) = .FAL  SE • 

?CC  READ  ( 5 , 1C CO)  TITLE 
WF I TE ( 6 , 1 0 1 0 ) TITLE 
WRITEC  6, 1020  ) 

ICO  RE  AO  ( 5,  1030  1 I OP T , ( A ( I I , I = 1 f 6 1 
IF(  IOPT  ,EQ,99999)  GO  T3  300 

GO  TO  ( 1,2, 3,4,5, 6,7, 8,9,10, 11 ,12,  13,14,15  l.ICPT 
WHITE <6,9000 > It'PT 
GO  TO  100 

3 ■>  0 I F ( # NO  T , OP  T ( 1 2 ) , A NO  • NLEVEL  ,G  T , 1 3 ) NLE  VE  L = 1 0 
IF (OPT  <111  GO  TO  30  1 
W R I TE ( 6, 900 1 1 
i'jPT  ( 25  1 = .TRUE  • 

RETURN 

301  I F ( OPT (711  RETURN 
WRI TE(6,90C21 

OPT ( 25  1= .TRUE  . 

RETURN 

C 

OPTION  1 READ  T MI  GRID  CONTROLS 
I OPT ( 1 1 =, TRUE. 

I F (OPT (12  11  WRITE! 6, 90  0 3 1 
I E ( OPT  (13  1)  WR  I TE(  6, 9004  1 
?PT (121 =. FALSE, 

PT ( 1 3 1*, FALSE , 

NCD=A<  l 1 
NCA=A ( 2 1 

IF (A(  3 l.tO.C.O)  GO  TO  109 

NGD  = A ( 3 1 

NG A= A ( 4 1 

IF (NGA  ,EQ,  l .OR.  NGD.EO.ll  UP T ( o 1 = . T W Ut . 
liO  I Tt( 6, 2CvC  1 IOPT,NCD,NCA 
v*«I  TE(  o,  20  l C 1 NGD  • NG  A 

SET  PARAMETERS  AND  CALCULATE  THE  NUMBER  OF  SHtETS 
1 NCS=129/NGA 

NSHE t T= ( \C A/N CS ! ♦ 1 

1 h ( MOO(  NC  A ,NC  S I ,EQ  ,0  1 NSHtF  T =N  SHfc  L T-  1 
IF (NSHEET ,FQ.  1 1 GU  TO  ICO 

IF  MiJRF  THAN  ONF  PAGE*  - SET  NUMBER  OF  CELLS  PER  SHEET 
112  P = ( N C S*  N sHF  FT  1 - f . C A 

NCRH  = NCr'/2 


GRDO 1419 
GRDO 1420 
GRDO 14 JC 
GRDO 1440 
GRDO 1450 
GRDO 1 460 
GRDO 1470 
GRDO 1 *7 1 
GRDO 1480 
GRDO 1490 
GRDO 1 500 
GRDO 1510 
GRDO l 520 
GRDO 1530 
GRDO 1540 
GRDO 1 550 
GRDO 1 560 
GRDO 1 570 
GRDO 1 580 
GRDO 1 590 
GRDO 1600 
GRDO 1610 
GRDO 1 620 
GRDO 1 630 
GRDO 1 640 
GRDO 1 650 
G ROO 166  0 
GRDO 1670 
GRDO 1 680 
GRDO 1690 
GRDO 1 700 
GRDO 1710 
GRDO 1 720 
GROO 1 730 
GROO 1 740 
GRDO l 750 
GROO 1 760 
GRDO 1 77C 
GRDO 1 780 
GRDO 1 790 
GROO 1 800 
GRDO l 8 10 
GRDO 1620 
GRDO 1830 
GRDO 1840 
GRDO 1 850 
GRDO 1 860 
GRDO 1870 
G RDO 1 600 
oRD  D 1 890 
GRDO 1900 
GRD01 91 0 
GRDO 1 920 
GRDO 1 930 
GRDO 1 940 
GRDO 1950 
GRDO I960 
GRDO 197C 
GRDO 1 960 
GRDO 1 99C 
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p **  r>  \J 

i I 


NCF=NCS-NCRH 
NCL=NCS-NCrtfNCRH 
GO  TO  100 
C 

2 0PT<2)  = A<  i ) • GT  *0 • 0 
OPT ( 20 ) = A( 1 l.LT.G.C 
1 F QRM=  A ( 1 ) 

NOD=A(  3 ) 

W R I T E ( 6 • 30L  0 ) IOPT.IFORM 

IF(NOD.GT.l)  WR ITEC6, 3020  > NOO 
GO  TO  100 

C 

OPTION  3 READ  THE  NUMBER  OF  LEVELS 

3 Nt_E  VEL  = A ( 1 ) 

WRITE(6*1C4G>  I3PT.NLEVEL 
GO  TO  100 

C jPT I CN  4 MINIMUM  VALUE 

4 VALMIN= A(  1 ) 

JPT  ( 4 ) = A ( 2) .EQ.O.  3 

I F ( OP  T ( 4 ) ) WR 1 TEC  6*  1 050  ) IOPT.VALMIN 
GO  TO  100 

c 

C OPTION  5 MAXIMUM  VALUE 

5 V AL  M AX  = A ( 1 ) 

OPT ( 5 ) - A ( 2) .EG.O.C 

I F < OP  T C S I I »ft  I TE(  o.  1060  I ZJPT.VALMAX 

GU  TO  103 

f OPTION  6 VALUE  SCALING 

6 OPT (6  >=A<  1 ) .GT  .0.0 

IF ( .NOT .OPT ( 6 ) ) GO  TO  ICO 

00  1 G 2 l=l«o 

1 0 2 RANGE! I ) = A< I > 

IF (NLEVEL .GT .6  I RFA0C5,1G3C)  I B * ( RANGE  ( I ) * I *7 . NLE VEL ) 
WRI  TE< 6* 1080 ) ICPT * I RANGE (I  ) . I = 1 .NLEVEL  > 

GO  TC)  100 

C 

OPTION  7 SYMBOLISM 

7 OPT  C 7 ) = . TRUE  . 

RE AO(S, 1C9L > SYMBOL 
WRITEC6*  1 ICO)  I OPT .SYMBOL 

GO  TO  100 

( 

(.  3PTICN  8 FLAG  POINT  SWITCH 

8 OPT (0)  = A(  1 ) .EG.  1 .0 

IF ( OPT ( 8 ) ) WR I TE(6, 1 180 ) I0PT 
C OPTION  9 HISTOGRAM  SWITCH 

v*0  TO  100 

» JP  T ( 2 1 )=A(  1 ) .GT  .0 
1PT (22 )*A( 2 ) .GT .0 

1 F ( OPT (21))  «RITE(6. 1230  ) 

I F ( CPT ( 22 ) ) WRI TE< fc . 1240 ) 

GO  TO  ICC 

C GPTION  10  TEXT 

12-  J FORMAT  (/ I 5.  3X  . I I HMAP  TEXT  IS./8X.1IH ) 

I ■-  WMT£<6#  12C0I  IOPT 
NTX*C 


GRD0200C 
G TOO  20  1 0 
GRD02C20 
GRDO  20  30 
GRD02040 
GROO 2050 
GR002060 
G ROO  20  7 0 
GROO  20  80 
GRD02090 
GRDO  210  0 
GROO  2110 
GRJ02120 
GROO  2 1 30 
GROO 2 140 
GRD02150 
GROO  2 1 60 
GROO  2 I 70 
GROO 2 180 
GRDO  2190 
GROO 2 200 
GRDO  22 1 0 
GROO  2220 
GR002230 
GROO  224  0 
GROO  2250 
GRD02260 
GRD02270 
G RDO  2260 
GRDO  2290 
GROO  2 30  0 
GROO  2 310 
G ROO  2 320 
GROO  2 3 30 
G ROO  2 34  0 
GR002350 
G 3 DO  2360 
GR002370 
G_R002380 
GROO  2 390 
G ROO  240  0 
GRDO 24 10 
GRDO  2 4 20 
ORD0  2430 
G RDO  24  4 0 
GRDO  2450 
GRO02460 
GRDO  24  7 0 
GRDO  24  80 
GROO  2490 
GROO 2500 
GR0025 10 
GROO  2520 
GRDO  2530 
GRDO  2540 
GRDO  2550 
GRD02  5b0 
GROO  2 5 70 
GRDO  2580 
GRDO  2590 
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I 


t 


>0  K 6 N=1.40 

GRUG26G0 

I BEG  = N T X ♦ 1 

G RDO 26  1 C 

I ENO=NT  x*2C 

GRDO  26  20 

READ  (5#lOrC)  ( TEXT( 1 ) , I*I0EG, IEND ) 

GROO  2630 

IF ( TEX T ( I BEG ) • tG • T XTfcND ) GO  TO  107 

GRD0264C 

NTX= iend 

GROO  2650 

1 6 

CCNT  INUt 

GRDC  2660 

■ HTH6,1  220  1 

GROO 26 70 

i oa 

READ  ( b. 1000)  TXTNO 

GROO  2600 

IF ( TXTND.EQ.TXTENG)  GO  TO  107 

GRD02690 

IF( TXTND.EQ.P ACEND ) GO  TO  300 

GRD02700 

GO  TO  100 

GRDC27 1 0 

l 0 7 

• RITE  (6,1210)  (Tc_XT(I  ),I=1,NTX) 

G RDO  2 72  0 

GJ  TO  ICO 

G RDO  2 7 30 

c 

GROO  2 74  0 

1 1 

00  103  1=1,3 

GROO  2 7 50 

10  3 

OPT ( 84 I ) = A < I ) • EG.  1 • C 

GROO  2 760 

I F ( OPT ( 9 ) ) *RITE( 6,  1 1 40  ) 

GRD02770 

I F ( OPT (10))  WR I TE(  6, 1 1 50  ) 

G ROO  2 78  0 

I F ( OP  T ( l l ) ) MR  I TE ( 6 « 1 1 60  ) 

GRD02790 

GO  TO  100 

G RDO  200  0 

C. 

GROO  26 1 0 

c j; 

T PAP  OPTION 

GRD02820 

1 2 

OPT ( 12  ) = A(  1 ) , E 0 • 1 ,0 

GROO  2030 

IF( .NOT.OPT ( 1 2 ) ) GO  TO  100 

G RDO 2840 

*R1 TEC  6, 1 1 10 ) 10PT 

GROO205O 

NLt V EL  = 20 

GRD02860 

I F ( OP  T ( 5 ) ) NL  E VEL  =19 

GROO  28  7 0 

I F ( NGD  • EQ  • 4 .AND.  NGA.F0.5)  GO  TCI  100 

GROO  2080 

NGD  = 4 

GRD02890 

NGA  = 5 

GRD02RU0 

GO  TO  150 

G ROO  2 91 0 

G ROO  29  20 

C GRID  NUMBERING  OPTION 

GROO  29  30 

1 3 

NUMCH= A ( 1 ) 

G RDO  2940 

OPT ( 13)  = A(  1 ) • GT  .0.0 

GROO  2950 

I F C «NO  T»CPT  ( 13)  ) GO  TO  ICO 

G*^DO 

IF (NUMCH.GT .1  ) REAOI 5#  40  02 ) ( NSC  I >,  1 = 1 ,NUMCH) 

GRD02970 

NCST= A( 2) 

6-(Ow2980 

N»RST  = A( 3) 

GRD02990 

IF (NCST.EO.O)  NCST=1 

GRDO  3000 

IF (NHST .EO.O ) NRST=NCD 

6^0030 1 0 

NR  I TE  ( t ,40o0  ) I PT,NCST,nPST 

GROO  30  2 0 

IF ( ( NCA4NGA 420 ) ,GT . ( NSHEET* ( 1 3 -NGA ) ) ) GO  TU  111 

GROO  30  30 

GO  TO  100 

G ROO  304  0 

1 1 1 

NSHEET=NSH6ET  4 1 

GROO  3050 

GO  TO  112 

GROO 3C 60 

C 

GRDO  30  70 

C Ph 

-tSCALl-D  DATA  CPTI  N 

G ROO  30  80 

1 4 

l PT ( 14 ) = A(  1 ) . FQ  . 1 .0 

GRDO  3C  90 

•JPT(22)  = .TPUE, 

GROO  310  0 

IFCCPTC14))  *Pl TT (6, 1 190 ) 

GROO  3 1 1 0 

GO  TO  ICO 

GROO  3120 

1 J 

OPT  ( i b )=  A(  1 ) , G T , 0 

GROO 31 30 

IFC.PTC15))  *WITE(6,l25o) 

GRD031 40 

GO  TU  L) 

G RDO  3150 

1 w 

FORMAT  ( 2 C A 4 ) 

GROO  3160 

l ie 

r M A T ( IX,  20  A 4 , / ) 

GROO  31  70 

1 r ?. 

FORMAT C//1X#27HELE CT IVES  USED  FOR  THIS  MAP,/IX,27H 

1 

— — 1/, 

GRO0319G 

1 
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FORMAT ( I5.5X.6F1U.0)  GRD0320C 

FORMAT  (/I  5. 3X,  12. 7H  LEVELS)  GRD032I0 

FORMAT  (/I5,  3Xf  33HTHE  MINIMUM  VALUE  IS  SPECIFIED  AS, E 10.2)  GRD0J220 

FORMAT! /I5. 3X, 33HTHE  MAXIMUM  VALUE  IS  SPECIFIED  AS.FiO.2)  GRD03230 

FORMAT I/I5.3X.36MT HE  RELATIVE  SIZE  OF  EACH  LEVEL  IS  -./I0FI0.2/)  GR003240 

FORMAT ( 25 A 1 ) GRD03250 

FORMAT </I 5. 3X.  1 7HTHE  SYMBOLS  ARE  - . 4 ( / I OX . 25 A 1 ) / ) GRD03260 

FORMAT (✓ l 5. 3X . 2CHDOT  MAP  SYMBOLS  USED)  GRD03270 

FORMAT (/ JX, 33h I 1 THE  VALUES  MAPPED  ARE  ulSTED)  GRD032B0 

FORMAT (/3X. A2M1 i THE  VALUES  MAPPfcO  ARE  STORED  ON  CARDS)  GRD03290 

FCRMAT  C/3X.42HI 1 THE  LEVELS  MAPPED  ARE  STORED  ON  CARDS)  GRD033C0 

FORMAT  !✓ 15, 3X , 27MNO  FLAG  POINTS  IN  GRID  CELL)  GRD033I0 

FCRMAT </3X. 40H14  THE  DATA  IS  ASSUMED  TO  BE  PRESCALED)  GRD03320 

FORMAT ( 2uX, 20 A4 ) GRO03330 

FORMAT <///20X , 22HTOO  MANY  LINES  OF  TEXT)  G3003340 

FORMAT !/4X ,' 9 HISTOGRAM  INCLUDED  AFTER  FREQUENCIES')  GRD03350 

FORMAT !/4X, • 9 NUMERIC  INFORMATION  SUPPRESSED')  GRD03360 

FORMAT  !/I5,  • TIME  FOR  EACH  MAP  PRINTED')  GRD03370 

FORMAT  (/I  5, 3X , 12HGR  ID  SIZE  ISVI4,I5H  CELLS  DUNN  AND, I4,13H  CELLS  AGRD03380 
LCROSS)  GRD03390 

FORMAT!  BX.J2HCELL  SIZE  IS.14.20H  CHARACTERS  DDKs  AND.I4.1BH  CHGRD03400 

) AR ACTE  RS  ACROSS)  GRD034  10 

FORMAT  !/ I 5, 3X , 20HTHE  FLEXIN  LOCATOR  =.13)  GRD03420 

FORMAT ( 8X , 1 9HTHE  OATA  IS  READ  IN.I3.5H  SETS)  GR003430 

FORM AT  ( / I 5. 3X . 24HGR ID  NUMBERING  BEGINS  AT.2I6)  GRD03440 

FURMAT ! 10  I 1 ) GRD03450 

FORMAT  (/IS.  3X.  '•»•****  INVALID  ELECTIVE  NUMBER  ♦Hmm')  GRD03460 

FORMAT !//37H  ELECTIVE  l NOT  SPECIFIED  - MAP  STOPS)  GRD03470 

FORMAT (✓/ 37H  ELECTIVE  7 NOT  SPECIFIED  - MAP  STOPS)  GRD03480 

FORMAT (//51H  ELECTIVE  12  IS  CANCELED  WHEN  ELECTIVE  1 IS  CHANGED)  GPD03490 

FORMAT !//5lH  ELECTIVE  13  IS  CANCELEO  »HtN  ELECTIVE  1 IS  CHANGLO)  GRD03500 

END  GRD035IO 
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- map 


IDF  NT  INDATA 

S U B MC  U T I N E INDATA  ( NCU .NCA ,ND,  1 A) 

COMMON/C/  PC1GC00) 

COMMON/B/  NLEVEL,  IFORM,NOO, I HI  S T , NT X f OPT C 2 b ) 

COMMON /L  IM/  VALMAX . VAL M 1 N . R ANGE ( 20 ) • VAL (2 1 ) .PR  ANGE ( 20  ) .NFRE  GC 2b) 
CUMMUN/STTS/  STATC7) 

LOGICAL  FIPSTfOPT 
DIMENSION  l AC  I FORM  I 
v,  INTEGER  2 

C 2=11 

I F ( OP T C 2 4 ) ) REWIND  9 
noat  A=  NCAANCD 
F InST*  .TRUE • 

L READ  IN  THE  OATA 

C NLTE  OPT  C 24  )=.TRUE • FOR  IRREGULAR  OUTLINES  AND  OPT < 2 ) = • TRUF • FUR 
L READING  A LOGICAL  FILE 
L=  1 
M = NCA 

DO  ICO  I=l,NCO 
IC  = NCA*(  I- 1 ) 

IF ( .NOT .OPT C 24  ) ) GO  TO  111 
R E AD ( 9 1 MC  B « MCE 

IF ( ( MCE+MCB ) .GE *NCA > GO  TO  170 

M = NC  A-MCF 

L-MCfl 

IFIL.EQ.O  ) GO  TO  101 
DO  102  J=  1 • L 
I D=  I C-F  J 

102  P( I D >=-999999, 0 
101  L=L F 1 

111  IF C .NOT .OPT C 2 ) ) GO  TO  150 
DO  103  J=L.M 
I D= I C F J 

CALL  FLEX  INC  I FORM, PC  ID)  »F  IRST , ND  ) 

10J  F IRST= .FALSE. 

GO  TO  151 
1 c 2 CONTINUE 
153  CONTINUE 

L lb o REALM  1 4,END=152>  ( P ( I C ♦ J ) , J=L , M ) 

151  IF(M.EQ.NCA)  GO  TO  100 

M=  M ♦ 1 

DO  105  J=M,NCA 
I D = I C ♦ J 

1 w5  PC  I D >=-999999.0 
1C~  CONTINUE 

IF  C OPT  C 2 > .ANO.  ND.GT.NOD)  REWIND  12 

C 

C WRITE  UR  PUNCH  THE  UN  SC ALE  Q DATA  VALUES 

I F C OPT  C 9 > ) wR  ITF.Co,  l CIO  ) 

I F C JPT  C 9 > ) wRI TEC6# 1000 > ( N • PC N ) , N = 1 , NDA T A > 

IF  C OPT  C 10  ) ) WRI TEC  7,  1C0  0 ) CN»P(N)tN*l*NOATA) 

IFC .NOT. OPT  4 14 ) )GO  TO  14 1 

r SET  SC AL I D DATA 

DO  140  N= 1 | NOAT  A 

IF  ( PC N >• F O .-999999.0 ) PCN)=24.0 
PC  N > =P  ( N > FI 

GO  TO  2e ) 

l‘  INC  THfc  MAXIMUM  OR  MINIMUM  CF  DATA 


GRDO  3519 
GRD03520 
G 300  35  30 
GRDO  354 0 
GRDO  3550 
GRDO  3560 
GRDO  357  0 
GRDO  3580 
G R DO  3590 
GROO 3600 
GRDO 36 10 
GRDO 3620 
GRDC 3630 
GRDO  364  0 
GRDO  3650 
GRD03660 
GRD03670 
GRDO  3680 
GRDO  3690 
GRDO  3700 
GR0037 1 0 
GRDO  3720 
GRDO  3730 
GRDO  3740 
GRDO  3750 
GRDO  37b0 
GRDO  37  70 
GRDO  3 780 
GRDO  3790 
GRDO  3800 
GRD03810 
GR003820 
GRDO  3830 
uRDO  384  0 
GRDO 3850 
GRDO  3860 
GRDO  38  7 0 
GRDO  38  7 1 
GRDO  3800 
GRD03890 
GRDO  390  0 
GRD03910 
GRDO  3920 
GRDO  3930 
GRDC  394  0 
G ROO  3950 
GRDO  3960 
GRDO 3970 
GRDC  3980 
GRDO  3990 
GRU0400Q 
GRDO  4 0 1 0 
GRDO  4020 
GRDO 40  30 
GRDO  4040 
GRD04  050 
GRDO  4 060 
GRD04070 
GRDO  4 0 80 
GRDO  4 090 
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141  CONTINUE 

GR004 100 

IFCNO.GT.il  GO  TO  90 1 

GRD04 1 1 0 

PN=  0 

GR004 120 

PH  I a —9  99  999 • 0 

GR004130 

PL  0=  999999.0 

GR004 140 

PX1  =0 

G R DO  4 1 SO 

PX2  = C 

GROG  4 160 

PX  3 = 0 

GRD04 1 70 

PX4  = 0 

GRD04 1 80 

901  CONTINUE 

GR004  1 90 

00  204  N= 1 * NO A T A 

GRD04200 

IF(P(N  ) .EQ. -999999. 0 1 GO  TO  204 

GRD04210 

PN  = PN  ♦ 1 

GRD04220 

20  5 PL  0=  AM  I N 1 ( PL  O fP  ( N ) 1 

GRD042 JC 

PHl=AMAXl  (PHI  • P ( N 1 ) 

GROO  4 2 4 0 

PX 1=PX 1 ♦PIN 1 

GR0042S0 

P X 2=  FX  2 + P ( N 1**2 

GROO  4 26  0 

PX3=PX3*P(N 1**3 

GROO 4 2 70 

PX4  = PX4  FP (N )**4 

GRD04280 

204  CONTINUE 

GROO  4 290 

IF(NC.GT.l)  GO  TO  902 

GRQ0430Q 

IF( .NOT .OPT (4 ) ) VALMIN=PLO 

G R 00  4 3 1 0 

IF ( .NUT .OPT ( 5 1 ) VALMAXaPHl 

GRO04  320 
G R PO  4 330 

SET  PARAMETERS  TO  SCALE  THE  DATA 

GR004  340 

20  l VAL  ( 1 1 = VALM  IN 

GRD0435C 

A RANGE  = V AL  MAX  — V AL  M IN 

GROO  4 360 

IF ( ARANGE .LT .0 .000000 1 1 GO  TO  171 

GRD04 370 

IFf  .NOT .OPT (6  1 1 GO  TO  118 

GROO  4 380 

TWANGE=0 .0 

GROO  4 390 

DO  121  1 = 1*  NL  E VEL 

GROO  44  0 0 

121  T P ANlE = TR ANGfc ♦RANGE ( l 1 

GRD044 1 0 

118  VAL INC = AWANGE/NLEVEL 

GRD0442O 

DO  122  1 = 1*  NLEVEL 

GR0044 30 

IF (OPT (6)1  VAL  I NC  = R ANGE ( I 1 * ARANGE/ TRANGE 

GRD04440 

PR  ANGE ( I 1 = V AL  I NC  * 1 00.0/ARANGE 

GRQQ44  5 0 

122  V AL  ( !♦  1 ) =VAL  ( I )*VAL  INC 

GRD04460 

VAL  (NLEVEL ♦ 1 1 = VAL MAX 

GR0044  70 

902  CONTINUE 

GROO  4 4 80 
GROO  4490 

CALCULATE  THE  SCALED  DATA  VALUE 

GR004500 

199  DC  123  N= 1 , ND  A T A 

GROO  4b 1 0 

IF (P(N ) .FQ. -999999. C ) GO  TO  124 

GRD04520 

IF (P(N I .LT . VALM  IN > GO  TO  127 

GR004530 

DO  12b  1=1. NLEVEL 

GRU04S40 

IF(P(N) .Lc.VAL ( !♦ 1 ) ) GO  TO  126 

GRD04550 

125  CONTINUE 

GR004S60 

P(N  1 = 23*0 

GHD04S70 

GU  TO  123 

GR004S80 

124  P(N)=25.0 

G 4D04S90 

GO  TO  123 

GROO  4 60  0 

126  P ( N ) = I 

GROO  4 6 1 0 

GJ  TO  123 

GRD04620 

12/  P( N 1 =2 1 *0 

G ROQ  4630 

IF (OPT (121)  PfN)*G»0 

GR004640 

123  CONTINUE 

GROO  4 650 

IF(CPT( 1 1 >) WRITE! 7,10001  <N.P(N),N=l.NOATA) 

GHOO  4660 
GR004670 

CALCULATE  frequencies 

GROO  4 680 

IF (NO. GT. 1 1 GO  TO 90 3 

GR004690 

76 


20C  DO  129  1*1.25 

IF ( I.LT.0)  ST  AT ( I >*0 
1 29  NFREOC I )»0 

903  CONTINUE 

I F | OPT (14))  GO  TO  9C4 

STAT I 1 )*PN 

STAT(2)=PMI 

STAT ( 3 )=PLO 

STAT  C 4 )=PXl/PN 

STAT15  1*(  SORT  i PN*PX2-PX1  »'«2  > )/PN 
STAT 16 ) =0 
ST  AT  ( 7 ^0 

904  CONTINUE 

DO  203  N* 1 iNDATA 
I * P C N > 

203  NFREOC  I 1* NFREOC  I ) ♦ 1 
RETURN 

171  WR I TE ( 6 9 1 030  ) 

OPT ( 25  > * . TRUE • 

RETURN 

170  MR I TEC  6,  1020 > 1 

• RI TE<0,9OOO  ) 

REMINO  8 
STOP 

1000  FORMAT ( I5.5XFI0.2) 

101C  FORMAT ( 1M 1 , 8< /)•  1 X. 20HUNSCALED  DATA  VALUE S •/ 1 X . 20H 

1 ) 

1020  FORMAT C2X* 39HERR0R  - TOO  MANY  CELLS  REMOVED  FRCM  ROM. 15) 
1033  FORM  AT ( • DATA  RANGE  IS  ZERO  - MAP  STOPS*) 

90  C 0 FORMAT <✓/✓///////, * *, 46X ,•**** ******  END  CF  MAP  PRINT  QyT 

1 *•*  * ) 

END 


GR004700 
GRD0471 0 
GRD04720 
GRD04730 
GRD04740 
GRD04750 
GRD04  760 
GRO04770 
GRD04780 
GR004790 
GRDO40OO 
GRD04810 
G WOO  4 82  0 
GW  DC  48  30 
GRD04840 
G ROO  4850 
GRD04860 
GRD04870 
GWDO  4880 
GW004890 
GRD04900 
GRD0491 0 
GRD04920 
G ROO  49  30 
GRD04940 
GRD04950 
G WOO  4 96  0 
G T0049  70 
GWD04980 
'GRO04990 
GRD05000 
GRD05010 
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I DENT  FLATON 

SUBROUTINE  Fl  ATON(  I SH,  NO  ) 

CGMMCN/C/  PI  ICOOO  ) 

COMMON/ A/  NCD, NCA, NGO, NGA, NSHEET ,NCS ,NCF  , NCL , NUMCH , NC ST  . NR  ST 
CJMMON/B/  NLEVfcL . I FORM  * NOD . IHIST.NTX.OPTI25) 

COMMON  /SVMBLS/  SYMBOL ( 2 5 . 4 ) . N SI  I 0 ) • T I TLE I 60 ) . TEX T< 80 0 ) 

CCMMGN  /OUTPUT/  MAPALL ( 1 29 . 4 I 
O I ME  NS  I ON  XI  10 . 129 ) • INUMI  129.4 ) 

INTEGER  BLANK  * SYMBOL. SVM 
LOGICAL  OPT  t F I N 

DIMENSION  NUMI 10) . STCEL 14.5). IR14) 

C INTEGER  2 

C Z-  1 1 

DATA  BLANK/1H  / . PL  US  / 1 H ♦/ , S YM  l / 1 H l / , S Y M/ 1 H C / 

DATA  NUM/1 HO . 1 HI » 1 H2 . 1 H3 , 1 H4 . 1 H5 . I H6 . 1 H 7 • l H8 • 1 H9 / 

DATA  STCEL/ 16.0. 6. 0.1b. 0.0. 0.18.0,9. 0,14.0,1  1. 0, 4. 0.17.0,1  .0,12.0 
1 2.0, 5.0,10.0, 19.0, 20 .0,7.0, 13. 0,3.0/ 


SETS  THE  PARAMETERS  FOR  EACH  SHEE T OF  THE  MAP 
F IN  = .F  ALSE. 

N=NC  A 
MNCP=0 

IFINSHEET .EO. 1 ) GO  TO  500 
N=NCF 

IFIISH.EQ.l)  GO  TO  5C0 
N-NCS 

MNCP=NCF4NCS4 I I SH~ 2 ) 

IF(  ISH.EO. NSHEET)  N=NCL 
500  JF  lRST  = 64-( NGA4N )/2 
I X= INGO/2 )♦ 1 
JZ=NGA/2* 1 
KK  = 1 
NCH=  J 

IF 1NUMCH.GT . 1 ) NCH=4 
XRST=NPST-NCO*IND- 1 ) 
eCRDER=SYM I 

IF ( .NOT. OPT ( I 3 ) ) GO  TO  l 00 
K Z = MNCP- 1 
OC  200  Kxl.N 
KZ=KZ* 1 

NUB*  NC  ST  MKZ/NUMCHl 
INUMIK , 1 )=NUB/ 100 

I NUM (K  , 2 ) = NUb/ 1 0-  I NUMI K , 1 | ♦ 10 
INUMIK , 3)  = NUB- J NUM IK.)  )♦  100— INUM( K ,2)  410 
IN  = MCD IK  Z.NUMCH) 

INUMIK , 4 )*NS1  IN*  I ) 

• Rl TEI 8,  . JCO ) I BLANK, BORDER  , I SYMBOL  <25. 1 ) • J*  1 • 129 ) .BORDER .1  = 1.3) 

IF  I .NCT .OPT  I 1 3 ) ) GO  TO  180 
DO  201  J* 1 , 1 29 
DO  201  K* 1.4 
MAPALL I J.K  > = SYMBUL  125.1) 

0RITEI 8, 1002) 

100  2 FORMAT  1 32 X , • BCDE F GH I JKL MNOPQRS T U ViX Y Z ABCD6FGH I JKL MNUPQRSTU VMX Y 2 At) 
lOtFGHI JKLHNOPQ* ) 

4F  I TF I 8 , 1003) 

1-0J  FORMAT  ( J 2 X • ' AAAAAAAAAAAAAAAAAAAAAAAA) 30BBB RBBHBBBBQUBHBBBBBBBdBC C 
2CCCCCCCCCCCCCC* I 

180  MRITEIB.1000) I BLANK .BORDER ,<  SYMBOL! 25 , l ) , J = 1 , 1 29 1 .BORDER ,1*1,3) 

IF (FIN)  WFTURN 
FIN* .TRUE. 


203 

100 


2)1 


GR00501 9 
GROO  50  20 
GROO  50  30 
GRD05040 
GRD05050 
GRD35060 
GRD05070 
GRDO  50  00 
GRDO 5090 
GRO05 100 
GRD051 10 
GRO051 20 
GRDO  5 1 30 
GRDO  514  0 
GRDO  5 1 50 
, GRDO  5160 
GRD05 170 
GRDObl BC 
GRD051 90 
GRD0S200 
GROOS21 0 
GRDO  5220 
GRD05230 
GRDO  52 4 0 
GRD05250 
GRDO  5260 
GRDO  52  7 0 
GRDO  5280 
GRDO  5290 
GRU05300 
GRD05J1 0 
GRDO  5320 
GRDO  5330 
GRDO  5 340 
GRDO  5350 
GRU05360 
GR005370 
GRDO  5380 
GRDO  5390 
GRDO  540  0 
GRD0S4  10 
GRDO  54  20 
GRDO  54 30 
GROO  54  4 0 
GRDO  545  0 
G R DO  54  6 0 
GRDO  54  7 0 
GRDO 5480 
GR00549C 
GROO  550  0 
GRDO  5510 
GRDO  5520 
CGR  005530 
GRDO  5540 
GRDO  5550 
CGRD05560 
GRDO  55  70 
GRDO  5580 
GR005590 
GRD05600 
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GRID  CELL 


181  DO  124  IRUM=1,NCD 

BLANKS  GUT  A ROW  TO  BACK  GROUNO  SYMBOLISM 

00  101  1*1, NGD 

DO  101  J*  1*129 

IFC  I .GT.4)  GO  TO  101 

MAP ALL  C J,  I > = SYMBOL!  25,1) 

101  X ( I , J) =25.0 

INITIALISE  THE  LOOP  FOR  EACH  ROM  OF  CELLS 
JB=  JF I RST 

ID* ( IPG*- 1 )*NC A4MNCP 

DO  130  K=1.N 

ID  = ID  ♦ 1 4 

IF (P(  I D I .fcQ  .25  .0  I GO  TO  130 

INSERTS  SYMBOLISM  KEY  FOR  A GRID  CELL 

DO  132  I=1.NGD 

00  132  JY=1,NGA 
J=JB4JY 

IFI0PTI12))  GO  TO  131 
X ( I , J > = P(  ID  ) 

IF  ( OPT  I 8 ) ) GO  TO  132 

1 F < I *N£  • I X .OR.  JY.N6.J2*)  GO  TO  132 
XI X, J)*X( I , J ) + 1 0 

IFCPC  10) #GT . 20  *0 ) X(  I . J)  = X(  I • J )-9 

GO  TO  132 

131  M AP  ALL  C J , I ) = BLANK 

I F | ST  CEL ( I. JY  ) .LE.PI  10))  MAPALLI J.  I )=SYM 
13?  CONTINUE 
130  J0=  J 84  NGA 


00  124  M= 1 , NGD 

IF  COPT  C 12 ) ) GO  TO  191 
CONVERTS  A LINE  TO  SYMBOLISM 
00  103  J=l, 129 
JSYM*X(M. J) 

00  103  K = 1 . 4 

103  M AP  ALL  ( J , K ) = SYMBOL ( JSYM,KI 
GO  TO 190 
191  KK=M 

190  IFC .NOT .OPT f 1 3 ) ) GO  TO  220 
IFCM.NE.IX)  GO  TO  220 

IFCISH.GT.l  .AND.  1 SH • L T . N SHEE T ) GO  TO  220 

NUH=XRST- C I ROM- 1 l/NUMCH 

I«<  1 )=NU8/1 00 

IRC 2)= NUB/ 1C- IRC  1)410 

IR  C 3 ) = NUB  — IRC  1 )* 1 00— IR  C 2 )4 1 0 

IN=NUMCH-M0  0C  C I ROM- | ) , NUMCh ) 

IRC  4 ) = NSC  IN ) 

JX= JF I qst-h 

00  223  J=  l # 2 

DO  224  1=1. NCH 

1 I* IRC  I )♦  1 

MAP ALL ( JXfKK )=NUMC  I 1 ) 

224  JXsJX+1 
223  JX=JF IRST ♦N4NGA46 
•RITES  A LINE 

220  MR  I TEC  8.  I0C0  ) BLANK. BORDER  # C MAPALLC J,KK)*J*| .1  29  ) . BCROER 
IF  COPT  C 12 ) ) GO  TO  124 
•RITES  OVERPRINT  LINES 


GR005610 
GROO  5620 
GROO  5630 
GROO 5640 
GROO  565  0 
G ROO  566  0 
GR005670 
GROO 5680 
GROO  5690 
GR005700 
GR0057 10 
GRD05720 
GRD05730 
GROO  57  4 0 
GROO  57  50 
GROO  5760 
GROO 57 70 
GRO05780 
G ROO  5790 
GR005800 
GROO  58 1 0 
GROO 5820 
GROO  5830 
GROO  584  0 
GR005850 
GR005660 
GROO 5870 
GROO  56  80 
GROO  5890 
GRD05900 
GROO  59 1 0 
GRD05920 
GRD0S930 
GROO  594  0 
GRD05950 
GROO  5960 
GR005970 
GRD05980 
GROO  5990 
GROO 6000 
GROO  60 1 0 
GRD0602  0 
G ROO 60  30 
G ROO  60  4 0 
GR006v50 
G ROO  60  6 0 
GRD06070 
GROO 60 80 
G ROO  60  90 
GROO 8 1 00 
GROO 6 I 10 
GROO  6 1 20 
G ROO  6 130 
GR006140 
GRD06 1 50 
GR006 160 
GRD061 70 
GR0061 60 
GR006190 
GROO  620  0 
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DO  123  1*2.4 
DO  120  J~l,129 

IF  f MAP ALL ( J • I ) • NE • BLANK)  CO  TO  121 

120  CONTINUE 
GO  TO  123 

121  «RITE(8V  1000)  PLUS  .BLANK  . ( MAP  ALL (J.I  ) «J=1«129)  .BLANK 

123  CONTINUE 

124  CONTINUE 
GO  TO  100 

l 000  FORMAT  < 1 32A1  ) 

1001  FORMAT (2A1.129X.A1 ) 

ENO 


GR0062 10 
GRD06220 
GRD06230 
GR00624G 
GRD06250 
GR006260 
GRD06270 
GRD06280 
GRD06290 
GR006300 
GRD06310 
GRD06320 


IDENT  HMEAO 

SUBROUTINE  NHEAD 

COMMON/ B/  NLEVEL . I FORM ,NOO,IHIST#NTX%OPT(25) 

COMMON/ LIM/  V AL MA X • VAL  MIN# RANGE  ( 20 ) • VAL  <21  ) .PRANGEf 20 > #NFREQ(25) 
COMMON/STTS/  STAT(7) 

COMMON  /SYMBLS/  SYMBOL ( 25, 4 ) ,NS( 1 0 ) . T I TLE( 60 ) . TEX T ( 800 ) 

COMMON  /OUTPUT/  MAPALL ( 1 29 # 4 ) 

U I MENS  ION  MAPL  ( 84 , 4 ) , JSYM(  12,4  » , AS YM (12,41  ,FMT ( 5 > , MAPL A ( 1291 
TQUi VALENCE  ( MAPL ( 1#  1 ) . MAPALL ( 1 • 1 ) ) • ( J S YM ( 1,1)  « MAPALL (60 • 3)  ) 
EQUIVALENCE  (MAPL A ( 1 ) , MAPALL ( 1 • l ) ) • ( A SY  M ( 1 • 1 ) .MAPALL ( 1.4)) 
DIMENSION  LON ( 10) , IGH< 1 1 > 

INTEGER  BLANK # SYM I # SYMBOL • ASYM 
LOGICAL  FIN, OPT 
C INTEGER  Z 

C /=  1 1 

DATA  FMT/* <2X#  • # >2A4# • • • l OX, • # • 1 OF  1 • , • 0.2)  '/ 

DATA  A/*M INI • ✓ , B/ • MAX  I #/,C/*MUM  •/,□/•  •/,E/#I0X,*/ 

OATA  LOK/'L' , •□•,•••  ,•  # r * V • , • A • , *L # , #U • , •£ • , • S •/ 

DATA  IGH/ "H* « • I • 9 'G* « "H1 « " • , • V * , • A • , • L » , • U • , • E • , • S • / 

OATA  BLANK/1H  /, SYM I / INI / , EQS/ 1H=/ 

NW0RK=NLEVEL* 1 
FMT  4 31 =D 

IF(0PT(4>)  FMT  < 3 ) = E 
DO  120  1*1, 129 

1 20  MAPLA(  I ) = EL  ANK 

IF (OPT (22)  ) GO  TO  130 

• RITE  (8,  1000)  NLP VEL , VALMI N, VALMAX ,ST AT ( 4 ) • ST  AT ( S ) 

I F ( • NO  T • OP  T ( 1 2 ) ) GO  TO  500 

• RITE(8#20  10  ) 

• PITE(  8, 20  30  ) ( N , V AL ( N ) , VAL ( N4  1 ) , NFREQ ( N ) , N* I ■ NL  F VLL  Q 

RE  TURN 

50  0 «RI TE ( 8 , 1 ◦ 10  ) 

• RI  TE( 8. FMT ) A,C,(VAL(  I ) , I * i , NLEVEL  ) 

• RITE(8#FMT  ) 8 • C • (V  4f_  ( I ) # I=2,N»ORK  ) 

•RITE  (8,1020) 

• RITfc(8,FMT)  0,0,  ( PRAN<*E<  I ) • 1 3 1 • NLEVEL  ) 

•RITE  (8. 1030 ) 

1 JO  J = 0 

1M  .NOT.OPT  ( 4 ) ) GO  TO  ?00 


G RD0  6 329 
GRD06330 
GRD06340 
GRD06350 
GRD06360 
GRD06370 
GRD06380 
GRD06390 
GRD06400 
GRD06410 
GRD06420 
GRD06430 
GRD06440 
GRD06450 
GRD06460 
GRD064  70 
GRD064 80 
GRD06490 
GRD06500 
GRD06510 
GRD06S20 
GHD06530 
GR006540 
GROO  6550 
G ROO  656  0 
GRD06570 
GRD06580 
GROO  6 590 
G R DO  660  0 
GRD066 I 0 
GRD06620 
GROO  6630 
GRPO  0640 
GRD06650 
GRD06660 
G RDO  66 7 0 
GRDO660O 
GRD06690 
GROOO  700 


J=  J*  1 

DO  2C 1 K * l , 4 

JSYM( J #K ) = SYMBOL(  2 1 • K ) 


201  4SYM(J,K)  = SYMBOL ( 2 2 • K ) 

DO  210  L=I, NLEVEL 

LL  = NL  E VEL ♦ 1-L 
210  NFREGILL 41  )=NFREQ(LL  ) 

NFRE  Of  1 ) = NFREQ (21  ) 

DO  211  1=1,10 

2 1 l MAPL A(  I ) = LC • ( I ) 

20  *>  DU  2C2  JJ  = 1 #NLEVEL 
J = J ♦ 1 

00  202  K=l#4 

JSYM( J,K ) = 5 YM  60L ( JJ,K) 

202  ASYM( J,K )=SYMBOL( JJ410,K ) 
IF(  #NO  T *OP  T ( 5 ) ) GC  TO  203 
J = J4  1 

DO  204  K * 1 , 4 

JSYM ( J ,K ) = SYMBOL ( 23,K ) 

204  ASYM(J,K ) = SYMBOL ( 24 • K ) 

NF  RE  Q ( J ) = NF  REQ  ( 2 3 ) 


GROO 6710 
GRD06720 
GRD06730 
GRD06740 
GHD06750 
GROO  6 76  0 
GROO  6 7 70 
GROO  6 700 
GROO  6 7 90 
GRDO60U  0 
GRD068 1 u 
GRD06820 
GRD06830 
GRD06840 
GR006850 
GRD06860 
GRD068  70 
GRD06880 
GRD06890 
GRD06900 
GR00691 0 
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JJ*  < J-  1 )*  10 
00  212  1 3 1 • II 
212  MAPLAC I+JJ )* IGM( I ) 

203  LL* J 
JaJ-1 

F I N=  .FALSE. 

NBORMLL*  10*1 

WRITE (8.  1 1 0 0 ) ( M APL A ( I ) .1=1  .NWORK) 

WR  1 TE  ( 8 . 11  10H  I.I=1.J) 

122  WR I TEC 8*  1 100  ) (EOS  ,l3l, NwORK) 

IF(F1N>  GO  TO  127 

123  00  124  J* 1.2 

WRITE  (8*1120)  ( < JSYM( L* 1 > ,M=l , 9) ,L=1 ,LL) 

00  124  K = 1 • 3 
KK=5-K 

124  SR  1 TE  (8*1130)  ( ( JSYM(L  *KK ) *M* 1 9 9 ) *L*  1 *LL  ) 

I F ( F IN  ) GO  TO  122 

WRITE! 8. 1 140  ) ( ( J SYM( L * 1 )*M=1*4)*  ASYM(L»1 

Cl*  L*l,LL> 

DO  126  K 3 1 * 3 
KK=5-K 

126  SR1TE( 8. 1 130 ) ( ( J SYM  ( L * KK  ) * M3  1 * 4 ) * A SYM ( L * KK ) * 

Cl*  L31«LL) 

F IN*. TRUE* 

GO  TO  123 

127  WRITE  (8*1160)  (NFREQI  I )*1  = 1*LL) 

IF(  .NOT .OPT ( 2 1 ) )RE TURN 
BLANK  PRINT  ARRAY 
00  100  I3  I * 4 
00  100  J= 1 . 84 

100  M APL  ( J * 1 ) 3 BL  ANK 
INITIALIZE  PRINT  ARRAY 
J=  1 

MFRE  Q3  0 
DO  104  I31*LL 
MAPL(J*1>  3 SYM I 
MF REG3 MFREQ*NFREQ( I ) 

J*«J4  1 

OC  103  K=  1 ,4 
00  101  L=’»2 
MAPL ( J.K  >3JSVM( 1 ,K ) 

101  J=J* 1 

M APL ( J , K ) 3 ASY M ( I * K ) 

J=  J ♦ 1 

OC  102  L=  1 • 2 

MAPLU.K  ) 33SYM  ( I * K ) 

102  J = J ♦ 1 

103  J= J-5 
J = J ♦ 5 

MAPI.  ( J.  1 ) 3S  VM  I 

104  J=J+ l 
MXFREQ=0 
00  105  I 3 1 * LL 
6IGN6R  = NFREQ(  I )*  100 
BIGN01  3 MFREO 
JIGN0R  3 ( b ibNBR/fi  IGNO  I ) ♦ l 
IF(NFREQ( Il.GT.O)  NFREQi I ) =3 I6N0K 

105  MXF REQ  = MAX0 (NFREQI  I ) .MXFREQ ) 

MXFREQ=MXFREQ— 1 
KiQRK=7«LL 


GRD06920 
GRD06930 
GROO  6940 
GRD06950 
G ROO  696  0 
GR006970 
GRO06980 
GROO  6990 
GRDO7C00 
GR0070 10 
GROO  70  20 
G ROO  70  30 
GROO  7040 
GROO  70  50 
GROO  7060 
GRD07070 
GROO  7080 

( JSY  M( L • 1 ) *M=1 .4GR007090 

GRD07 100 
GROO  7 1 1 0 
GROO  7 1 20 

( JSY  M(L*KK)  *M31  * 4 GROO  7 1 3G 
GROO  7 140 
GROO  71 50 
GROO  7 16  0 
GROO  7 1 70 
GROO  71 80 
GR007 1 90 
GRD07200 
G ROO  7 2 l 0 
GRD07220 
GROO  7230 
GROO  7240 
GROO  7250 
GRDO7260 
GROO  72  70 
GROO  7 280 
GROO  7 290 
GROO  7300 
GROO  7 3 10 
GROO  7 320 
GROO  7 330 
GRDO  7340 
GROO  7 350 
G ROO  7 360 
GRDO  7370 
GRDO  7 3 80 
GRDO  7 390 
GROO  74  00 
GRDO  7 4 1 0 
GRDO  7420 
GROO  74  30 
GROO  7440 
GRDO  744 1 
GRDO  7442 
GRDO  74  4 3 
GROO  7450 
GROO  7460 
GROO  74  7 0 
GROO  7480 


u u u u 


c 

c 


DO  110  I = 1 i MXFREQ 

GROO  74  9C 
GRDO /600 

BLANK  USED  BARS 

GRDO  7b 1 0 

DO  107  J=1.LL 

GRDO  7520 

IF  ( NFREQ ! J ) .NE.  I- 1 ) GG  TO  107 

GRDO  7530 

KSTART=7*J-6 

GROO  7540 

KSTOP=  74 J 

GRDO  7550 

DO  106  K^KSTART, KSTOP 

GRDO  7 5o0 

DO  106  L— 1.4 

GRDO  7570 

106 

MAPL  <K.L)=9LANK 

GROO  7580 

107 

CONTINUE 

GRDO  759C 
GROO  7600 

•RITE  LINE 

GR00761 0 

WRITE  (8,  1220)  I , ( MAPL  !K  , 1 ) , K= 1 , K WOR K ) 

GRDO  7620 
GRDO  7630 

WRITE  OVERPRINT  LINES 

GRDO  7640 

DO  109  3=2.4 

GRDO  7650 

DO  108  K>lyKWORK 

GRDO  7660 

IF  ! MAPL  (K.J)  • N E • BLANK)  GO  TO  109 

GRDO  7670 

108 

CONT INUE 

GRDO  7680 

GO  TO  110 

GROO  7690 

109  WRITE  ( 8 • 1200  1 (MARL  < K . J ) • K= I . KwORK ) 

110  CONTINUE 

WRITE!  0,  1230  ) NFRE0I25) 

RETURN 

1000  FORMAT  ! //  / • DATA  MAPPED  IN  • • 1 3 . • LEVELS  BETWEEN  EXTREME  VALUES 
C.F10.2.*  AND • .FIQ.2.*  MEAN  a'fF10.2, ' ST.  DEV.  »*fF10.2l 
1010  FORMAT  1///44H  ABSOLUTE  VALUE  RANGE  APPLYING  TO  EACH  LEVEL) 
1020  FORMAT  1///64H  PERCENTAGE  OF  TOTAL  ABSOLUTE  VALUE  RANGE 
10  EACH  LEVEL) 


GRDO  7 700 
GRD0771 0 
GRDO  7 7 20 
GRDO  7 730 
OF • GRDO 7 740 
GRDO  7 750 
GRDO  7760 
APPLYING  T GRDO  7 7 7 0 
G ROO 7 78  0 


1030  FORMAT  C///58N  FREQUENCY  DISTRIBUTION  OF  DATA  POINT  VALUES  IN 
l LEVEL  ) 

1100  FORMAT! 1 IX. 121 A1 ) 

1110  FORMAT!*  LEVELS* ,9X9 *0* 9 1 1!8X, 12)  ) 

1120  FORMAT  ! 1 1X9  1 2!  IX, 9A 1 ) ) 

1130  FORMAT  ! 1M+.  1 0 X . 1 2 ! 1 X . 9A  1 ) f 
1140  FORMAT!*  SYMBOLS  *•  12!1X.9A1)) 

1 1 6 C FORMAT!*  FR E QUENCY *,  I 7 , l l I 1 0 ) 

1 2 OC  FORMAT  !IH*9  9X9  l 2! 3X,  7A  1 ) ) 

1220  FORMAT  ! 4X.  16.12!  3X.7A1  ) ) 

123)  FORMAT ( 10 X. • HI STOGRAM  EXPRESSED  AS  PERCENT  OF  ALL  NON-BACKGROUND 
CELLS*/ 10X9 • NUMBER  OF  BACKGROUND  CELLS  =*,I5) 

20  10  FORMAT  !6H  LE VEL , 9X , 7MM I N I MUM,  eX , 7HMAX  I MUM  , 8X , 9HFRE QUENCY  ,/  » 

PO  JO  FORMAT IIS.7XF10.2.5XF10.2.9XIS) 

END 


EACHGRDO  7790 
GRDO  7 80  0 
GRDO  781 0 
GRDO  7820 
GRDO  7830 
GRDO  7 84  0 
GRDO  7850 
GRDO  7860 
GROO  7©  70 
GRDO  788  0 
CGRDO  7890 
GRDO  7900 
GRDO  7910 
GRDO  7920 
GRDO  7930 
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t 
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I Dt  NT  FLE  X IN  GHD07939 


SUBROUT  INF.  FLE  X I N ( I FORM,  T,F  IRST  ) 

GROO  7 940 

■ 

DIMENSION  SO I L { 22) 

GRD07950 

OATA  SO  IL/4  « i 5 • >6*  » 25*  • 2 £«  *59*  • 6 0 • «63#  *64* 

.76. *1 13..1 14. . 1 33. . 134 .GRD07960 

1 . 1 38. • 168*  « 170  •• 20  6. *207. . 2 08.  | 2 I 7 • » 7 1 • / 

GRDO  7970 

NTYPE =22 

GROO  7900 

READ! 11.11)  S. V 

GROO  7990 

1 1 

FORMAT (3X.F3.0.30X.F2.0) 

GRD08000 

T=  10  . 

GRDO0O 10 

DO  300  J=t. NTYPE 

GR008020 

IF < S.EQ.SO IL ( J ) ) GO  TO  450 

GRD080 30 

300 

CONTINUE 

GR008040 

I 

IF<  V.EO. 12.)  T = 1 . 

GROG  8050 

V 

IFlV.EQ.d.)  T=2. 

GRD08060 

IF< V.EQ. 11.)  T =3. 

GR008070 

I F ( V .EG • 3 • ) T = 4. 

GROO  80  8 0 

IFCV.EQ.3l. ) T*4. 

GRD08090 

GO  TO  500 

GRD08100 

450 

(FCV.NE.12.)  T=5. 

GRD081 1 C 

i 

IF(V.NE.8.)  T * 5 • 

GRD08120 

IF  C V.NE  .11.)  T = S . 

GROO  8 1 30 

I F C V «NE • 3 . ) T * 5 • 

GR0C81 40 

I F ( V .NE • 3 1 • ) T = 5 • 

GRD081 50 

lMv.EQ.12.)  T = 6 . 

GROC  8 160 

IMV.EQ.8.)  T=7. 

GRD08I 70 

IF C V.EQ. 11.)  T 2 8 . 

GRD08180 

1 

IFCV.EQ.3.)  T=<5. 

GROOQ 190 

1FIV.EQ.31. ) T*9. 

GR008200 

500 

CONT INUE 

GRD082 I 0 

> 

RETURN 

GRQO022O 

• 

ENO 

G R008230 

84 
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APPENDIX  B 

VANDENBERG  IRREGULAR  OUTLINE 
FOR  GRID 


AD-A040  462  SAN  DIE60  STATE  UNIV  CALIF  CENTER  FOR  REGIONAL  ENVIR— ETC  F/G  6/6 

EC0L06ICAL  ASSESSMENT  OF  VANDENBERG  AIR  FORCE  BASE.  CALIFORNIA.— ETC(U) 
SEP  76  R M REILLY.  F P STUTZ.  C F COOPER  F29601-75-C-0116 

UNCLASSIFIED  AFCEC-TR-76-15-V0L-3  NL 

^0*0462  K32H 

b ■■HliillllilllllH 

■ END 

DATE 

FILMED 

7-77 


T 


Column 

11111111 


12345678901234567 

1 

15 

13 

1 

15 

12 

2 

16 

12 

1 

16 

13 

1 

17 

13 

1 

17 

14 

1 

18 

14 

2 

18 

15 

2 

18 

17 

1 

17 

18 

2 

17 

19 

1 

17 

24 

3 

16 

24 

1 

15 

34 

2 

15 

35 

1 

14 

35 

1 

14 

34 

1 

14 

33 

1 

13 

32 

2 

13 

30 

1 

12 

31 

1 

12 

32 

1 

11 

33 

1 

11 

34 

1 

11 

32 

1 

10 

30 

1 

10 

28 

2 

9 

28 

3 

8 

26 

2 

7 

26 

2 

6 

26 

4 

5 

26 

1 

5 

27 

1 

5 

28 

1 

4 

29 

1 

4 

32 

2 

4 

34 

1 

2 

33 

1 

4 

32 

1 

5 

30 

1 

5 

23 

1 

5 

21 

1 

5 

20 

1 

6 

16 

1 

8 

16 

87 
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APPENDIX  C 

SUBROUTINE  FLEXIN  FOR  GRID  OUTPUT  DISPLAYING 
SOILS  WITH  HIGH  EROSION  POTENTIAL 


SUBROUTINE  FLE XIN < IFORM, T ,F I RST» 

REAL  S0ILU83I, LEVEL!  183)  _ ,,  ,,  „ 

OAT  A SOIL/2* , 3 . » 9 • , 16  • , 18  » ♦ 23«,29.,25. , 2 7. ,28 . « 29 • , 30»»31.» i 33» 
l«39.»35.,  37.,3  8.,93.»55.»56.»57.|58.»59.»60.,62.»66.«67.»7.l.,79.»7 
2 5*.  76 *.78.,  79., 80. .61. ,91. .101*. 107. ,109*. 1 1 2.. 1 13.  ,114.,  115 .,  1 18. 

7..  261.. 263. .270.. 2 72. .2 73.. 2 75. .277. .278.. 281. .282. .283.. 289. .2 86. 

8. 28 7. . 288.  ,290., 2 91. ,292.  ,2  18  . , 2 2 7 . , 5 . . 99.  ,96  . ,98  . . 52  . , 63  . , 69  . , 68  . 

9.89.. 52..93.  , 102.  ,103.  ,106.  .LOB.  , 110.  , t 30  .,  170. ,172.,  177..  .79.,  186 

11..  187.. 190. .2 19. .220. .229.. 230. .239. .290. .291. .251. .252. . 253 .. 255. 
12, 256.  ,269.  ,265.  , 267.  , 2 7 l • « 2 /9»,276.r,?79..280.,285.,209»«6.,65.,  Ifa9 

13..  188.,99.,99.,ll6.  , 1 59.  , 160.  , 1 89  . , 09  . ,207»,208»,257.«L,3.,  1 79  • , 9 

199.. C./  „ , 

DATA  LEVEL/ 11 7*1. ,2*2. ,96*3. ,9*9. , 19*5. / 

NT  'PE  = 183 

READ! 9,101  T 
10  FORMAT! 3X.F3.0} 

DO  3C0  J=  l.NTYPE 
IF( T.EQ.S0IL1 J1J  GO  TO  950 
300  CONTINUE 
T=  0 

GO  TO  960 
950  T*LEVEL1J» 

960  CONTINUE 
RE  r URN 
ENO 


89 


APPENDIX  D 

SUBROUTINE  FLEXIN  FOR  GRID  PROGRAM  DISPLAYING  AREAS  OF 
HIGH  EROSION  POTENTIAL  BASED  ON  SOILS  AND  VEGETATION 


SUriKOLTINt  FLEXlNUFOR^.T, FIRST) 

REAL  SC  I L ( 183)  , LEVEL!  133) 

DATA  Suit  / 2 . ,3.,4.,l6.«l3.,23.,24.,25.,27.,2c.,29.,  30., 31.,  i 2 • , ^ • 
l,  34. ,35.,  27«,2i:.,43.,53.  ,5c.  » j / • » 3d  • » ‘j'?  » tit?»i71»i7^»t  7 

25..  16..73. , 79. ,60. ,81. ,91 . ,101.  , 107. , 109. , l 12. , llj.  , l 14. , 1 13  .,  I U. 

3. 115. . 123. .126. .127. .129.  ,131. ,132. , 133. ,13 4. ,136. ,137. ,133., It  1.. 
4 16  2 .,  166.,  167.  , 168. , 1 75.  , 176.  ,135. , 196. ,20^. , 2 0 6 . ,2  1 1.,  212. , 2 1 3 ., 2 

514. .  2 15. ,216. ,417. ,226. ,223. ,231. ,232.,  233.  ,45,. , 2 3 j . , 2 1 6 • , 257 . , 2 5 

68.,  242.  ,243., 244.  ,245.  ,246.  ,2*,  7.  ,243.,249.,  230.,2*>‘,.,23o.,29  5«,  260 

7.. 261..263..27u.,272.  ,273.  ,275.,  277.,  27  3.  ,251.,  232.  , 253. ,2c*,.,  7 ij. 
8, 267.  ,263. ,290. ,291. ,292.  ,2li.  ,227.  , 5 . , 44  . , 4o  . , 43  . , *J  2 . , 5 « , 6 4 « , o . 
9,39.  ,52. ,93. , IC2.  ,103.  ,106.  ,103.  ,110.,  130**1  7^., 174. ,177.,  I 7 9 • , I >n 

11. , 167. , 150. ,215. ,220.  ,229. ,230.  ,239. ,2^0., 241. ,251.  ,232. ,233  .,  > . 

1 2 , 2 5 6 . , 2 6 4 . , 2 6 5 » , 2 6 7 . , 2 7 1 » , 2 7 4 • , 2 7 6 • , 2 7 ) • , 2 8 0 . , 2 8 5 . , 2 39  • , c*  • , 6 5 . » 1 7-  9 

13.  , 168.  ,54.  ,99.,  116.,  I 5 9 • ,160.  ,184.,204.  ,207«,2o6.,257«,  173.,  17*.. ,4 

149.,  C./ 

DATA  LEVEL/117*!. ,2*2.  ,46  ♦3.  ,*»*4.  , 14*5./ 

NTVPE=183 
Rl AU ( 5, l l 1 S,V 
11  FORMAT! 3X,E 3.0, 30X.F2.  01 
T = S 

DO  300  J=1,NTVPl 
IF!  T.EU.SjIH  Jl)  GO  TC  45C 
300  Cun t l nue 
T = 0 

GO  TO  4oC 
450  T=Lfc VEL 1 J) 

46  0 CUN  T I NUE 

IF!  V.E0.16.)  f-0  TC  461 


46  1 


IF ( V.NE . 16. 1 GO  TO  470 

IF!  T . r 


> Eu . I . ) T = 6 
IF1T.EC.2.)  T=  7. 
IF! T.L0.3.)  T=6. 
IF!  T.EJ.4.I  T =5. 
IF ! T .F  J .5)  T= 1C. 
470  CUN  T I NOt 
F E TORN 
END 
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RETURN 

END 
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APPENDIX  F 


SUBROUTINE  FLEXIN  FOR  GRID  MAP  DISPLAYING 
AREAS  OF  SUITABLE  HABITAT  FOR  THE 
CALIFORNIA  LEGLESS  LIZARD  ON  VANDENBERG  AFB 


SUBROLTINE  F LE X ! N(  I FORM » T ,F  I R $ T J 

REAL  SOIK22I  . . , „ . 

DATA  SO  lL/4.t5.*6.*25«*28.*59.*&0.*65.*64.*76.*  i 1.  3.*114.*1j3.*134. 
1*138. *168. *170. *206. *207.  *408*  *217.  *71./ 

NTYPE=22 
REA0(9,IU  S,V 
11  FORMAT!  3X*F3.0.30X,F2. 71 
T*  10. 

00  300  J = l , N T V P E 

IF ( S.EO.SOlLi J >1  GO  TO  450 
300  CONTINUE 

I F ( V.E0.12. ) T= 1 . 

IFiv.FU.8.)  T = 2. 

IFlV.EQ.ll.)  T = 3. 

1 F ( V.EU.3.)  T =4. 

IF1V.E0.31.)  T =4 . 

GO  TO  500 

450  IFIV.NE.12.)  T=5. 

1 F l V . Nf  . 8.  ) T = 5.. 

IFIV.Nfc.ll. I T = 5. 

V IFIV.NE.3.)  T = 5 . 

* IFIV.NE.3l.)  T =5 . 

IF< V.fcy.12.)  T ~6. 

IFIV.fcO.8.)  T = 7 . 

IFiv.FO.ll.)  Tr e. 

IF( V.FQ.3.)  T=9. 

I FI  V. EO  .3  l. I T=9. 

500  CUN  T INUfc 
PF  TORN 
END 
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APPENDIX  G 

SEARCH/COUNT  PROGRAM  AS  ADAPTED  TO  THE  BURROUGHS  3500  COMPUTER 

l COMPILE  SNFT10  XFORTN  LIB  DATA  CARDS 


1 LSTl 
I DENT  MAIN 

PILE  U=AnFTOT,UNIT=TAPE*RECORO=SO*BLOCKIn6=1. UNLABELED* FIXED 
C FREQUENCY  COUNTING  PROGRAM 

c old  runray  and  ner  construction  surrounding  it  (43  cells  ) 

REAL*®  C(43>’CX 
integer  S,EX*EL*D.IO 
INTEGER  V<9> 

integer  valuesuoo>*valuexuoo)*valuel<ioo)  * valuev <5o » 

REAL  FREQSdOO)  * FREQX ( 100 > *FREQLC  100)  »^REQV ( 50 ) 

REAL  ACS(100)*ACX(100)*ACEL(100)*aCV (100) 

DATA  C/*  A0O86* • • BB086* . • CB086* » *2X085' * ’ *8005' * • B0O85' . 'CgOsS' » • 00 
1085'  * 'XAOb**'  , • YA081*'  * 'ZA084* , * A8O04*  , *88084* , * C30B1*'  * »DBO04f , 1 YAO0 
23*  * • 2*083* * » AB083*  * *BB083* * *CB083* * *05083*  * *EB083* * • ZA082*  * • AB0B2* 
3* *80082** *C0O82' * *08082' * ' EB082* * ♦ AB08I • # * 8B081 * * * CB081 * * *08081'** 
4Eb081*»*Bb080**'C8080'»  *08080'*  *E3080' * 'CB079' * *00079* * *£3079' * *FB 
5079* . 'DB078*  * *£8078*/ 

NN=43 

JFREQS=0 
JFREQXSO 
JFREQL=0 
JFREQVSO 
DO  20  KS1.NN 
20  VALUES(K) rO 
00  30  Ksl,NN 
30  VALU£X(K)sO 
DO  40  Ksl*NN 
40  VALUELUOsO 
DO  50  Ksl,NN 
50  VALUEVOOrO 
DO  60  Ksl » NN 
60  FREQS(K>=0 
DO  70  KslrNN 
70  FRE8X ( K) SO 
DO  80  X:l*NN 
80  FREQL(K)=0 
DO  90  Ksl.NN 
90  FREQV(K):o 
DO  270  K*1*NN 
270  ACS(K>ro.o 
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DO  271  K=1*NN 

271  ACX(K)S0,0 

DO  272  K=1»NN 

272  ACcLOOsO.O 
DO  273  K=1*NN 

273  ACV(K)=0t0 

8 *EAO(U'll)  S*EX;EL»(VCI)#I  = 1»9**D»CX*1D 

11  F0RMAT<3X,I3,4X«I2»4X»I1»1X.9C2X,I2>,UX,I«»,1X,A5,1X,I4) 

DO  300  JSl'NN  “ 

IF(CX.EQ.CU))  50  TO  450 
300  CONTINUE 

IF < 1D.LT . 4646>  50  TO  8 
50  TO  3005 

450  NRlTE(6»l2>  S»EX»EL#(V(I)*Is1*9)»0*CX*10 

12  FORMAT l *0*r *S=,*l3»iX»'EX='»I2*lX, •EL=*»H*lX» • V=* * 9< 12  * 2X ) . 9X » 14, 
UX’AS'IX'I*) 

IFUFREQS.E8.0)  So  TO  110 
DO  100  U=l, JFRE9S 
IF ( S*NE« VALUES (U) ) SO  TO  100 
FREQS<J)sfRE&S(J)*1 
GO  TO  200 
100  CONTINUE 
110  JFREOSsJFREaSU 
VALUES! JFREQS I -S 
FREQS(JFREQS)=1 

200  continue 

IF(JFREQX.EQ.O)  So  TO  610 
DO  600  J=1»JFREQX 
IF!EX.NE.VALUEXU)  1 GO  TO  600 
FREax(J)SFREQX(J)«>l 
60  TO  800 
600  CONTINUE 
610  UFREOXsjFREaxn 

VALUEX ( JFREOX ) sEX 
FREQX( JFR£0X)=1 
800  CONTINUE 

iF(JFREOL.Ea.O)  60  TO  1050 
DO  1000  J=1,JFRE8L 
IF(Eu.NE«VALuEL(U>)  So  to  1000 
FRE8l<J)=FREqL(J)*1 
50  TO  1200 

iooo  continue 
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1050  JFREaL=.JFREQl_n 

VAUIEl(JFRE3L>=EL 
FREQL(  JFREQDS1 
1200  CONTINUE 

00  3000  K=1'9 

iFlJFREaV.Ea.O}  60  TO  2050 
00  2000  Jr  1 > jFREQv 
IF(V(K) ,NE.VALUEV( J)>  SO  TO  2000 
FR£8V< J)=FREQV( J)*l 
SO  TO  2200 
2C00  CONTINUE 
2u50  JFRE3yrjFREQV*i 

VALUEV (JFREQy)sV(K) 

FRE8V(JFREaV)=l 
2200  CONTINUE 
3000  CONTINUE 

IF(ID.LT.4646>  SO  TO  8 
3005  CONTINUE 

00  250  L=1»Nn 

250  ACS(L)=FREfiS<L)*22.957 
DO  251  L=1»NN 

251  ACX(l)=FREQX(L)*22,957 
DO  252  L=1»Nn 

252  ACEL<L)=FRE9L<L)*22.957 
DO  253  L=1»NN 

253  ACV«l)=FRE3V(L)*2.55 
*RIT£<6»3010) 

3010  FORMAT1///1X, ’SOIL  TYPE  FRE3.  EST»  ACRES’) 

WRITE<6.3015>  f VALUES(N) »FRESS(N) , ACS<N) # N=1»NN> 

3015  FORMAT ( 3X» I 3»7X*F4.0»5X*F5.0) 

«RIT£<6,3020> 

3020  F0«MAT(///1X, ’EXPOSURE  Frcq,  EST,  ACRES’) 

tfRlT£(6,3025)  < VALUEX<N> *FRESX(N) , ACx<N) * Nri»NN) 

3025  FORMAT(3X,I2#8X*F4.0*5X»F5.0) 

VRlTfc  (6*3030 ) 

3030  F0RMAT(///1X, ’ELEVATION  FR£Q.  EST*  ACRES’) 

RRIT£(6<  3035)  (VALUEL(N) *FREQL(N) ,AC£L(N) » N=1»NN> 

3035  FORMAT <3X, I 1,9X»F4.0»5X*F5.0> 

MR  1 T£( 6,3040) 

Toro  format ( ///ix *’ vegetation  fres. (sjbcells)  est.  acres*) 
WRl T£ ( 6*  3045)  ( VALUEV ( N) * FREQV ( N> , AC V ( N) . N=1»NN> 

3045  FORMAT (3X» 12, 10XpF5.0»11X*F5.0) 

STOP 

ENO 


APPENDIX  H 


SMALL  CELL  COUNT  PROGRAM  LISTING  AS  ADAPTED  TO  THE  BURROUGHS  3500  COMPUTER 

1 

2C0MPILE  SNFT10  XFORTN  LIB  DATA  CARDS 
3 
11 

3LST1 

8 

IDENT  MAIN 

FILE  11=ANFT0T,UNIT=TAPE,REC0RD=80, BLOCKINGS ,UNLABELED , FIXED 
C SMALL  CELL  COUNT  PROGRAM 

INTEGER  V(9) , VALUEV(35),  FREQV(35) 

JFREQV=0 

8 READ( 1 1 ,11)  (V(I),I=1 ,9),  ID 
11  FORMAT (18X,9(2X, 12) ,22X,I4) 

DO  3000  K=1 ,9 

IF(JFREQV.EQ.O)  GO  TO  2050 
DO  2000  J=1 ,JFREQV 
I F ( V ( K ) .NE. VALUEV (J ) ) GO  TO  2000 
FREQV(J)=FREQV(J)  + 1 
GO  TO  2200 
2000  CONTINUE 
2050  JFREQV  = JFREQV  + 1 
VALUEV( JFREQV)  = V(K) 

FREQV(JFREQV)  = 1 
2200  CONTINUE 
3000  CONTINUE 

IF( ID. LT. 4646)  GO  TO  8 
WRITE(6,300) 

300  F0RMAT(//1X, 'VEGETATION  FREQ.') 

WRITE(6 ,400)  (VALUEV(N) ,FREQV(N) , N=l,35) 

400  FORMAT (5X , I2,7X , 14) 

STOP 

END 


2END 

3 
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(The  reverse  of  this  page  is  blank) 


Hq  USAF/PREE 

Hq  USAF/PREVP 

Hq  USAF/PREVX 

Hq  USAF/RDPS 

Hq  USAF/SAFOI 

Hq  USAF/SGPA 

OSAF/SAFIL 

AFSC/DE 

AFSC/DEV 

AFSC/SGB 

AFSC/DLCAM 

ATC/DEPV 

CINCSAC/DEV 

CINCSAC/DEPA 

CINCSAC/DEPV 

CINCSAC/SGPA 

TAC/DE 

AFRES/DEEE 

\ AFIT/DEM 

‘ AUL 

AFOSR/NL 

AMRL/DAL 

OEHL/CC 

, OEHL/OL-AA 

OEHL/OL-AB 

AFWL/SUL  (Tech  Library) 
USAFSAM/EDEO 
USAFSAM/EPS 
AFRPL/Library 

* AMRL/THE 

SAMTEC/SEH 

> SAMSO/DE 

SAMSO/DEP 

SAMSO/LVRO 


INITIAL  DISTRIBUTION 

1 Defense  Res  & Engr/AD  (E&LS) 

2 Dir,  USA  WW  Exp  Sta 
2 Chief  of  Naval  Ops 

1  Environ  Protection  Div,  OP-45 

1 Technology  Transfer  Staff  (EPA) 

2 National  Science  Foundation 

2  U.S.  Fish  $ Wildlife  Services 

1 Tetra  Tech,  Inc 

2 McDonnel-Douglas  Aircraft  Co 

1 Aerospace  Corp 

2 Tech  Applications  Ctr 

1 San  Diego  State  Univ 

3 San  Diego  State  Univ/ Library 

2 DMA/Aerospace  Center  (PPRN) 

1 ISTRAD/CC 

1 4392  Aerosg/CC 
1 4392  Aerosg/DE 
1 4392  Aerosg/ACD 
1 Scripps  Inst  of  Oceanography 
1 Library 

1 Univ  of  Calif,  Geography  Dept 

1 AFRCE/Eastern  Region 

3 AFRCE/Central  Region 

2 AFRCE/Western  Region 

2 S.  Central  Coast  Commission 
1 Environmental  Analysis  Systems 

1 U.S.  Forest  Service 

2 Los  Padres  National  Forest 
1 Cal  State  Polytechnic  Univ 
1 Dept  of  Landscape  Arch 

1 Univ  of  Calif,  Library 

3 Calif  Polytechnic  State  Univ 
1 Library 

1 Ctr  for  Regional  Environmental 


SAMSO/SGX 

SAMSO/WE 

SAMSO/PP 

SAMSO/JA 

SAMSO/OI 

SAMSO/RS 

SAMSO/DECP 

AMD/RD 

ADTC/CSV 

ADTC/DLOSL 

USAF  Regn  Civ  Engr 

1 Med  Service  Wg/SGB 

DDC/TCA 


1 Studies  2 

1 San  Diego  State  Univ/Biology  Dept  6 
1 Univ  of  Calif  Natural  Land  $ 

1 Water  Reserves  System  1 

1 Ctr  for  Regional  Environmental 

1 Studies/San  Diego  State  Univ  1 

3 San  Diego  State  Univ/Computing 

2 Center  1 

1 San  Diego  State  Univ/ Geography 

1 Dept  3 

2 Harvard  Univ/ Lab  for  Computer 

1 Graphics  5 Spatial  Analy  1 

12  Univ  of  Cal/Water  Resources  Ctr  1 
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INITIAL  DISTRIBUTION  (CONCLUDED) 


Stanford  Univ/Dept  of  Civil  Engr  1 
Santa  Barbara  Museum  of  Natural 

History  1 
Naval  Weapons  Ctr  1 
AFCEC/SU  1 
AFCEC/SUL  1 
AFCEC/WE  1 
AFCEC/EV  9 


